Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,356,2010TG54.5% 54.6 / 89.7 66T304P (ACC→CCC) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (8/22);  new base G (36/0);  total (44/22)
Fisher's exact test for biased strand distribution p-value = 3.21e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.09e-07

CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCG  >  NC_000913/1356054‑1356349
                                                                                                                                                      |                                                                                                                                                      
cGTTGTCATCGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTTC‑GGGGtt                                                                                                                                                     >  1:412870/1‑151 (MQ=255)
    gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGGGTTTAAt                                                                                                                                                 >  1:1523245/1‑151 (MQ=255)
     tCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAAt                                                                                                                                                 <  2:578815/150‑1 (MQ=255)
      cATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGc                                                                                                                                               <  1:1080362/151‑1 (MQ=255)
      cATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGc                                                                                                                                               <  1:452150/151‑1 (MQ=255)
            cccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTGC‑GGGGTTTAAAGCCAGATa                                                                                                                                         >  1:1175250/1‑151 (MQ=255)
            cccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGCCAGATa                                                                                                                                         <  1:1047380/151‑1 (MQ=255)
              cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTGC‑GGGGGTTAATGCCAGATAgg                                                                                                                                       >  1:381947/1‑151 (MQ=255)
                         cgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTGC‑GGTGGTTAATGCCAGAAAGGGGAC‑GTTCAt                                                                                                                            >  2:1486418/1‑150 (MQ=255)
                           gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTGC‑GGGGGTTAAAGCCAGAAAAGCGAC‑ggttattc                                                                                                                          >  2:1502619/1‑143 (MQ=255)
                           gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATcc                                                                                                                          <  1:938504/150‑1 (MQ=255)
                           gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑CTTTGG‑GGTGGTTAATGCAAGATAGGCAAC‑GCTAATcg                                                                                                                          >  2:1615000/1‑149 (MQ=255)
                            tAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATccc                                                                                                                         <  2:914587/150‑1 (MQ=255)
                            tAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATccc                                                                                                                         <  1:1000587/150‑1 (MQ=255)
                              aaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGG‑TTTGGG‑GGGGTAAAAAGCCAAAAAggggcgtttaacccccg                                                                                                                        >  1:1242124/1‑136 (MQ=255)
                                aTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGAc                                                                                                                    <  2:712854/151‑1 (MQ=255)
                                         gccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTAGC‑CGGGGTAAATGACCAAAAGTGGAC‑GTTTTACCCAGGACGGAAC‑GTTa                                                                                                          >  2:1383030/1‑151 (MQ=255)
                                            gTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTGC‑GGGGGTAAAAGCCAGAAAAGCGAC‑GTTCATCCCAGGAGGCAGC‑GTTtaaa                                                                                                       >  1:469610/1‑150 (MQ=255)
                                             tttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACg                                                                                                      <  2:1259556/151‑1 (MQ=255)
                                             tttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACg                                                                                                      <  1:33648/151‑1 (MQ=255)
                                                cAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTGC‑GGGGTTTAATGCCAGAAAGGCGGC‑GTTCCTCCCCGGGCGCAGC‑GTTAAACGCaa                                                                                                   >  2:528050/1‑151 (MQ=255)
                                                  gCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGG‑GGGGTTAAAAGCCAAATAGGCGAC‑GGTCATCCCCAGACGCAGC‑GGTAAAAGCAAg                                                                                                  >  1:82351/1‑150 (MQ=255)
                                                   ccGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGGGGTAAAAGCCACATAGGGGGC‑GTTTCTCCCCGGACGCACC‑GTTTAACGCAAgcg                                                                                                >  2:132919/1‑151 (MQ=255)
                                                     gTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGG‑TTTGGG‑GGGGTTAAAT‑CCCGAAAGGGGGC‑GTTTATCCCCAGAAGCAGCGGTTAAAC‑CAAGcgcg                                                                                              >  2:1317229/1‑151 (MQ=255)
                                                      tACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTGC‑GGGGTTAAAAGCCCGATAGGGGAC‑GTTTATCCCAGGACGCAGC‑GTTTAACGCAAGCGCgg                                                                                             >  1:991005/1‑151 (MQ=255)
                                                       aCCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGG‑GGGGGTAAATTGCAGAAAGGGGGC‑GGTCATTCCCGGACGCCGCGGTTAAACG‑AAACGCggg                                                                                            >  2:309299/1‑151 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTGG‑GGGGTTTAATGCCCGATAAGGGAC‑GTTTATTCCAGGACGCCGC‑GGTTAACCGAAGCGCgggg                                                                                           >  2:1220636/1‑150 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTGG‑GGGGGTAAATGCCAGAAAGGGGACGGTTCTTCCC‑GGGCGCCGC‑GTTTAACGCAAACGCgggg                                                                                           >  1:1610588/1‑150 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGG‑GGGGTTAAATGCCAGATAGGCGAC‑GTTCATCCCCGGACGCAGC‑GGTAAAAGCAAACGCgggg                                                                                           >  1:1100170/1‑150 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGG‑GGGGTTAAATGCCAAATAGGGGGC‑GGTCATCCCCGGAAGCAGC‑GGTAAAAGCAAGCGCgggg                                                                                           >  1:1221373/1‑150 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGG‑GGGGGTAAAAGCCAGAAAGGCGAC‑GGTCATCCCAGGGCGCAGC‑GGTAAACGCCAGCGCgggg                                                                                           >  2:1429346/1‑150 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGCGGGGGTTAAATG‑CAGAAAAGGGGC‑GTTTATCCCCGGACGCCGC‑GGTAAAAGccaaagccggg                                                                                           >  2:40092/1‑141 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGGGTTTAATGCCAAATAAGGGGC‑GTTTATCCCCGGGCCCCGC‑GTTAAACGCAAGAGCcggg                                                                                           >  1:9751/1‑150 (MQ=255)
                                                          ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTGG‑GGGGGTAAATGCCAGATAGGCGAC‑GGTTATCCCAGGACGCAGC‑GTTAAAAGCAAACGCGGggcg                                                                                         >  2:897132/1‑151 (MQ=255)
                                                          ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGG‑GGGGGTAAAAGCCAGAAAGGGGAC‑GGTCACCCCCCGGCGCAGC‑GGTAAAACCCAACGCgggggg                                                                                         >  1:1453992/1‑148 (MQ=255)
                                                          ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGGGTTAAATGGCAGAAAAGCGAC‑GGTCATCCCCGGAAGCAGC‑GGTAAACGCAAGCGCGGggcg                                                                                         >  1:210107/1‑151 (MQ=255)
                                                          ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGGGGTTAAAGCCAGATAGGGGGC‑GGTCCTCCCCCGGCGCCGC‑GGTTAACCCCAGCGCGGggcg                                                                                         >  1:1081699/1‑151 (MQ=255)
                                                           taataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTGC‑GGTGGTTAAAGCCAGCTAGGGGAC‑GTTCCTCCCCGGACGCAGC‑GTTAAAAGCAAACGCGGggcgt                                                                                        >  2:1267565/1‑151 (MQ=255)
                                                              taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑GTTGGG‑GGGGTTAAAATGCAGATAAGCGGC‑GCTCCTTCCAGGGCCCAAC‑GCTTAAAACCAGGGCGGggcgtcc                                                                                      >  1:1359818/1‑149 (MQ=255)
                                                                   ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTGC‑GGTGTTAAAAGCCCGAAAAGCGAC‑GTTCCTCCCCGGGCGCCGC‑GGTTAAAGCAAACGCGGGGCGTCACGCaaa                                                                                >  1:70873/1‑151 (MQ=255)
                                                                        gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGc                                                                                                                                               <  1:1105383/85‑1 (MQ=255)
                                                                        gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑GTTAGC‑GGGGTTAAAAGc                                                                                                                                               >  2:1104861/1‑85 (MQ=255)
                                                                          aTCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGG‑GGGGGTCAACACCCGATAGGGGAC‑GCTCCTCCCACGAGGCAGA‑GTTAAAAACCAGGGCGGGGCGGCCCCCAAAAGGGGTa                                                                         >  1:588325/1‑151 (MQ=255)
                                                                                cgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctgg                                                                    <  1:1501807/150‑1 (MQ=255)
                                                                                    ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTGG‑GGGGTTAAAAGCCAGAAAGGGGAC‑GGTCATCCCCGGAAGCAGG‑GGTAAAAGCAAGCGCGGGGCGGCACGCAAAAAGGATAGctggcgggg                                                               >  1:1072350/1‑150 (MQ=255)
                                                                                    ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTGG‑GGGGGTAAAAGCCAGAAAGGCGAC‑GTTTATCCCAGGACGCAGC‑GTTAAAAGCAAACGCGGGGCGGCACGCAAAAAGGAAAgccggcgtgg                                                               >  2:931282/1‑147 (MQ=255)
                                                                                    ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGG‑GGGGTTAAAAGCCAGATAGGGGAC‑GTTTATCCCCGGAAGCAGC‑GTTTAAAGCCAGCGCGGGGCGTCACGCAAAAAGGATAgccgggggg                                                                >  1:897108/1‑146 (MQ=255)
                                                                                    ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGG‑GGGGGTTAATGGCAGATAGGCGAC‑GTTCATCCCCCGGCGCAGC‑GCTAAACACAAGCGCGGGTGGCCCCCCGAAATGGATAGccggcggg                                                                >  2:201255/1‑150 (MQ=255)
                                                                                    ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑GTTGGC‑GGGGTTAAATGCCAGATAGGGGAC‑GTTTATCCCCGGAAGGAGG‑GTTAAACGCCAGCCCGGGTCGTCACGCAAAATGGGAAGctggcggg                                                                >  2:1494937/1‑150 (MQ=255)
                                                                                     gTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGGGTTAAATGCCAGAAAGGCGAC‑GGTCCTTCCAGGAAGCCGC‑GTTTAACGCAAACGGGGGGGGGCAAGCAAAAAGGATAGatggcgggga                                                              >  1:1242705/1‑149 (MQ=255)
                                                                                     gTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑CTTGGC‑GGTGTTTAAATCCAGATAGGCGGC‑GGGCGTCCCAAGATGCCGG‑GCTTAAAGCAAGCGCGGGGCCTCACCCAAAAAGGGTAgtgcggttgga                                                              >  1:248485/1‑141 (MQ=255)
                                                                                       tattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcagc                                                            <  2:1483307/151‑1 (MQ=255)
                                                                                            aTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggcc                                                       <  2:128904/151‑1 (MQ=255)
                                                                                                   gTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGG‑GGGGGTAAATGCCAGATAGGCGGC‑GTTCATCCCCGGACGCCGC‑GGTTAAAGGAAGCGCGGGGGGGCCCGCAAAAAGGATAGCTGGCTGGGAGCAGgccaggccag                                                >  1:1576809/1‑151 (MQ=255)
                                                                                                    ttccAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAAGGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                  <  2:507359/146‑1 (MQ=255)
                                                                                                    tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                  <  2:562966/148‑1 (MQ=255)
                                                                                                    tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                  <  2:1167482/148‑1 (MQ=255)
                                                                                                    tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                  <  2:428639/148‑1 (MQ=255)
                                                                                                         gcgcATGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTTGGG‑GGGGTTAAAT‑CCAGATAGGGGGC‑GGTCATCCCCGGGCGCCGG‑GTTAAAAGCAAACGCGGGGCGGCACCCAAAAAGGGAAGCTGGCCGGGAGCaggcccggccccgtct                                            >  1:1295385/1‑143 (MQ=255)
                                                                                                         gcgcATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                  <  2:1294778/143‑1 (MQ=255)
                                                                                                                gCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTTAATGCCCGATAGGCGGC‑CGTCTTTCCCCGGCCCAGC‑GGTTAACGCCAGCGCGGGGCGTCCCACCAAAAGGAAAGCTGGGTGGCCGCccggcaggcccggaaaccgccttt                                    >  1:1370691/1‑136 (MQ=255)
                                                                                                                                       gc‑gg‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGAc             <  1:201306/150‑1 (MQ=255)
                                                                                                                                        c‑gg‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACg            <  1:1527518/150‑1 (MQ=255)
                                                                                                                                           g‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGGTTCGACAGACGTc          >  2:1367714/1‑150 (MQ=255)
                                                                                                                                              ttGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCTCATTCCCCGGTCAGGAGTTTCGACAGACGTCCgg       <  1:986313/151‑1 (MQ=255)
                                                                                                                                                ggC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTg     >  1:1362271/1‑151 (MQ=255)
                                                                                                                                                     gtgtTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccg  >  1:956851/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCG  >  NC_000913/1356054‑1356349

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