Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,589,125 T→C 5.8% F511L (TTT→CTT)  acrD → aminoglycoside/multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,589,1250TC 5.8% 384.4 / 10.0 121F511L (TTT→CTT) acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base T (60/53);  new base C (4/3);  total (64/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

Only 100 of 126 total aligned reads displayed.

CACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGCCG  >  NC_000913/2588977‑2589271
                                                                                                                                                     |                                                                                                                                                  
caccacCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCt                                                                                                                                                    >  3:442638/1‑149 (MQ=255)
   caccGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCtttt                                                                                                                                                 >  4:180936/1‑149 (MQ=255)
    accGGTGCCATCTATCGCCAGTTCTCAATTACCATTGTCGCGGCGATGGTGCTGTCAGTACTGGTATCGATGATCCTCACTCCGGCTCTGTGTGCCACTCTACTTAAGACACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCttttt                                                                                                                                                >  3:801179/1‑149 (MQ=255)
        gTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGc                                                                                                                                             <  1:398995/148‑1 (MQ=255)
         tGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCt                                                                                                                                           <  3:547584/149‑1 (MQ=255)
          gCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCt                                                                                                                                           <  1:186632/148‑1 (MQ=255)
               cTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGATTTTTTGCCTGGttt                                                                                                                                      >  2:372880/1‑148 (MQ=255)
                tATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTaa                                                                                                                                    <  4:853018/149‑1 (MQ=255)
                 aTCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTaa                                                                                                                                    <  2:471925/148‑1 (MQ=255)
                 aTCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                                                                                                   <  3:477543/149‑1 (MQ=255)
                   cGCCAGCTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCTACATGGGCAAGAAGGCTTTTTTGCCTGGTTTAAcc                                                                                                                                  >  1:326067/1‑148 (MQ=255)
                    gCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAg                                                                                                                                <  3:819468/149‑1 (MQ=255)
                     ccAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAg                                                                                                                                <  2:941172/148‑1 (MQ=255)
                     ccAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAg                                                                                                                                >  1:877851/1‑148 (MQ=255)
                          tctcTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCTGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGtt                                                                                                                           >  1:132296/1‑148 (MQ=255)
                          tctcTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGtt                                                                                                                           >  1:279132/1‑148 (MQ=255)
                            tctATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTa                                                                                                                         <  1:351376/148‑1 (MQ=255)
                            tctATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTa                                                                                                                         <  2:132296/148‑1 (MQ=255)
                             ctATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTaa                                                                                                                        <  2:242036/148‑1 (MQ=255)
                             ctATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAAc                                                                                                                       <  3:227275/149‑1 (MQ=255)
                              tATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCTACTGAAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAAc                                                                                                                       >  3:589540/1‑149 (MQ=255)
                              tATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAAc                                                                                                                       >  2:834989/1‑148 (MQ=255)
                               aTTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCg                                                                                                                     >  4:524451/1‑149 (MQ=255)
                                ttACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCttt                                                                                                                                                  <  3:680923/119‑1 (MQ=255)
                                ttACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCttt                                                                                                                                                  >  4:680923/1‑119 (MQ=255)
                                ttACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCg                                                                                                                     >  2:341829/1‑148 (MQ=255)
                                    cATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGt                                                                                                                                        >  2:37105/1‑125 (MQ=255)
                                    cATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGt                                                                                                                                        <  1:37105/125‑1 (MQ=255)
                                    cATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAAc                                                                                                                 >  1:413330/1‑148 (MQ=255)
                                     aTTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACg                                                                                                                <  1:834989/148‑1 (MQ=255)
                                      ttgttgCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGACTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGcc                                                                                                              >  3:189376/1‑149 (MQ=255)
                                            cggcgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACg                                                                                                         <  2:712914/148‑1 (MQ=255)
                                             ggcgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCt                                                                                                       >  3:175653/1‑149 (MQ=255)
                                              gcgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCt                                                                                                       <  2:94142/148‑1 (MQ=255)
                                                  tGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGa                                                                                                   >  1:57463/1‑148 (MQ=255)
                                                     tGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaaaa                                                                                                <  2:965516/148‑1 (MQ=255)
                                                      gCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaaaaa                                                                                               <  2:645586/148‑1 (MQ=255)
                                                        tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgg                                                                                             <  2:413330/148‑1 (MQ=255)
                                                        tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgg                                                                                             >  1:509011/1‑148 (MQ=255)
                                                          tCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgggg                                                                                           >  2:14635/1‑148 (MQ=255)
                                                            aGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTg                                                                                         <  2:877851/148‑1 (MQ=255)
                                                                    tAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAt                                                                                 >  2:68596/1‑148 (MQ=255)
                                                                      gCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtct                                                                              >  4:578644/1‑149 (MQ=255)
                                                                       cGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATtct                                                                              <  1:213507/148‑1 (MQ=255)
                                                                        gatgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTcc                                                                            >  3:12020/1‑149 (MQ=255)
                                                                        gatgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTcc                                                                            >  3:187367/1‑149 (MQ=255)
                                                                         atgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTcc                                                                            >  1:295093/1‑148 (MQ=255)
                                                                          tgatCCTCACTCCGGCTCTGTGTGCCACACTACTCAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCAc                                                                          >  4:698036/1‑149 (MQ=255)
                                                                             tCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGa                                                                                                            <  1:924522/112‑1 (MQ=255)
                                                                               cTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAg                                                                     >  4:301059/1‑149 (MQ=255)
                                                                                tCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAg                                                                     <  1:20323/148‑1 (MQ=255)
                                                                                 cACTCCGGCTCTGTGTGCCACACTACTTAAGCCAATG‑AAAAAAGGTGAGCATCATAGGCAAAAAAGCCTTTTTGcccatg                                                                                                                                        <  1:668043/80‑5 (MQ=37)
                                                                                  aCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTaa                                                                                                                        <  4:638805/95‑1 (MQ=255)
                                                                                  aCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTaa                                                                                                                        >  3:638805/1‑95 (MQ=255)
                                                                                   cTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTg                                                                 <  4:12020/149‑1 (MQ=255)
                                                                                      cGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcg                                                               <  2:806600/148‑1 (MQ=255)
                                                                                      cGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcg                                                               <  1:154491/148‑1 (MQ=255)
                                                                                       ggCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt                                                             >  4:582627/1‑149 (MQ=255)
                                                                                        gCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt                                                             >  1:808579/1‑148 (MQ=255)
                                                                                            tgtgtgCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGAtt                                                        <  3:131529/149‑1 (MQ=255)
                                                                                                tgCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATtgtg                                                     <  2:860019/148‑1 (MQ=255)
                                                                                                  ccACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAt                                                   <  2:660222/148‑1 (MQ=255)
                                                                                                   cacaCTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAtt                                                  >  2:244630/1‑148 (MQ=255)
                                                                                                    acacTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAATAGGCTTTTTTGCCTGGTTTAACCAGCTGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTATCCTGCGCTGGATTGTGATTTa                                                <  4:801179/149‑1 (MQ=255)
                                                                                                    acacTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGAGAATTTTCGCCAc                                                                          >  3:741713/1‑123 (MQ=255)
                                                                                                     cacTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTAt                                               <  3:391605/149‑1 (MQ=255)
                                                                                                     cacTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGAAGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTAt                                               <  4:518266/149‑1 (MQ=255)
                                                                                                      cctactTAAGCCACTG‑AAAAAAGTTCAGCATCAGGGGCAAAAAGGC‑TTTTTGCCTAGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGt                                             <  3:102181/148‑1 (MQ=255)
                                                                                                      actactTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTAt                                               >  1:819947/1‑148 (MQ=255)
                                                                                                      actactTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGACGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTAt                                               >  2:630194/1‑148 (MQ=255)
                                                                                                        tactTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAggg                                                                                            <  4:256860/101‑1 (MQ=255)
                                                                                                        tactTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAggg                                                                                            >  3:256860/1‑101 (MQ=255)
                                                                                                          ctTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCct                                          >  4:187236/1‑149 (MQ=255)
                                                                                                            tAAGCCACTG‑AAAAAAGGTGATCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCACCGCCGCACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCctg                                         <  2:21659/148‑1 (MQ=255)
                                                                                                              aGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCctgctg                                      <  4:305572/149‑1 (MQ=255)
                                                                                                                  aCTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTa                                                                                                                         >  2:722865/1‑62 (MQ=255)
                                                                                                                  aCTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTa                                                                                                                         <  1:722865/62‑1 (MQ=255)
                                                                                                                  aCTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTg                                  >  4:197506/1‑149 (MQ=255)
                                                                                                                   cTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTTGATGGTAATTTATTTCCGGCTTCTTg                                  <  2:140072/148‑1 (MQ=255)
                                                                                                                     g‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTggc                                >  1:660390/1‑148 (MQ=255)
                                                                                                                         aaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCa                            >  3:324546/1‑149 (MQ=255)
                                                                                                                           aaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATg                          <  3:119098/149‑1 (MQ=255)
                                                                                                                            agttGAGCATCATGGGCAAAAAGTCTTTTTTGCCTGTTTTAACCAGATGTTTAACCGCAACGCCGTACGCTACGAAAAAGGTGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTCATTTATGTCCGGATGCTTGGCGGCATgg                         <  4:702170/146‑1 (MQ=255)
                                                                                                                                 aGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttc                    <  4:164652/149‑1 (MQ=255)
                                                                                                                                 aGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttc                    <  4:114452/149‑1 (MQ=255)
                                                                                                                                     tcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctg                 <  1:711534/148‑1 (MQ=255)
                                                                                                                                      catGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCGGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgtt               >  3:417151/1‑149 (MQ=255)
                                                                                                                                      accaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgtt               >  4:792284/5‑149 (MQ=255)
                                                                                                                                         gggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgttcc             <  1:202926/148‑1 (MQ=255)
                                                                                                                                          ggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTgttcctgttcct            >  1:567017/1‑148 (MQ=255)
                                                                                                                                          ggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTgttcctgttcct            >  1:154946/1‑148 (MQ=255)
                                                                                                                                          ggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTgttcctgttcct            >  1:88563/1‑148 (MQ=255)
                                                                                                                                           gCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGttcctgttcctg           <  2:660390/148‑1 (MQ=255)
                                                                                                                                           gCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGc          <  3:578644/149‑1 (MQ=255)
                                                                                                                                              aaaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgtt                     >  1:158477/1‑135 (MQ=255)
                                                                                                                                              aaaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgtt                     <  2:158477/135‑1 (MQ=255)
                                                                                                                                              aaaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGt        >  1:932761/1‑148 (MQ=255)
                                                                                                                                               aaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGtt       <  1:269220/148‑1 (MQ=255)
                                                                                                                                                aaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGttt      >  1:54163/1‑148 (MQ=255)
                                                                                                                                                   gCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGCCg  >  4:363705/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGCCG  >  NC_000913/2588977‑2589271

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑