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

CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCG  >  NC_000913/1356054‑1356349
                                                                                                                                                      |                                                                                                                                                      
cGTTGTCATCGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTTCGGGGtt                                                                                                                                                     >  1:412870/1‑151 (MQ=255)
    gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAAt                                                                                                                                                 >  1:1523245/1‑151 (MQ=255)
     tCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAAt                                                                                                                                                 <  2:578815/150‑1 (MQ=255)
      cATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGc                                                                                                                                               <  1:1080362/151‑1 (MQ=255)
      cATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGc                                                                                                                                               <  1:452150/151‑1 (MQ=255)
            cccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTTAAAGCCAGATa                                                                                                                                         >  1:1175250/1‑151 (MQ=255)
            cccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATa                                                                                                                                         <  1:1047380/151‑1 (MQ=255)
              cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAATGCCAGATAgg                                                                                                                                       >  1:381947/1‑151 (MQ=255)
                         cgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGTGGTTAATGCCAGAAAGGGGACGTTCAt                                                                                                                            >  2:1486418/1‑150 (MQ=255)
                           gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTTAAAGCCAGAAAAGCGACggttattc                                                                                                                          >  2:1502619/1‑143 (MQ=255)
                           gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATcc                                                                                                                          <  1:938504/150‑1 (MQ=255)
                           gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGCTTTGGGGTGGTTAATGCAAGATAGGCAACGCTAATcg                                                                                                                          >  2:1615000/1‑149 (MQ=255)
                            tAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATccc                                                                                                                         <  2:914587/150‑1 (MQ=255)
                            tAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATccc                                                                                                                         <  1:1000587/150‑1 (MQ=255)
                              aaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGGGGGGGTAAAAAGCCAAAAAggggcgtttaacccccg                                                                                                                        >  1:1242124/1‑136 (MQ=255)
                                aTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGAc                                                                                                                    <  2:712854/151‑1 (MQ=255)
                                         gccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTAGCCGGGGTAAATGACCAAAAGTGGACGTTTTACCCAGGACGGAACGTTa                                                                                                          >  2:1383030/1‑151 (MQ=255)
                                            gTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGGTAAAAGCCAGAAAAGCGACGTTCATCCCAGGAGGCAGCGTTtaaa                                                                                                       >  1:469610/1‑150 (MQ=255)
                                             tttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACg                                                                                                      <  2:1259556/151‑1 (MQ=255)
                                             tttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACg                                                                                                      <  1:33648/151‑1 (MQ=255)
                                                cAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTTAATGCCAGAAAGGCGGCGTTCCTCCCCGGGCGCAGCGTTAAACGCaa                                                                                                   >  2:528050/1‑151 (MQ=255)
                                                  gCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTAAAAGCCAAATAGGCGACGGTCATCCCCAGACGCAGCGGTAAAAGCAAg                                                                                                  >  1:82351/1‑150 (MQ=255)
                                                   ccGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTAAAAGCCACATAGGGGGCGTTTCTCCCCGGACGCACCGTTTAACGCAAgcg                                                                                                >  2:132919/1‑151 (MQ=255)
                                                     gTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGGGGGGGTTAAATCCCGAAAGGGGGCGTTTATCCCCAGAAGCAGCGGTTAAACCAAGcgcg                                                                                              >  2:1317229/1‑151 (MQ=255)
                                                      tACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGGGTTAAAAGCCCGATAGGGGACGTTTATCCCAGGACGCAGCGTTTAACGCAAGCGCgg                                                                                             >  1:991005/1‑151 (MQ=255)
                                                       aCCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTAAATTGCAGAAAGGGGGCGGTCATTCCCGGACGCCGCGGTTAAACGAAACGCggg                                                                                            >  2:309299/1‑151 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTTAATGCCCGATAAGGGACGTTTATTCCAGGACGCCGCGGTTAACCGAAGCGCgggg                                                                                           >  2:1220636/1‑150 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTAAATGCCAGAAAGGGGACGGTTCTTCCCGGGCGCCGCGTTTAACGCAAACGCgggg                                                                                           >  1:1610588/1‑150 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTAAATGCCAGATAGGCGACGTTCATCCCCGGACGCAGCGGTAAAAGCAAACGCgggg                                                                                           >  1:1100170/1‑150 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTAAATGCCAAATAGGGGGCGGTCATCCCCGGAAGCAGCGGTAAAAGCAAGCGCgggg                                                                                           >  1:1221373/1‑150 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTAAAAGCCAGAAAGGCGACGGTCATCCCAGGGCGCAGCGGTAAACGCCAGCGCgggg                                                                                           >  2:1429346/1‑150 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTTAAATGCAGAAAAGGGGCGTTTATCCCCGGACGCCGCGGTAAAAGccaaagccggg                                                                                           >  2:40092/1‑141 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTTAATGCCAAATAAGGGGCGTTTATCCCCGGGCCCCGCGTTAAACGCAAGAGCcggg                                                                                           >  1:9751/1‑150 (MQ=255)
                                                          ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTAAATGCCAGATAGGCGACGGTTATCCCAGGACGCAGCGTTAAAAGCAAACGCGGggcg                                                                                         >  2:897132/1‑151 (MQ=255)
                                                          ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTAAAAGCCAGAAAGGGGACGGTCACCCCCCGGCGCAGCGGTAAAACCCAACGCgggggg                                                                                         >  1:1453992/1‑148 (MQ=255)
                                                          ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATGGCAGAAAAGCGACGGTCATCCCCGGAAGCAGCGGTAAACGCAAGCGCGGggcg                                                                                         >  1:210107/1‑151 (MQ=255)
                                                          ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGGTTAAAGCCAGATAGGGGGCGGTCCTCCCCCGGCGCCGCGGTTAACCCCAGCGCGGggcg                                                                                         >  1:1081699/1‑151 (MQ=255)
                                                           taataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGTGGTTAAAGCCAGCTAGGGGACGTTCCTCCCCGGACGCAGCGTTAAAAGCAAACGCGGggcgt                                                                                        >  2:1267565/1‑151 (MQ=255)
                                                              taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTGGGGGGGTTAAAATGCAGATAAGCGGCGCTCCTTCCAGGGCCCAACGCTTAAAACCAGGGCGGggcgtcc                                                                                      >  1:1359818/1‑149 (MQ=255)
                                                                   ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGCGGTGTTAAAAGCCCGAAAAGCGACGTTCCTCCCCGGGCGCCGCGGTTAAAGCAAACGCGGGGCGTCACGCaaa                                                                                >  1:70873/1‑151 (MQ=255)
                                                                        gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGc                                                                                                                                               <  1:1105383/85‑1 (MQ=255)
                                                                        gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTAGCGGGGTTAAAAGc                                                                                                                                               >  2:1104861/1‑85 (MQ=255)
                                                                          aTCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTCAACACCCGATAGGGGACGCTCCTCCCACGAGGCAGAGTTAAAAACCAGGGCGGGGCGGCCCCCAAAAGGGGTa                                                                         >  1:588325/1‑151 (MQ=255)
                                                                                cgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctgg                                                                    <  1:1501807/150‑1 (MQ=255)
                                                                                    ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGTTAAAAGCCAGAAAGGGGACGGTCATCCCCGGAAGCAGGGGTAAAAGCAAGCGCGGGGCGGCACGCAAAAAGGATAGctggcgggg                                                               >  1:1072350/1‑150 (MQ=255)
                                                                                    ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTAAAAGCCAGAAAGGCGACGTTTATCCCAGGACGCAGCGTTAAAAGCAAACGCGGGGCGGCACGCAAAAAGGAAAgccggcgtgg                                                               >  2:931282/1‑147 (MQ=255)
                                                                                    ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGTTAAAAGCCAGATAGGGGACGTTTATCCCCGGAAGCAGCGTTTAAAGCCAGCGCGGGGCGTCACGCAAAAAGGATAgccgggggg                                                                >  1:897108/1‑146 (MQ=255)
                                                                                    ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTTAATGGCAGATAGGCGACGTTCATCCCCCGGCGCAGCGCTAAACACAAGCGCGGGTGGCCCCCCGAAATGGATAGccggcggg                                                                >  2:201255/1‑150 (MQ=255)
                                                                                    ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTGGCGGGGTTAAATGCCAGATAGGGGACGTTTATCCCCGGAAGGAGGGTTAAACGCCAGCCCGGGTCGTCACGCAAAATGGGAAGctggcggg                                                                >  2:1494937/1‑150 (MQ=255)
                                                                                     gTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGGGTTAAATGCCAGAAAGGCGACGGTCCTTCCAGGAAGCCGCGTTTAACGCAAACGGGGGGGGGCAAGCAAAAAGGATAGatggcgggga                                                              >  1:1242705/1‑149 (MQ=255)
                                                                                     gTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGCTTGGCGGTGTTTAAATCCAGATAGGCGGCGGGCGTCCCAAGATGCCGGGCTTAAAGCAAGCGCGGGGCCTCACCCAAAAAGGGTAgtgcggttgga                                                              >  1:248485/1‑141 (MQ=255)
                                                                                       tattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcagc                                                            <  2:1483307/151‑1 (MQ=255)
                                                                                            aTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGcaggcc                                                       <  2:128904/151‑1 (MQ=255)
                                                                                                   gTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGGGGGGGTAAATGCCAGATAGGCGGCGTTCATCCCCGGACGCCGCGGTTAAAGGAAGCGCGGGGGGGCCCGCAAAAAGGATAGCTGGCTGGGAGCAGgccaggccag                                                >  1:1576809/1‑151 (MQ=255)
                                                                                                    ttccAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAAGGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                  <  2:507359/146‑1 (MQ=255)
                                                                                                    tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                  <  2:562966/148‑1 (MQ=255)
                                                                                                    tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                  <  2:1167482/148‑1 (MQ=255)
                                                                                                    tGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                  <  2:428639/148‑1 (MQ=255)
                                                                                                         gcgcATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGGGGGGTTAAATCCAGATAGGGGGCGGTCATCCCCGGGCGCCGGGTTAAAAGCAAACGCGGGGCGGCACCCAAAAAGGGAAGCTGGCCGGGAGCaggcccggccccgtct                                            >  1:1295385/1‑143 (MQ=255)
                                                                                                         gcgcATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCaggccaggcc                                                  <  2:1294778/143‑1 (MQ=255)
                                                                                                                gCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTTAATGCCCGATAGGCGGCCGTCTTTCCCCGGCCCAGCGGTTAACGCCAGCGCGGGGCGTCCCACCAAAAGGAAAGCTGGGTGGCCGCccggcaggcccggaaaccgccttt                                    >  1:1370691/1‑136 (MQ=255)
                                                                                                                                       gcggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGAc             <  1:201306/150‑1 (MQ=255)
                                                                                                                                        cggTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACg            <  1:1527518/150‑1 (MQ=255)
                                                                                                                                           gTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGGTTCGACAGACGTc          >  2:1367714/1‑150 (MQ=255)
                                                                                                                                              ttGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCTCATTCCCCGGTCAGGAGTTTCGACAGACGTCCgg       <  1:986313/151‑1 (MQ=255)
                                                                                                                                                ggCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTg     >  1:1362271/1‑151 (MQ=255)
                                                                                                                                                     gtgtTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccg  >  1:956851/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGCAGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCG  >  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: