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

CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCC  >  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)
                             aaaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTTAAAGCCAGATAGGCGGCG‑GTCA‑Tccc                                                                                                                           >  1:111833/1‑149 (MQ=255)
                             aaaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCa                                                                                                                          <  1:604329/150‑1 (MQ=255)
                                  agagGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAAAGCCCGATAGGGGACG‑GTCC‑Tccc                                                                                                                           >  1:1443786/1‑144 (MQ=255)
                                           cgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTaaa                                                                                                           <  2:1495949/151‑1 (MQ=255)
                                           cgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGGTAAATGCCAGAGAAGCGACG‑GTCA‑CCCCAGGACGCCGCGTTa                                                                                                             >  2:535609/1‑149 (MQ=255)
                                            gTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTAATACCAAATAGATGACG‑TTCATTCCCAgccggccggttaaa                                                                                                            >  2:1239422/1‑137 (MQ=255)
                                             tttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTTAAAGCCAGATAGGGGACG‑TTCA‑TCCCAGGACGCAGCGTTAAAc                                                                                                          >  2:1036618/1‑150 (MQ=255)
                                              ttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTTAATGCCAGATAGGCGACG‑TTCA‑TCCCCGGACGCAGCGGTaaa                                                                                                           >  1:686116/1‑148 (MQ=255)
                                              ttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACg                                                                                                         <  2:193389/150‑1 (MQ=255)
                                                    cGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAATGCCAAATAGGCGACG‑TTTA‑TCCCCAGACGCAGCGGTTAAAGCAAacg                                                                                                   >  2:967110/1‑150 (MQ=255)
                                                    cGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAgcgc                                                                                                  <  2:835437/151‑1 (MQ=255)
                                                     ctACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTAAAAGCCAGATAGGCGGCG‑TTCC‑TCCCAGGGCGCAGCGTTTAAAGCaa                                                                                                      >  1:1161803/2‑146 (MQ=255)
                                                      tACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTAAAAGCCCGATTCGCGGAG‑GTCA‑TTCCCGGGCGCCAGGGTCAACGCAAGAGCg‑g                                                                                               >  1:1206653/1‑151 (MQ=255)
                                                       aCCATAATAGATGGATTGCATCAGGCGCTGGGTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGAGTATTGAGCGGCGGTTTGGCGGTGTTAAATGCCCGATAGGCGCCG‑TGCA‑TCCCAGGACGCAGCGTTAAACGCAATcgcg                                                                                                 <  2:514636/149‑1 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑Ggt                                                                                             <  1:638044/151‑1 (MQ=255)
                                                             ataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGATTGGCGGGGTTTAATGCCCGATATGCGACGTTTCG‑CCCCAGG‑CGCCGCGTTAAAAGCAAGCGCg‑g                                                                                               >  1:37962/1‑144 (MQ=255)
                                                              taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCa                                                                                        <  2:269672/150‑1 (MQ=255)
                                                                gATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAATGCCAGATAAGCGACG‑TTTA‑TCCCCGGACGCAGCGGTTAACGCAAGCGCG‑GGGCGTCAc                                                                                       >  1:1538477/1‑149 (MQ=255)
                                                                  tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTTAATGCCAGATAGGCGGCG‑GTCA‑TTCCAGGGCGCAGCGGTTAAAGCAAGCGCG‑Gggcgtccc                                                                                       >  1:602799/1‑145 (MQ=255)
                                                                     aTTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAAAGCCAGATAGGCGGCG‑TTCC‑TCCCCGGACGCCGGGGTTAAAGCAAACGCG‑GGGCGGCACGCaaa                                                                                  >  1:864132/1‑149 (MQ=255)
                                                                        gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGa                                                                             <  2:140466/151‑1 (MQ=255)
                                                                         cATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGa                                                                             <  1:1039602/150‑1 (MQ=255)
                                                                           tCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTTAATGCCAGATAAGCGGCG‑TTTA‑TCCCCGGACGCAGCGTTAAACGCAAGCGCG‑GGGGGGCACGCaaaa                                                                                 >  1:784327/1‑144 (MQ=255)
                                                                             aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAgct                                                                        <  1:1521716/151‑1 (MQ=255)
                                                                             aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAAAGCCAGATAGGCGACG‑TTCA‑TCCCCGGACGCCGCGTTAAAAGCAAACGCG‑GGGCGGCACGCAAAATGGATa                                                                           >  1:829729/1‑148 (MQ=255)
                                                                             aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGCGGGGGTTAATGGCCAATAAGCGACG‑GTCA‑TCCCACGACGCCGCGGTTAACGCAAGCGCG‑GGTGGTTCCGCGAAATGGAAAgc                                                                         >  1:1120551/1‑150 (MQ=255)
                                                                                    ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAATGCCAGATAAGCGACG‑GTCA‑TCCCCGGACGCAGCGTTAAAAGCAAGCGCG‑GGGCGTCACGCCAAAAGGATAGcttgctgg                                                                  >  1:1100066/1‑150 (MQ=255)
                                                                                     gTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTAAAAGCCAGAAAGGGGAAG‑TTCA‑TCCCCGGACGCACAGGTAAAAGGAAGCGCG‑GGGCGGTACGCAAAATGGGTa                                                                           >  2:904181/1‑140 (MQ=255)
                                                                                     gTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGgca                                                                <  1:1569251/151‑1 (MQ=255)
                                                                                       tattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGgcag                                                               <  2:109569/150‑1 (MQ=255)
                                                                                        attaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGgcagc                                                              <  1:1059635/150‑1 (MQ=255)
                                                                                           aaTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTTAAAGGCAGAAAAGCGGCG‑TTTA‑TCCCCAGGCGCAGCGGTAAACGCCAAGGCG‑GGGCGTCACGCAAAAAGGATAggtgg                                                                      >  2:699010/1‑139 (MQ=255)
                                                                                             tCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTAAAAGACAGAAAAGGGAAG‑GTTA‑ACCCAGGGCGCCGCGCTAAACGCAAACGCG‑GGGCGTCCCGCAAAAAgggaaagtgggcgg                                                                  >  2:1460238/1‑129 (MQ=255)
                                                                                             tCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGCGGGGGTAAAACACCGAGAGGCGCAG‑TTCC‑TCCCAGGACCGAGCGTTAAACACAAACGCC‑GGTAGCCCCGCAAAAAAGTTAGGCTGCTGGcagcagc                                                           >  2:914859/1‑147 (MQ=255)
                                                                                              cGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggccag                                                       <  1:1150396/151‑1 (MQ=255)
                                                                                               gCCACTGCCAGCGCAGGGCGGACATCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAACAGCCAGATAGGCGACG‑TTCA‑TACCAGGACGCAGCGTTGAACGCATGCGCA‑GTGCGTCGCGTACAATGGATAGCTGTCTCGCAGcaggcca                                                        <  2:205914/149‑1 (MQ=255)
                                                                                                   gTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GTTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGgcaaggacag                                                  <  2:569722/151‑1 (MQ=255)
                                                                                                        aGCGCATGGCGGACATCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGGTAGCTGGCTGGGAGCAGGCCAGGCCAGAACGt                                             <  2:81062/151‑1 (MQ=255)
                                                                                                         gcgcATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTc                                            <  1:1037457/151‑1 (MQ=255)
                                                                                                         gcgcATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTc                                            <  2:1101555/151‑1 (MQ=255)
                                                                                                            cATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCa                                         <  1:290067/151‑1 (MQ=255)
                                                                                                             aTGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGTGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCAt                                        <  2:367438/151‑1 (MQ=255)
                                                                                                               ggCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTc                                      <  1:383858/151‑1 (MQ=255)
                                                                                                                     cAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGt                                <  1:205433/151‑1 (MQ=255)
                                                                                                                       gCGGGATTATTTAGCGGCGGTTTGGTGGCGTTATATACCCAATAGGGGTCC‑TTCA‑TTCCCGGACGGAGCGGTAAACGCAAGAGCGCGGTGGGC‑CCCAAAATGGGGAGCTGGCTGGCAGCAGGCCAGGGCCGAACgcccgccccccccggc                               >  2:103687/1‑135 (MQ=255)
                                                                                                                        cGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAg                             <  1:1377602/151‑1 (MQ=255)
                                                                                                                            attattTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCGCGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGt                         <  1:704612/151‑1 (MQ=255)
                                                                                                                            attattTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGt                         <  1:440225/151‑1 (MQ=255)
                                                                                                                            attattTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGt                         <  2:908766/151‑1 (MQ=255)
                                                                                                                            attattTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGt                         <  1:1327699/151‑1 (MQ=255)
                                                                                                                               attTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGttt                       <  2:919697/150‑1 (MQ=255)
                                                                                                                                 ttAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCg                     <  1:210264/150‑1 (MQ=255)
                                                                                                                                  tAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGa                    <  1:292721/150‑1 (MQ=255)
                                                                                                                                   aGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACa                  <  1:20260/151‑1 (MQ=255)
                                                                                                                                    gcggcgGTTTGGCGGGGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGGCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGAc                   >  1:182330/1‑149 (MQ=255)
                                                                                                                                     tggcggTCTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGAc                   >  2:566588/2‑148 (MQ=255)
                                                                                                                                     cggcggTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGa                    >  2:1519765/1‑147 (MQ=255)
                                                                                                                                         ggTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTc            >  2:1327970/1‑151 (MQ=255)
                                                                                                                                          gTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTcc           <  1:1506727/151‑1 (MQ=255)
                                                                                                                                             tGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑GCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCgg         <  2:873292/150‑1 (MQ=255)
                                                                                                                                              ggCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGt        <  1:1633604/150‑1 (MQ=255)
                                                                                                                                              ggCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTg       <  1:908705/151‑1 (MQ=255)
                                                                                                                                                 ggTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccg    >  1:1228127/1‑151 (MQ=255)
                                                                                                                                                  gtgtTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgc   >  2:257179/1‑151 (MQ=255)
                                                                                                                                                   tgtTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgcc  >  1:362795/1‑151 (MQ=255)
                                                                                                                                                   tgtTAAATGCCAGATAGGCGACG‑GTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccg    <  2:1248764/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                         
CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACG‑TTCA‑TCCCAGGACGCAGCGTTAAACGCAAGCGCG‑GGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCC  >  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: ‑