Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,356,2010TG41.1% 184.4 / 33.3 108T304P (ACC→CCC) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (9/54);  new base G (44/0);  total (54/54)
Fisher's exact test for biased strand distribution p-value = 1.90e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

GTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCC  >  NC_000913/1356055‑1356351
                                                                                                                                                    |                                                                                                                                                        
gTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTa                                                                                                                                                      <  2:1538611/151‑1 (MQ=255)
   gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTaaa                                                                                                                                                    <  2:1930682/150‑1 (MQ=255)
   gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAAt                                                                                                                                                   <  1:539680/151‑1 (MQ=255)
     cATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGc                                                                                                                                                 <  2:843824/151‑1 (MQ=255)
        aGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCa                                                                                                                                               <  1:1366795/150‑1 (MQ=255)
           cccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATa                                                                                                                                           <  1:172278/151‑1 (MQ=255)
            ccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTAAATCCAGATa                                                                                                                                           >  1:1137629/1‑150 (MQ=255)
            ccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAATGCCAGa                                                                                                                                             >  2:846473/1‑148 (MQ=255)
              aCGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGc                                                                                                                                        <  1:1475355/151‑1 (MQ=255)
                gAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGa                                                                                                                                      <  1:2104893/151‑1 (MQ=255)
                         ggTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATGCCAGAAAAGGGACGTTCATcc                                                                                                                            >  1:846949/1‑151 (MQ=255)
                            aaaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCa                                                                                                                          <  1:2091747/150‑1 (MQ=255)
                                  gagGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCa                                                                                                                   <  2:1309395/151‑1 (MQ=255)
                                   agGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGTGTTTAATGACAGAGAgg                                                                                                                                         >  2:1911903/1‑129 (MQ=255)
                                      ccgccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAAAGCCCAATAGGCGACGGTCATTCCCGGACGCAGCg                                                                                                                >  2:1550706/1‑150 (MQ=255)
                                       cgccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAATGGCAGAAAAGCGGAGTTTCTCCCAgggagccgcggt                                                                                                              >  2:2064886/1‑141 (MQ=255)
                                       cgccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTAAATGCCAGATAAGGGACGTTCATCCCCGGGCGCAGc                                                                                                                 >  2:1221036/1‑148 (MQ=255)
                                        gccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGCGGGGGTAAAAGACAGGAAGGCGAACCTCATTCCAAGGCGCAgg                                                                                                                 >  2:189425/1‑146 (MQ=255)
                                              tCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGc                                                                                                        <  1:372389/150‑1 (MQ=255)
                                               cAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTTAAAGCCCGATAAGGGGCGTTCAACCCCGGACGCAGCGTTAAAAGc                                                                                                        >  2:1822089/1‑149 (MQ=255)
                                               cAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAGCGCaa                                                                                                      <  1:1169179/151‑1 (MQ=255)
                                               cAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCa                                                                                                       <  1:284135/150‑1 (MQ=255)
                                                 gCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTAAAAGCCAGATAGGCGCCGTTTATCCCAGGACGCCGCGTTaaa                                                                                                           >  2:1582976/1‑144 (MQ=255)
                                                 gCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGGGGGGGTAAATGCACGAAAGGGGAAGTTTCTCCCAGGGCGCAGCGGTaaa                                                                                                           >  2:2095658/1‑144 (MQ=255)
                                                   cGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAAATCCAGAAAGGCGACGTTTATCCCCGGACGCCGCGGTAAAAGc                                                                                                        >  2:1423347/1‑145 (MQ=255)
                                                    gTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTAAATACCAGATAGGGGGCGGTCATACCCGGAGGCAGCGGTTAAAGCAAgcg                                                                                                  >  1:1916710/1‑149 (MQ=255)
                                                      aCCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCggg                                                                                              <  1:1887407/151‑1 (MQ=255)
                                                      aCCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCgg                                                                                               <  1:739134/150‑1 (MQ=255)
                                                       ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgt                                                                                             <  2:1464679/151‑1 (MQ=255)
                                                        cATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAATGCCAGATAGGGGACGTTCATCCCAGGACGCAGCGTTAAAAGCAAGCGCggg                                                                                              >  1:2042391/1‑149 (MQ=255)
                                                         ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAATGCCAGAAAAGCGACGTTCATCCCCAGACGCAGCGTTAAAAGCAAGCGCGGggcg                                                                                           >  1:1584832/1‑151 (MQ=255)
                                                         ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAAAGCCAGATAGGCGACGGTCCTCCCCAGACGCAGCGTTAAACGCCAACGCgggg                                                                                             >  2:1062608/1‑148 (MQ=255)
                                                         ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTTAAAGCCAGAAAGGCGGCGGTTATCCCCAGACGCCGGCGTAAAAGCCAACGCggg                                                                                              >  2:1007859/1‑148 (MQ=255)
                                                         ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGGGGGGTTAAAATCCCGATAAGGGACCTTCCTCCCCCGACACCGCGGTAAAAACAAGGGCggg                                                                                              >  2:1339220/1‑148 (MQ=255)
                                                           aataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgtcgtc                                                                                         <  2:743661/151‑1 (MQ=255)
                                                            ataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCa                                                                                        <  2:2023107/151‑1 (MQ=255)
                                                            ataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAATGCCAGAAAGGGGGCGTTCATCCCCGGGCGCCGCGGTAAAACCAACCGCggg                                                                                              >  2:812944/1‑145 (MQ=255)
                                                             taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTAAATGCCAGATAGGCGGCGGTCATCCCCGGACGCAGCGGTAAAAGCAAACGCggg                                                                                              >  1:1794531/1‑144 (MQ=255)
                                                             taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAAAGCCCGATAAGGGGCGGTCATCCCCGGACGCCGCGGTAAAAGCAAAGGCGGGGCGGCAc                                                                                       >  1:1748276/1‑151 (MQ=255)
                                                              aGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACg                                                                                      <  2:1218250/151‑1 (MQ=255)
                                                                aTGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCa                                                                                    <  2:1545458/151‑1 (MQ=255)
                                                                  ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCaaa                                                                                  <  1:1982025/151‑1 (MQ=255)
                                                                  ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGGTAAATACCAGATAGGCGACGGTCATCCCCGGGCGCAGCGTTAAAAGCAAACGCGGGGCGGCACGCaa                                                                                   >  2:435316/1‑150 (MQ=255)
                                                                   gATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCaaaa                                                                                 <  1:610089/151‑1 (MQ=255)
                                                                   gATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCaaaa                                                                                 <  2:732078/151‑1 (MQ=255)
                                                                     ttGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGGGGGGTTAAAAGGCAGATAGGGGACGTTTATTCCCAGACGCAGCGTTAAAAGCAAGCGCGGGGGGTCACGCAAAATg                                                                               >  2:1049864/1‑151 (MQ=255)
                                                                      tGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAATGCCAGATAGGGGACGGTCATCCCAGGAAGCAGCGGTGAAAGCAAACGCGGGGCGGCACGCaaaaa                                                                                >  2:111537/1‑148 (MQ=255)
                                                                       gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGa                                                                             <  2:1513474/151‑1 (MQ=255)
                                                                         aTCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATa                                                                           <  2:889186/151‑1 (MQ=255)
                                                                            aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATGCCAGATAGGCGACGGTCCTCCCAGGACGCCGCGGTAAAAGCAAACGCGGGGCGTTACGCAAAAAgg                                                                              >  2:1246426/1‑145 (MQ=255)
                                                                             ggCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctg                                                                       <  2:458210/151‑1 (MQ=255)
                                                                              gcgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctgg                                                                      <  2:1725057/151‑1 (MQ=255)
                                                                              gcgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATGCCAGACAGGCGACCTTCTTCCCAGGACGCAGCGGTTAACGCAAACGCGGGGCGGCACGCAAAAAgggaag                                                                          >  2:643825/1‑143 (MQ=255)
                                                                               cgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctggc                                                                     <  1:1052793/151‑1 (MQ=255)
                                                                                gcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAATGCCAGATAGGGGACGGTCATCCCAGGACGCAGCGGTAAACGCAAACGCGGGGCGTCACGCAAAATGGATAactgg                                                                      >  1:1784960/1‑149 (MQ=255)
                                                                                  tGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAAAGCCAGAAAGGCGGCGTTCATCCCAGGACGCAGCGGTAAACGCAAGCGCGGGGCGGCACGCCAAATGGATAGctgggtgg                                                                  >  1:375798/1‑151 (MQ=255)
                                                                                  tGGCTATTAACCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTGGCCGTGTTTAATGCCCGATAGGCGACGGTCATCCCAAGACACAGCGTTAAAACCAAACGGGGGGCGGCACGCaaaaaaggaaagtggg                                                                     >  2:2142503/1‑136 (MQ=255)
                                                                                     ttattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAATGCCAGATAGGCGACGTTCATCCCCGAGCGCCGCGGTAAACGCCACCGCGGGGCGGCACGCCAAAAGGGAAActggctgg                                                                  >  2:485974/1‑148 (MQ=255)
                                                                                       attaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcagca                                                             <  2:2092375/151‑1 (MQ=255)
                                                                                       attaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcagca                                                             <  1:1620598/151‑1 (MQ=255)
                                                                                         taATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGcagcag                                                            <  2:2079185/150‑1 (MQ=255)
                                                                                          aaTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTTAAAGCCAGATAAGCGGCGTTCATCCCCGGACGCAGCGTTAAACGCAAGCGCGGGGCGTCCCGCAAAAAGGGTAGctggctgg                                                                  >  1:25041/1‑143 (MQ=255)
                                                                                           aTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTGGGGGGGGTTAAAGGCCGATAGGCGGGGGTCATCCCCGGACGCAGCGTTAAAAGCAAGCGCGGGGCGGCACGCCAAATGGGTAGCGGGCTGGCAGcaggcc                                                         >  1:544898/1‑151 (MQ=255)
                                                                                            tCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGCTTGGCGGCGGTAAATGCCCCAAATGCGACGCTTATTCCCCGGCCCCGCGTTAAACACACGAGCGGGGCGGCACGCAAAAAAGATAactggggtcgggcccgcc                                                         >  2:1754751/1‑137 (MQ=255)
                                                                                             cGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggcca                                                        <  1:1972011/150‑1 (MQ=255)
                                                                                                cAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAATGCCAGATAGGGGACGATTATCCCCGGACGCAGCGCTTAAAGCCAGCGCGGGGCGGCACGCAAAAAAGAGAGCTGGCTGGCAGccgggcagg                                                      >  2:2112288/1‑149 (MQ=255)
                                                                                                cAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                    <  1:1963353/151‑1 (MQ=255)
                                                                                                 aGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAATGCCCGATAGGCGGCGGTCATCCCCGGACGCAGCGTTAAAAGCAAGCGCGGGTCGGCACGCAAAAAGGATAGCTGGCTGGCAGCagggcaggcc                                                    >  1:2144213/1‑150 (MQ=255)
                                                                                                 aGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAAAGCCAGATAAGCGGCGGTCATCCCAGGACGCCGCGTTAAAAGCAAACGCGGGGCGGCACGCAAAAAGGATAGCTGGCTGGCAGCaggccaggcc                                                    >  1:579819/1‑150 (MQ=255)
                                                                                                    gCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGa                                                 <  2:457600/150‑1 (MQ=255)
                                                                                                     ccAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAATGCCAGATAGGCGACGTTCATCCCAGGAAGCAGCGGTAAAAGCAAGCGCGGGTCGTCACGCAAAATAGATAGCTGGCTGGCGACCGGGCAGGCCCGaa                                                >  1:386099/1‑150 (MQ=255)
                                                                                                     ccAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAAAGCCAGATAGGGGACGGTCATCCCCCGACGCAACGGTTAACGCAAGCGCGGGGTGTCACGCAAAAAAGATAGCTGGCTGGGGGGAAGGCAGGCCAGaa                                                >  1:699726/1‑150 (MQ=255)
                                                                                                      cAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAATGCCAGATAGGCGACGGTCATCCCCGGACGCAGCGTTAAAAGCAAGCGCGGGGCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCCGGCCAGaa                                                >  1:1701127/1‑149 (MQ=255)
                                                                                                       aGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGTGGTAAAAGCCAGATAGGGGACGTTCCTCCCCGGACGCAGCGTTAAAAGCAAAGGCGGGTCGCCACGCAAAAAGGAAAGctggctgg                                                                  >  2:931815/1‑130 (MQ=255)
                                                                                                        gcgcATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAATGGCAGGAAGGGGACGTTTCTCCCCAGGCCGAGCGGTAAAAGCAAACGGGGGGCGGCCCGCAAAAAGGGAAGGTGGCTGGGAGCcgggccggcccgga                                                >  2:1415531/1‑143 (MQ=255)
                                                                                                          gcATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGc                                          <  1:1421448/151‑1 (MQ=255)
                                                                                                           cATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCa                                         <  1:1454148/151‑1 (MQ=255)
                                                                                                             tGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGTGGGGTTTAAAGCCCGGGAGGGGACGTTGATCCCCGGGCCCAGCGTTCAACGCAAACGCGGGGCGGTACGCAAAAAGGATAACTCGCCGGGCGCAGGGCCGGCCaaaa                                                >  2:185218/1‑142 (MQ=255)
                                                                                                              ggCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTTAATGCCAGAGAGGGGGCGGTCATCCCAGGACGCCGGGTTAAAAGCAAGAGCGGGGCGGCCCGCAAAAAGGGTAGACGGGTGGGAGGAGGCCCGGCCCCAACGgcgc                                          >  2:1743785/1‑147 (MQ=255)
                                                                                                              ggCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAAAGCCAGATAGGCGGCGGTCATCCCAGGACGCAGCGGTAAAAGCAAGCGCGGGGCGTCACGCAAAAAGGATAGCTGGCTGGGAGCAGGCCAGGCCCGAACGTCGCAtt                                       >  1:1695945/1‑150 (MQ=255)
                                                                                                               gCGGACAGCGGGATTATTTAGCGGCGGTTTTGAGGTGGTAAAAATCAGATAGGAGACGCTCAACCCAGGGGGCCGCGTTAAAAGCAAGACCGGGGCCGCCACCAAAAAGGGAAAGGCGCTGGCAGCCGGCCCGGCCACAACgaccccttc                                      >  1:763681/1‑142 (MQ=255)
                                                                                                               gCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTcc                                     <  2:749994/151‑1 (MQ=255)
                                                                                                                     aGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTc                               <  2:201096/151‑1 (MQ=255)
                                                                                                                     aGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAAAACCAGATAGGCGACGGTCATCCCCGGACGCCGCGTTAAAAGCAAGGGCGGGGGGGCACGCAAAAATGATAGCTGGGTGGCAGCAGGCCCGGCCCGAAAGTCGCATTcccccg                                 >  1:1834711/1‑147 (MQ=255)
                                                                                                                      gCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCa                              <  2:1474781/151‑1 (MQ=255)
                                                                                                                       cGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATGCCCGATAAGCGACGTTCATCCCCAGACGCAGCGTTAAAAGCAAGCGCGGGGCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGGCGCATTcccccg                                 >  1:99831/1‑145 (MQ=255)
                                                                                                                                  aGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACa                  <  1:890831/151‑1 (MQ=255)
                                                                                                                                    cggcggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGa                <  2:959547/151‑1 (MQ=255)
                                                                                                                                    cggcggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGa                <  2:1137165/151‑1 (MQ=255)
                                                                                                                                     ggcggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGAc               <  2:347917/151‑1 (MQ=255)
                                                                                                                                      gcggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACg              <  2:1345615/151‑1 (MQ=255)
                                                                                                                                      gcggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACg              <  1:787068/151‑1 (MQ=255)
                                                                                                                                         ggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTc            >  2:678395/1‑151 (MQ=255)
                                                                                                                                              tGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCgg         <  1:1603765/150‑1 (MQ=255)
                                                                                                                                              tGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGt        >  1:103922/1‑151 (MQ=255)
                                                                                                                                               ggCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTg       <  1:1051756/151‑1 (MQ=255)
                                                                                                                                                 cGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgcc     <  1:157299/151‑1 (MQ=255)
                                                                                                                                                  ggTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccg    >  1:883781/1‑151 (MQ=255)
                                                                                                                                                  ggTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgcc     <  2:2198065/150‑1 (MQ=255)
                                                                                                                                                    tgtTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgcc  >  2:1001112/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                        
GTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCC  >  NC_000913/1356055‑1356351

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: