Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,356,2010TG37.7% 149.1 / 10.9 78T304P (ACC→CCC) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (5/43);  new base G (29/0);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 2.10e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCC  >  NC_000913/1356054‑1356351
                                                                                                                                                   |                                                                                                                                                         
cGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGgtg                                                                                                                                                          <  1:543419/149‑1 (MQ=255)
cGTTGTCATAGGCCCACGAGGCGAGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGgtg                                                                                                                                                          >  2:436195/1‑149 (MQ=255)
  ttGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGtt                                                                                                                                                        >  1:821799/1‑149 (MQ=255)
   tGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTaa                                                                                                                                                      <  1:518674/150‑1 (MQ=255)
   tGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTaa                                                                                                                                                      <  1:382549/150‑1 (MQ=255)
    gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTaaa                                                                                                                                                     <  2:960334/150‑1 (MQ=255)
      cATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGc                                                                                                                                                  <  1:225146/151‑1 (MQ=255)
        tAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTTAATggc                                                                                                                                                 >  1:1228116/1‑148 (MQ=255)
          ggCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGa                                                                                                                                              <  2:908417/151‑1 (MQ=255)
          ggCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAAGCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGTCGGTGTTAAATGCCAGa                                                                                                                                              <  2:1633319/151‑1 (MQ=255)
           gCCCACGAGGCGCGCTGTAAAATAGATGCCTCCGTTTCAGCCGTACCAGAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGTACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGAt                                                                                                                                             <  2:1277636/151‑1 (MQ=255)
           gCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTAAAAGcc                                                                                                                                                 >  1:561163/1‑147 (MQ=255)
              cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAAAGCCAGATa                                                                                                                                            >  1:407746/1‑149 (MQ=255)
              cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAATGCCAGATa                                                                                                                                            >  1:580307/1‑149 (MQ=255)
                  aGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATGCCAGa                                                                                                                                              >  1:1405688/1‑143 (MQ=255)
                       cgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGCGGGGGTTAATGCCAGATAAGCGACggtt                                                                                                                                  >  2:299409/1‑147 (MQ=255)
                            tAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTTAATAACAGAAAgggggaggttaacccc                                                                                                                            >  2:1160887/1‑137 (MQ=255)
                             aaaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTTAAAGCCAGATAGGCGGCGGTCATccc                                                                                                                           >  1:111833/1‑149 (MQ=255)
                             aaaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCa                                                                                                                          <  1:604329/150‑1 (MQ=255)
                                  agagGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAAAGCCCGATAGGGGACGGTCCTccc                                                                                                                           >  1:1443786/1‑144 (MQ=255)
                                           cgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTaaa                                                                                                           <  2:1495949/151‑1 (MQ=255)
                                           cgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGGTAAATGCCAGAGAAGCGACGGTCACCCCAGGACGCCGCGTTa                                                                                                             >  2:535609/1‑149 (MQ=255)
                                            gTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTAATACCAAATAGATGACGTTCATTCCCAgccggccggttaaa                                                                                                            >  2:1239422/1‑137 (MQ=255)
                                             tttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTTAAAGCCAGATAGGGGACGTTCATCCCAGGACGCAGCGTTAAAc                                                                                                          >  2:1036618/1‑150 (MQ=255)
                                              ttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTTAATGCCAGATAGGCGACGTTCATCCCCGGACGCAGCGGTaaa                                                                                                           >  1:686116/1‑148 (MQ=255)
                                              ttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACg                                                                                                         <  2:193389/150‑1 (MQ=255)
                                                    cGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAATGCCAAATAGGCGACGTTTATCCCCAGACGCAGCGGTTAAAGCAAacg                                                                                                   >  2:967110/1‑150 (MQ=255)
                                                    cGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAgcgc                                                                                                  <  2:835437/151‑1 (MQ=255)
                                                     ctACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTAAAAGCCAGATAGGCGGCGTTCCTCCCAGGGCGCAGCGTTTAAAGCaa                                                                                                      >  1:1161803/2‑146 (MQ=255)
                                                      tACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTAAAAGCCCGATTCGCGGAGGTCATTCCCGGGCGCCAGGGTCAACGCAAGAGCgg                                                                                               >  1:1206653/1‑151 (MQ=255)
                                                       aCCATAATAGATGGATTGCATCAGGCGCTGGGTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGAGTATTGAGCGGCGGTTTGGCGGTGTTAAATGCCCGATAGGCGCCGTGCATCCCAGGACGCAGCGTTAAACGCAATcgcg                                                                                                 <  2:514636/149‑1 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgt                                                                                             <  1:638044/151‑1 (MQ=255)
                                                             ataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGATTGGCGGGGTTTAATGCCCGATATGCGACGTTTCGCCCCAGGCGCCGCGTTAAAAGCAAGCGCgg                                                                                               >  1:37962/1‑144 (MQ=255)
                                                              taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCa                                                                                        <  2:269672/150‑1 (MQ=255)
                                                                gATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAATGCCAGATAAGCGACGTTTATCCCCGGACGCAGCGGTTAACGCAAGCGCGGGGCGTCAc                                                                                       >  1:1538477/1‑149 (MQ=255)
                                                                  tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTTAATGCCAGATAGGCGGCGGTCATTCCAGGGCGCAGCGGTTAAAGCAAGCGCGGggcgtccc                                                                                       >  1:602799/1‑145 (MQ=255)
                                                                     aTTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAAAGCCAGATAGGCGGCGTTCCTCCCCGGACGCCGGGGTTAAAGCAAACGCGGGGCGGCACGCaaa                                                                                  >  1:864132/1‑149 (MQ=255)
                                                                        gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGa                                                                             <  2:140466/151‑1 (MQ=255)
                                                                         cATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGa                                                                             <  1:1039602/150‑1 (MQ=255)
                                                                           tCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTTAATGCCAGATAAGCGGCGTTTATCCCCGGACGCAGCGTTAAACGCAAGCGCGGGGGGGCACGCaaaa                                                                                 >  1:784327/1‑144 (MQ=255)
                                                                             aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgct                                                                        <  1:1521716/151‑1 (MQ=255)
                                                                             aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAAAGCCAGATAGGCGACGTTCATCCCCGGACGCCGCGTTAAAAGCAAACGCGGGGCGGCACGCAAAATGGATa                                                                           >  1:829729/1‑148 (MQ=255)
                                                                             aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGCGGGGGTTAATGGCCAATAAGCGACGGTCATCCCACGACGCCGCGGTTAACGCAAGCGCGGGTGGTTCCGCGAAATGGAAAgc                                                                         >  1:1120551/1‑150 (MQ=255)
                                                                                    ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAATGCCAGATAAGCGACGGTCATCCCCGGACGCAGCGTTAAAAGCAAGCGCGGGGCGTCACGCCAAAAGGATAGcttgctgg                                                                  >  1:1100066/1‑150 (MQ=255)
                                                                                     gTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTAAAAGCCAGAAAGGGGAAGTTCATCCCCGGACGCACAGGTAAAAGGAAGCGCGGGGCGGTACGCAAAATGGGTa                                                                           >  2:904181/1‑140 (MQ=255)
                                                                                     gTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgca                                                                <  1:1569251/151‑1 (MQ=255)
                                                                                       tattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcag                                                               <  2:109569/150‑1 (MQ=255)
                                                                                        attaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcagc                                                              <  1:1059635/150‑1 (MQ=255)
                                                                                           aaTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTTAAAGGCAGAAAAGCGGCGTTTATCCCCAGGCGCAGCGGTAAACGCCAAGGCGGGGCGTCACGCAAAAAGGATAggtgg                                                                      >  2:699010/1‑139 (MQ=255)
                                                                                             tCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTAAAAGACAGAAAAGGGAAGGTTAACCCAGGGCGCCGCGCTAAACGCAAACGCGGGGCGTCCCGCAAAAAgggaaagtgggcgg                                                                  >  2:1460238/1‑129 (MQ=255)
                                                                                             tCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGCGGGGGTAAAACACCGAGAGGCGCAGTTCCTCCCAGGACCGAGCGTTAAACACAAACGCCGGTAGCCCCGCAAAAAAGTTAGGCTGCTGGcagcagc                                                           >  2:914859/1‑147 (MQ=255)
                                                                                              cGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggccag                                                       <  1:1150396/151‑1 (MQ=255)
                                                                                               gCCACTGCCAGCGCAGGGCGGACATCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAACAGCCAGATAGGCGACGTTCATACCAGGACGCAGCGTTGAACGCATGCGCAGTGCGTCGCGTACAATGGATAGCTGTCTCGCAGcaggcca                                                        <  2:205914/149‑1 (MQ=255)
                                                                                                   gTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGTTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGgcaaggacag                                                  <  2:569722/151‑1 (MQ=255)
                                                                                                        aGCGCATGGCGGACATCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGGTAGCTGGCTGGGAGCAGGCCAGGCCAGAACGt                                             <  2:81062/151‑1 (MQ=255)
                                                                                                         gcgcATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTc                                            <  1:1037457/151‑1 (MQ=255)
                                                                                                         gcgcATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTc                                            <  2:1101555/151‑1 (MQ=255)
                                                                                                            cATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCa                                         <  1:290067/151‑1 (MQ=255)
                                                                                                             aTGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGTGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCAt                                        <  2:367438/151‑1 (MQ=255)
                                                                                                               ggCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTc                                      <  1:383858/151‑1 (MQ=255)
                                                                                                                     cAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGt                                <  1:205433/151‑1 (MQ=255)
                                                                                                                       gCGGGATTATTTAGCGGCGGTTTGGTGGCGTTATATACCCAATAGGGGTCCTTCATTCCCGGACGGAGCGGTAAACGCAAGAGCGCGGTGGGCCCCAAAATGGGGAGCTGGCTGGCAGCAGGCCAGGGCCGAACgcccgccccccccggc                               >  2:103687/1‑135 (MQ=255)
                                                                                                                        cGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAg                             <  1:1377602/151‑1 (MQ=255)
                                                                                                                            attattTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCGCGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGt                         <  1:704612/151‑1 (MQ=255)
                                                                                                                            attattTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGt                         <  1:440225/151‑1 (MQ=255)
                                                                                                                            attattTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGt                         <  2:908766/151‑1 (MQ=255)
                                                                                                                            attattTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGt                         <  1:1327699/151‑1 (MQ=255)
                                                                                                                               attTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGttt                       <  2:919697/150‑1 (MQ=255)
                                                                                                                                 ttAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCg                     <  1:210264/150‑1 (MQ=255)
                                                                                                                                  tAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGa                    <  1:292721/150‑1 (MQ=255)
                                                                                                                                   aGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACa                  <  1:20260/151‑1 (MQ=255)
                                                                                                                                    gcggcgGTTTGGCGGGGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGGCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGAc                   >  1:182330/1‑149 (MQ=255)
                                                                                                                                     tggcggTCTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGAc                   >  2:566588/2‑148 (MQ=255)
                                                                                                                                     cggcggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGa                    >  2:1519765/1‑147 (MQ=255)
                                                                                                                                         ggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTc            >  2:1327970/1‑151 (MQ=255)
                                                                                                                                          gTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTcc           <  1:1506727/151‑1 (MQ=255)
                                                                                                                                             tGGCGGTGTTAAATGCCAGATAGGCGACGTTCAGCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCgg         <  2:873292/150‑1 (MQ=255)
                                                                                                                                              ggCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGt        <  1:1633604/150‑1 (MQ=255)
                                                                                                                                              ggCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTg       <  1:908705/151‑1 (MQ=255)
                                                                                                                                                 ggTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccg    >  1:1228127/1‑151 (MQ=255)
                                                                                                                                                  gtgtTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgc   >  2:257179/1‑151 (MQ=255)
                                                                                                                                                   tgtTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgcc  >  1:362795/1‑151 (MQ=255)
                                                                                                                                                   tgtTAAATGCCAGATAGGCGACGGTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccg    <  2:1248764/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                         
CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCC  >  NC_000913/1356054‑1356351

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