Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,589,122 G→A 5.5% G510S (GGC→AGC)  acrD → aminoglycoside/multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,589,1220GA 5.5% 412.3 / 8.2 125G510S (GGC→AGC) acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (63/55);  new base A (4/3);  total (67/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 127 total aligned reads displayed.

CGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTT  >  NC_000913/2588974‑2589267
                                                                                                                                                     |                                                                                                                                                 
cggcACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAg                                                                                                                                                   <  4:819780/149‑1 (MQ=255)
 ggcACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAgg                                                                                                                                                  <  3:720815/149‑1 (MQ=255)
   caccacCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACGCTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAGAGGc                                                                                                                                                 >  2:608333/1‑148 (MQ=255)
   caccacCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCt                                                                                                                                                >  3:442638/1‑149 (MQ=255)
      caccGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCtttt                                                                                                                                             >  4:180936/1‑149 (MQ=255)
          ggTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATAATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGc                                                                                                                                         <  3:458849/149‑1 (MQ=255)
           gTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGcc                                                                                                                                        <  3:594718/149‑1 (MQ=255)
           gTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGc                                                                                                                                         <  1:398995/148‑1 (MQ=255)
             gCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCt                                                                                                                                       <  1:186632/148‑1 (MQ=255)
                  cTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGttt                                                                                                                                  >  1:860019/1‑148 (MQ=255)
                  cTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGATTTTTTGCCTGGttt                                                                                                                                  >  2:372880/1‑148 (MQ=255)
                   tATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTaa                                                                                                                                <  4:853018/149‑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                                                                                                                            <  1:424393/148‑1 (MQ=255)
                        ccAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAg                                                                                                                            <  2:941172/148‑1 (MQ=255)
                        ccAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAg                                                                                                                            >  1:877851/1‑148 (MQ=255)
                        ccAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAg                                                                                                                            >  1:75897/1‑148 (MQ=255)
                             tctcTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCTGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGtt                                                                                                                       >  1:132296/1‑148 (MQ=255)
                             tctcTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGtt                                                                                                                       >  1:279132/1‑148 (MQ=255)
                               tctATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTa                                                                                                                     <  2:132296/148‑1 (MQ=255)
                               tctATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTa                                                                                                                     <  1:351376/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‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAAcc                                                                                                                  >  1:96139/1‑148 (MQ=255)
                                  aTTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCg                                                                                                                 >  4:524451/1‑149 (MQ=255)
                                   ttACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCttt                                                                                                                                              >  4:680923/1‑119 (MQ=255)
                                   ttACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCg                                                                                                                 >  2:341829/1‑148 (MQ=255)
                                       cATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGt                                                                                                                                    <  1:37105/125‑1 (MQ=255)
                                       cATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGt                                                                                                                                    >  2:37105/1‑125 (MQ=255)
                                         ttgttgCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGc                                                                                                           >  1:614543/1‑148 (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)
                                                       gTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaaaa                                                                                            >  4:391605/1‑149 (MQ=255)
                                                         gCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGaaaaa                                                                                           <  2:645586/148‑1 (MQ=255)
                                                           tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgg                                                                                         >  1:509011/1‑148 (MQ=255)
                                                           tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgg                                                                                         <  2:413330/148‑1 (MQ=255)
                                                            gTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgggg                                                                                       <  3:180936/149‑1 (MQ=255)
                                                             tCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgggg                                                                                       >  2:14635/1‑148 (MQ=255)
                                                               aGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTg                                                                                     <  2:877851/148‑1 (MQ=255)
                                                                gTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGc                                                                                   >  3:750310/1‑149 (MQ=255)
                                                                 tACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCg                                                                                  >  3:305572/1‑149 (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:187367/1‑149 (MQ=255)
                                                                           gatgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTcc                                                                        >  3:12020/1‑149 (MQ=255)
                                                                            atgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTcc                                                                        >  1:295093/1‑148 (MQ=255)
                                                                            atgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCa                                                                       >  4:125015/1‑149 (MQ=255)
                                                                             tgatCCTCACTCCGGCTCTGTGTGCCACACTACTCAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCAc                                                                      >  4:698036/1‑149 (MQ=255)
                                                                               atCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCg                                                                    >  4:20641/1‑149 (MQ=255)
                                                                                tCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGa                                                                                                        >  2:924522/1‑112 (MQ=255)
                                                                                tCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGa                                                                                                        <  1:924522/112‑1 (MQ=255)
                                                                                 ccTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGt                                                                   <  2:75897/148‑1 (MQ=255)
                                                                                  cTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAg                                                                 >  4:301059/1‑149 (MQ=255)
                                                                                   tCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAg                                                                 <  1:20323/148‑1 (MQ=255)
                                                                                    cACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATAGGCAAAAAAGCCTTTTTGcccatg                                                                                                                                    >  2:668043/1‑76 (MQ=38)
                                                                                    cACTCCGGCTCTGTGTGCCACACTACTTAAGCCAATG‑AAAAAAGGTGAGCATCATAGGCAAAAAAGCCTTTTTGcccatg                                                                                                                                    <  1:668043/80‑5 (MQ=37)
                                                                                     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)
                                                                                                      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‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTa                                            >  2:540610/1‑148 (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                                                                                        >  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‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTg                              >  4:197506/1‑149 (MQ=255)
                                                                                                                      cTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTTGATGGTAATTTATTTCCGGCTTCTTg                              <  2:140072/148‑1 (MQ=255)
                                                                                                                        g‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTggcg                           >  3:114452/1‑149 (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:114452/149‑1 (MQ=255)
                                                                                                                                      catcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcct              <  4:750310/149‑1 (MQ=255)
                                                                                                                                        tcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctg             <  1:711534/148‑1 (MQ=255)
                                                                                                                                         catGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCGGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgtt           >  3:417151/1‑149 (MQ=255)
                                                                                                                                          ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgttc          >  4:615838/4‑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)
                                                                                                                                                  aaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGtt   <  1:269220/148‑1 (MQ=255)
                                                                                                                                                   aaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGttt  >  1:54163/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTG‑AAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTT  >  NC_000913/2588974‑2589267

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