Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,589,122 G→A 5.2% 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.2% 370.7 / 5.6 115G510S (GGC→AGC) acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (69/40);  new base A (4/2);  total (73/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 118 total aligned reads displayed.

CGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGCCG  >  NC_000913/2588974‑2589271
                                                                                                                                                    |                                                                                                                                                      
cggcACCACCGGAGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGGGAGCATCATGGGCAAAAAgg                                                                                                                                                       <  2:1615914/150‑1 (MQ=255)
 ggcACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGc                                                                                                                                                      <  1:977631/150‑1 (MQ=255)
   caccacCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCtt                                                                                                                                                    >  1:706973/1‑150 (MQ=255)
   caccacCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGcct                                                                                                                                                    >  2:1202450/1‑148 (MQ=255)
       accGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCcttttt                                                                                                                                                <  2:684297/150‑1 (MQ=255)
          ggTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGcc                                                                                                                                             <  1:1785636/150‑1 (MQ=255)
           gTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTCGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGccc                                                                                                                                            <  2:191022/150‑2 (MQ=255)
             gCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCGCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTgg                                                                                                                                          <  2:856564/150‑1 (MQ=255)
              ccATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGt                                                                                                                                         <  1:463342/150‑1 (MQ=255)
               cATCTATCGCAAGGTCTATATTTACATTGTGGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGtt                                                                                                                                        <  2:287238/150‑1 (MQ=255)
                aTCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGttt                                                                                                                                       <  1:307992/150‑1 (MQ=255)
                 tCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTa                                                                                                                                      >  1:810719/1‑150 (MQ=255)
                   tATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                                                                                                    >  2:1316271/1‑150 (MQ=255)
                   tATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                                                                                                    >  1:1627288/1‑150 (MQ=255)
                      cGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAGGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAg                                                                                                                                 >  2:1836457/1‑150 (MQ=255)
                             tctcTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTa                                                                                                                          >  1:635531/1‑150 (MQ=255)
                                 tATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCg                                                                                                                      >  2:2046531/1‑150 (MQ=255)
                                     aCCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAAc                                                                                                                  <  2:2021026/150‑1 (MQ=255)
                                     aCCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAAc                                                                                                                  >  2:830101/1‑150 (MQ=255)
                                      ccATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTCCCTGGTTTAACCAGATGTTTAACCGCAACg                                                                                                                 >  1:2211648/1‑150 (MQ=255)
                                        aTTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGcc                                                                                                               >  1:202612/1‑150 (MQ=255)
                                         ttgtttCGGCGATGGTGCTGGCAGTACGGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGAGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCg                                                                                                              <  1:348392/150‑1 (MQ=255)
                                         ttgttgCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGATAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCg                                                                                                              >  1:90814/1‑150 (MQ=255)
                                           gttgCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGaa                                                                                                            >  1:1810541/1‑150 (MQ=255)
                                            ttgCGGCGATGGAGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAAc                                                                                                           <  2:1497801/150‑1 (MQ=255)
                                             tgCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACg                                                                                                          >  2:225734/1‑150 (MQ=255)
                                             tgCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACg                                                                                                          >  1:281397/1‑150 (MQ=255)
                                               cggcgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCt                                                                                                        >  2:1255128/1‑150 (MQ=255)
                                                  cgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACg                                                                                                     >  1:1422879/1‑150 (MQ=255)
                                                      ggTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑aaaa                                                                                                <  1:1235321/150‑1 (MQ=255)
                                                        tGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAg                                                                                              <  2:1024754/150‑1 (MQ=255)
                                                           tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAgggg                                                                                           <  2:1920240/150‑1 (MQ=255)
                                                           tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAgggg                                                                                           >  1:237689/1‑150 (MQ=255)
                                                           tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACAACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAgggg                                                                                           >  1:1187632/1‑150 (MQ=255)
                                                                gTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCg                                                                                      >  2:1806876/1‑150 (MQ=255)
                                                                  aCTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGaa                                                                                    >  1:210314/1‑150 (MQ=255)
                                                                   cTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGaaa                                                                                   >  1:1238594/1‑150 (MQ=255)
                                                                     ggTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAAt                                                                                 >  1:1036990/1‑150 (MQ=255)
                                                                     ggTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAAt                                                                                 >  1:2210153/1‑150 (MQ=255)
                                                                       tAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATtc                                                                               <  1:33917/150‑1 (MQ=255)
                                                                         gCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATtctc                                                                             >  1:527287/1‑150 (MQ=255)
                                                                           gatgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAAGGGGTGGCGAAAATTCTcc                                                                            >  2:1797608/1‑150 (MQ=255)
                                                                           gatgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTccc                                                                           >  1:1121100/1‑149 (MQ=255)
                                                                           gatgatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTccc                                                                           >  1:526014/1‑149 (MQ=255)
                                                                                tCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTa                                                                      <  1:192975/150‑1 (MQ=255)
                                                                                 ccTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAg                                                                     >  1:1832538/1‑150 (MQ=255)
                                                                                 ccTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAAAGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAg                                                                     >  2:1175703/1‑150 (MQ=255)
                                                                                  cTCACTCCGGCTCTGTGTGCCACACTACTTAACCCACTGAAAAAAGGAGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGc                                                                    <  1:568448/150‑1 (MQ=255)
                                                                                   tCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGcc                                                                   <  1:1869413/150‑1 (MQ=255)
                                                                                    cACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCCCCGTAGCCt                                                                  >  2:1628923/1‑150 (MQ=255)
                                                                                     aCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTg                                                                 >  2:1054522/1‑150 (MQ=255)
                                                                                     aCTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATCGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTg                                                                 <  1:1255151/150‑1 (MQ=255)
                                                                                     aCTACGGCTCTGTGTGCCACACTACTTAAGACACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACAAGATGTTTAACCGCAACGCCGAACGCTACT‑AAAAAGGGGTGGCTAAAATTCTCCACCGTAGACTg                                                                 >  2:1893605/1‑150 (MQ=255)
                                                                                         cGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt                                                             >  1:1682549/1‑150 (MQ=255)
                                                                                         cGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt                                                             >  2:2101335/1‑150 (MQ=255)
                                                                                         cGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCt                                                             >  2:232925/1‑150 (MQ=255)
                                                                                          ggCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCAGCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTg                                                            <  2:1151718/150‑1 (MQ=255)
                                                                                          ggCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCAGCATGGGCAAAAAGGCTTTTTTGCCGGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTg                                                            <  1:406833/150‑1 (MQ=255)
                                                                                               tgtgtgCCACACTACTTAAGCCACTTAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATtg                                                       <  1:1156880/150‑1 (MQ=255)
                                                                                                   tgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAt                                                   >  1:1208296/1‑150 (MQ=255)
                                                                                                    gCCACACTACTTAAGCCTCTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAtt                                                  <  2:1015028/150‑1 (MQ=255)
                                                                                                      cacaCTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAtttt                                                >  2:123796/1‑149 (MQ=255)
                                                                                                      cacaCTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTa                                                <  1:1583067/150‑1 (MQ=255)
                                                                                                       acacTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTAt                                               >  2:1771059/1‑150 (MQ=255)
                                                                                                           tactTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTcc                                           <  1:1660880/150‑1 (MQ=255)
                                                                                                           tactTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTcc                                           <  1:1234907/150‑1 (MQ=255)
                                                                                                             ctTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCctg                                         >  2:1035677/1‑150 (MQ=255)
                                                                                                              ttAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCctgc                                        <  1:1391812/150‑1 (MQ=255)
                                                                                                               tAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCctgct                                       >  2:726219/1‑150 (MQ=255)
                                                                                                                  gCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCGGGATTGTGATTTATGTCCTgctgct                                    <  1:2040675/150‑1 (MQ=255)
                                                                                                                   ccACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCtt                                   >  1:1810249/1‑150 (MQ=255)
                                                                                                                     aCTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCACCGCCGAACGCTACG‑AAAACGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTgg                                 >  1:81901/1‑150 (MQ=255)
                                                                                                                     aCTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTgg                                 >  1:1305871/1‑150 (MQ=255)
                                                                                                                     aCTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTgg                                 <  2:794532/150‑1 (MQ=255)
                                                                                                                      cTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTggc                                <  1:2041375/150‑1 (MQ=255)
                                                                                                                      cTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTggc                                <  1:1054522/150‑1 (MQ=255)
                                                                                                                       tGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTggcg                               >  1:1996946/1‑150 (MQ=255)
                                                                                                                        gAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTggcgg                              >  1:1014332/1‑150 (MQ=255)
                                                                                                                        gAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTggcgg                              >  2:542149/1‑150 (MQ=255)
                                                                                                                          aaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCa                            >  1:1141764/1‑150 (MQ=255)
                                                                                                                            aaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATg                          >  1:954762/1‑150 (MQ=255)
                                                                                                                               ggTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGgtg                       >  2:1470088/1‑150 (MQ=255)
                                                                                                                                gTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGgtgt                      <  1:1316271/150‑1 (MQ=255)
                                                                                                                                 tGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgtt                     >  2:537529/1‑150 (MQ=255)
                                                                                                                                  gAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttc                    >  1:724331/1‑150 (MQ=255)
                                                                                                                                    gCATCATGGGCAAAACGGCGTTTTTGCCTGGGTTAACCAGATGTTTAACCGCACCGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcct                  <  1:1375703/150‑1 (MQ=255)
                                                                                                                                      atcatGGGCAAAAAGGCTTTTTTGCCTGGGTGAACCAGATGTTTAACCGCAACGCCGCACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgt                <  1:349028/150‑1 (MQ=255)
                                                                                                                                         ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgttcc             >  1:2170954/4‑150 (MQ=255)
                                                                                                                                          tGGGCAATAAGGCTTTTTTGCCTGGTTTAACCAGCTGTTTAACCGCAACGCCGAACGCTACG‑ATAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTgttcctgttcct            >  1:1431722/1‑150 (MQ=255)
                                                                                                                                          tGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTgttcctgttcct            >  2:1207137/1‑150 (MQ=255)
                                                                                                                                          caggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCATGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTgatcctgttcct            >  1:441659/3‑150 (MQ=255)
                                                                                                                                            ggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGc          >  2:1232325/1‑150 (MQ=255)
                                                                                                                                            ggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGc          >  1:693203/1‑150 (MQ=255)
                                                                                                                                              cAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGt        >  1:942311/1‑150 (MQ=255)
                                                                                                                                                aaaaGGC‑TTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTg     >  1:680527/1‑150 (MQ=255)
                                                                                                                                                 aaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGATACG‑ATAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTCCTTGGCGGCATGGTGTTCATGTTCATGCTTTTg     >  2:335885/1‑150 (MQ=255)
                                                                                                                                                  aaGGCTTTTTTGGCTGGTTTAACCAGATGTTTATCCGCTCCGCCGAACGCTACG‑AAAATGGGGGGGCGAAAATTCGCGACCGTAGCCTGCGCTGGATGGTGAGTTATGGCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGc    <  1:222560/150‑1 (MQ=255)
                                                                                                                                                   aGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGcc   <  1:537729/150‑1 (MQ=255)
                                                                                                                                                    ggCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGCCg  >  1:1335220/1‑150 (MQ=255)
                                                                                                                                                    ggCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGCCg  >  2:1725962/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CGGCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACG‑AAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGCCG  >  NC_000913/2588974‑2589271

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑