Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913573,4090CA41.3% 166.7 / 39.3 93S13R (AGC→AGA) rusADLP12 prophage; endonuclease RUS
Reads supporting (aligned to +/- strand):  ref base C (48/6);  new base A (1/37);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 2.91e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCAC  >  NC_000913/573259‑573556
                                                                                                                                                           |                                                                                                                                                     
acgaAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGc                                                                                                                                                       <  4:1097798/151‑1 (MQ=255)
 tgaAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGc                                                                                                                                                       >  1:26319/2‑150 (MQ=255)
    aaTCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATAACTACAGCATCACATTACCCTGG‑CCTCC‑GAGCaat                                                                                                                                                    >  2:1508955/1‑150 (MQ=255)
     aTCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCaata                                                                                                                                                   >  1:202576/1‑150 (MQ=255)
      tCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCaataat                                                                                                                                                 >  3:423974/1‑151 (MQ=255)
        gACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATc                                                                                                                                                >  1:415469/1‑150 (MQ=255)
        gACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCg                                                                                                                                               >  2:50590/1‑151 (MQ=255)
        gACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCg                                                                                                                                               >  1:249255/1‑151 (MQ=255)
          ccgccgCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCt                                                                                                                                             >  2:1231377/1‑151 (MQ=255)
           cgccgcACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCt                                                                                                                                             >  2:890964/1‑150 (MQ=255)
            gccgcACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑t                                                                                                                                          >  1:222306/1‑151 (MQ=255)
             ccgcACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑tt                                                                                                                                         <  1:144887/151‑1 (MQ=255)
             ccgcACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑tt                                                                                                                                         >  4:708301/1‑151 (MQ=255)
                 acgcaTTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCg                                                                                                                                     >  4:1585237/1‑151 (MQ=255)
                    caTTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGcc                                                                                                                                   >  2:1224058/1‑150 (MQ=255)
                         tgttgACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAAt                                                                                                                              >  2:1392741/1‑150 (MQ=255)
                           ttgACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATcg                                                                                                                            >  3:44141/1‑150 (MQ=255)
                             gACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCgg                                                                                                                         >  2:1552012/1‑151 (MQ=255)
                              aCGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCggg                                                                                                                        >  2:840080/1‑151 (MQ=255)
                                gtggaatataaaaaaGAATGCGCGCTGGAAGGTAGGGAAAAAACACAGGTTTTCGGGCTGAAAGGGGGGT‑TTATT‑AG‑GGCGTGAATACCTACAGCACCACATTCCCCGGG‑CCCCG‑GAGAAATAATCTCTT‑TTCCCGCCATAATCGCgg                                                                                                                         <  2:1163713/137‑1 (MQ=255)
                                gCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGgcg                                                                                                                      >  3:164430/1‑151 (MQ=255)
                                   ggATATAAAAAAGAAGGCGGGTGGGAAGTTAGGGGGAAAACACAGGTT‑TTTGGGTTAAAGGGGGGGTTTATT‑AG‑GGCGTAAATCCCTACAGCATCACATTCCCCGGC‑CCCCC‑GAGAAATAACCGCTA‑TTACCGCCATAATCGCGGgcgc                                                                                                                     <  4:758437/149‑1 (MQ=255)
                                   gattaaaaaaaaGAAGGCGCGGGAAAAGGTAGGGCGAAAACACGGGTTATCGGGAAGAAAGAGGGGT‑TTATA‑AA‑GGCGTAAAT‑CCCACCACCCCCCATTCCCCCGGCCCCCC‑GACCAATAACCCCTA‑TTACCGCCATAATCGCgg                                                                                                                         <  2:107098/137‑1 (MQ=255)
                                    gATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGcgca‑c                                                                                                                  >  3:269080/1‑150 (MQ=255)
                                        taaaaaaGAAGGCGCGCGGGAGGGTATGGGGAAAACACGGTTTATCTGGCTAAAAGGGGGGT‑TTATT‑AA‑GGGGTGAATCCCTACAGCATCACATTCCCCTGG‑CCCCC‑GAGAAATAATCGTTT‑TTACCGCCATAATCGCGGGcgca‑cgc                                                                                                                <  3:222321/145‑1 (MQ=255)
                                        tGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGca‑cgcacg                                                                                                             >  1:415152/1‑151 (MQ=255)
                                        tGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGca‑cgcacg                                                                                                             >  4:310755/1‑151 (MQ=255)
                                         gCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGt                                                                                                            >  4:1437146/1‑151 (MQ=255)
                                          aaaaaaaAGGCGCGCGGGAAGGTGGGGTGAAAACACTGTTTTTCGGGTGGAAAGGGGGGT‑TTTTT‑AA‑GGAGGAAATACTTCCAGCACCACATTCCCCTGT‑CCCCA‑AAGAAATACACGTTA‑TTCCCGCCATAATCGCGGGcgca‑cgc                                                                                                                <  2:1310678/145‑1 (MQ=255)
                                          aaaaaGAAGGCGGGCGGGAAGGTAGGGGGAGAACCCAGTTTATCTGGCGAAAAGGGGGGT‑TTATT‑AG‑GGCGTAAATCCCTACAGCATCACATCCCCCTGC‑CCCCG‑GAGAAATATCCGTTT‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTc                                                                                                           <  3:97821/150‑1 (MQ=255)
                                           aaaaGAATGCGCGTTGGAAGGTATGGCGAGAACACGGTTTATGTGGCTGAAGGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACACCATCACATTACCCTGG‑CCTCC‑GACCAATAATCGCTA‑TTACCGCCATAATCGCGGGcgca‑cg                                                                                                                 <  2:1573165/144‑1 (MQ=255)
                                           aaaaGAATGCGCGGGGAAAGGTAGGGGGAAAACACAGTTTATCTGGCTGAAGGGGGGGT‑TTATA‑AG‑GGCGTAAATCCCTCCACCACCACATCCCCCTGC‑CCTCC‑GAGCAAAAATCTTTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCa                                                                                                          <  1:442836/151‑1 (MQ=255)
                                           aaaaGAAGGCGGGCGGAAAGGAAGGGAAAAAACCCAGTTTATCTGG‑GAAAAGGGGGGG‑TTTTTGAA‑GGGGTAAATCCCCACAGCATCACATTCCCCGGC‑CCACG‑GGAAAATAATCGATA‑TTGCCGCCATAATCGCGGGcgca‑cgca                                                                                                               <  2:840117/146‑1 (MQ=255)
                                            aaaGAATGCGCGCTGGAAGGTAGGGCGAGAACCGAGGTTGTCGGGAAGAAAGGGGGGT‑TTATT‑AG‑GGCGTTAATCCCTTCAG‑ATTATCCTCCCCCCG‑CCGCCAGAGAAAAAATGGCTT‑TCACCGCCATAATCGCGGGCGCA‑CGCACGTCAg                                                                                                         <  2:1018410/151‑1 (MQ=255)
                                                aaTGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAgcgc                                                                                                      >  2:1447652/1‑150 (MQ=255)
                                                        cgggAAGGTAGGGCGAAAACCGAGGTTATGGGGAAAAAAGGGGGGT‑TTATA‑AG‑GGAGAAAATTCTTACAGCAACTCTTCCCCTTGC‑CCTCC‑GAGAAATAATCGTTT‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGggc                                                                                              <  3:202590/148‑1 (MQ=255)
                                                           gAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGggcaggc                                                                                          >  3:138716/1‑151 (MQ=255)
                                                             aGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCAt                                                                                        >  4:1647265/1‑151 (MQ=255)
                                                              ggTATGGCGAAAACACGGTTTATCTGGCTGAAAGGGGGGG‑TTATA‑AA‑GGCGTAAATCCCTCCACCATCACATTCCCCTGC‑CCTCA‑GACAAAAAACCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGggcaggc                                                                                          <  4:1593113/148‑1 (MQ=255)
                                                                tATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GTGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATAcc                                                                                     >  1:469820/1‑151 (MQ=255)
                                                                tATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATAcc                                                                                     >  1:272382/1‑151 (MQ=255)
                                                                   ggggagaACACAGTTTATGTGGCTGAAAGGGGGGT‑TTATT‑AA‑GGGGGAAATCCCTACAGCACCACTTCCCCCGGC‑CCTCC‑GACCAATAATCGTTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATa                                                                                       <  2:1512090/143‑1 (MQ=255)
                                                                   ggagaaaccacaGTTTATCGGGTGGAAAGGGGGGT‑TTATT‑AA‑GGCGTAAATCCCTACAGCACCACATTCCCCGGC‑CCCCC‑GAAAAAAAATCCTTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATAc                                                                                      <  4:630366/139‑1 (MQ=255)
                                                                     cGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGAt                                                                                >  2:720328/1‑151 (MQ=255)
                                                                          acacAGGTTATCGGGCAAAAGGGGGGGG‑TTATA‑AG‑GGCGTAAACCCCCACACCACCACATCCCCCGGC‑CCCCC‑GAGAAAAAACCTCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAAc                                                                             <  2:150406/149‑1 (MQ=255)
                                                                            acaGGTTATCTGTATGAAAGGGGGTT‑TTATA‑AA‑GGCGTGAATCCCCACAACATCCCTTTCCCCTGC‑CCCCG‑GAGAAATACTCGCTT‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATa                                                                               <  2:1399356/145‑1 (MQ=255)
                                                                             caGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGc                                                                        >  2:1162594/1‑151 (MQ=255)
                                                                                 ttctgcggCTGAAAGGGGGTG‑TTTAT‑AA‑GGTGTGAATCCCGAATCTATCACTTCCCCGTCT‑CCACA‑AAGCAATAAAAGGTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATCCCGCGATAACGTc                                                                          <  1:144578/139‑1 (MQ=255)
                                                                                 ttatttggggaaaaGGGGGGG‑TTAAT‑AG‑GCGGTAAATCCCCACACCACCACTTCCCCGTGG‑CCCCC‑GACAAATATGCCCTA‑TCGCCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGtt                                                                          <  2:1595049/134‑2 (MQ=255)
                                                                                 ttATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGccc                                                                      >  2:1620642/1‑149 (MQ=255)
                                                                                 ttATCTGGCTAAAAGAGGGGG‑TTATT‑AG‑GGCGTAAATCCCTACAGCATCACATTCCCCTGG‑CCCCG‑GAGAAATAATCGTTT‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGccc                                                                      <  3:26321/149‑1 (MQ=255)
                                                                                  tATCTGGCTGATAGCGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGTGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGa                                                                    >  2:377838/1‑150 (MQ=255)
                                                                                   aTCGGGCTAAAAG‑GGGGG‑TTTTT‑AAGGGCGAGAAACCCCACAGCACCACTTCCCCCGGG‑CCCCG‑GACAAATAGCCGTTT‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAAt                                                                  <  2:950433/151‑1 (MQ=255)
                                                                                       gtgtGAAAGAGGGGT‑TTATA‑AG‑GGCGTAATTCCCCACAGCAACCCATCCCCCTGC‑CCCCG‑GAGAAAAAACATTTA‑TTTCCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATc                                                                 <  4:1280587/145‑1 (MQ=255)
                                                                                       ggCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCAtt                                                              >  2:945237/1‑151 (MQ=255)
                                                                                         cTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTa                                                             >  4:1015029/1‑150 (MQ=255)
                                                                                          tGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTaa                                                            >  2:768196/1‑150 (MQ=255)
                                                                                            aaaGGGGGGT‑TTATT‑AG‑GGCGTAAATCCCTACAACATCACATTCCCCTGC‑CCCCC‑GACAAATAATCGTTT‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTa                                                             <  3:272398/147‑1 (MQ=255)
                                                                                            aaaGGGGGGG‑TTATT‑AA‑GCCGTAAAACCCTACAACACCCCATTCCCCTGG‑CCTCC‑GAGAAAAAATCGTTT‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTaa                                                            <  3:379741/148‑1 (MQ=255)
                                                                                               ggggggt‑TTATA‑AA‑GGCGTAAATCCCTACAGCATCCCATTCCCCTGC‑CCCCC‑GAGAAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACg                                                       <  2:1009545/148‑1 (MQ=255)
                                                                                               gagGGGT‑TTATT‑AA‑GGCGTAAATCCCAACACCATCACATTCCCCTGC‑CCCCC‑GACAAATAACCCTTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACg                                                       <  4:524104/150‑1 (MQ=255)
                                                                                                 ggggt‑TTATT‑AG‑GGCGTAAATCCCAACAGCATCACATTCCCCGGG‑CCCCC‑GAAAAAAAACCGCTA‑TTCCCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCa                                                     <  1:228356/150‑1 (MQ=255)
                                                                                                 ggggt‑TTATT‑AA‑GGCGAAAACCCCTACACCACCACTTCCCCCTGC‑CCCCC‑GAGAAATAACCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCaa                                                    <  2:1107810/151‑1 (MQ=255)
                                                                                                  gggt‑TTATA‑AG‑GGGGTAAATCCCTACAGAACCACATCCCCCGGC‑CCCCC‑GAGAAAAAAGCCTTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAAt                                                   <  2:707351/151‑1 (MQ=255)
                                                                                                  gggt‑TTATA‑AA‑GGGGTAAACCCCACCAGCATCACTTCCCCCGGG‑CCTCG‑GAGAAATAACCGTTA‑TTCCCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAAt                                                   <  4:1128560/151‑1 (MQ=255)
                                                                                                       ttatg‑aA‑GGGGTGATGCCCAACACCATCACATCCCCCTGC‑CCCCC‑GAATAATAGTCGCTC‑TTGCCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATTCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATg                                                  <  2:1448877/148‑1 (MQ=255)
                                                                                                       ttata‑aG‑GGCGTGAATCCCTACACCATCACATTCCCCTGC‑CCCCC‑GAGAAATAATCGTTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCt                                                <  3:279176/150‑1 (MQ=255)
                                                                                                         att‑aG‑GGGGTAAATCCCTACACCACCAC‑TTCCCCCGG‑CCCCCGGAAAAATAATCGTTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTgg                                              <  2:35249/150‑1 (MQ=255)
                                                                                                         att‑aA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTgg                                              >  1:17159/1‑150 (MQ=255)
                                                                                                           t‑aA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGat                                            >  2:1526238/1‑150 (MQ=255)
                                                                                                             aa‑GGCGTAAATACCTACACCATCACATCCCCCTGC‑CCCCC‑GAGAAATAATCGTTT‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTgg                                              <  3:146158/147‑1 (MQ=255)
                                                                                                              aggggtaaaaTCCCTAAAGCATCACATCCCCCGGC‑CCCCG‑GAAAAATACTCCTTT‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTgg                                              <  2:1583283/139‑1 (MQ=255)
                                                                                                                 gCGTAAATACCTACAGCACCACTTTCCCCTGG‑CCTCG‑GAGAAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCgg                                       <  3:245809/151‑1 (MQ=255)
                                                                                                                  cGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGc                                      >  3:434015/1‑151 (MQ=255)
                                                                                                                   gaaattcccTACAGCACCACTTTCCCCGGC‑CCCCC‑GAGAAATAACCGTTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGcc                                     <  4:1372509/144‑1 (MQ=255)
                                                                                                                   gTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGcc                                     >  2:441700/1‑151 (MQ=255)
                                                                                                                    tGAATCCATACAGCATCACATTCCCCGGC‑CCCCC‑GAGAAATAATCCCTA‑TACCCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGc                                      <  4:18038/149‑1 (MQ=255)
                                                                                                                    tGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCt                                    >  4:973257/1‑151 (MQ=255)
                                                                                                                      aacccccacacCACCACACTCCCCTCC‑CCCCC‑GACAAATAATCACTA‑TACCCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCg                                        <  2:1165313/138‑1 (MQ=255)
                                                                                                                      aaTCCCTACAGCACCACATCCCCCTGC‑CCCCC‑GAGAAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTgg                                  <  4:1520112/151‑1 (MQ=255)
                                                                                                                       aTACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGc                                 >  3:42074/1‑151 (MQ=255)
                                                                                                                        tACCTACAGCATCACATTACCCT‑G‑CCTCCGGAGAATTAAT‑GCTATTTACCGCATTATTCG‑GGGCGGACCGCACGTCAGCGCGGGGGGGCAGGAAACCCGCGTTAAGGCCGCCCGAATCATTAAAAAGGCAAGGCTGGAAATCGGCTTg                                   >  2:1114021/1‑148 (MQ=255)
                                                                                                                          ccaacacCATCACCTCCCCCGCG‑CCCCC‑GAGAAATAATCGTTT‑TTCCCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTAt                              <  2:1144449/148‑1 (MQ=255)
                                                                                                                              cAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCt                          >  2:87910/1‑151 (MQ=255)
                                                                                                                               aGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCtg                         >  2:172004/1‑151 (MQ=255)
                                                                                                                                ccATCACATTCCCCGGC‑CCTCC‑AAGAAATAATCGCTT‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCtgt                        <  2:1435348/150‑1 (MQ=255)
                                                                                                                                   tCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGa                      >  2:1508844/1‑150 (MQ=255)
                                                                                                                                      caTTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGaaaa                   >  3:246373/1‑150 (MQ=255)
                                                                                                                                       aTTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATc                 >  2:520144/1‑151 (MQ=255)
                                                                                                                                        tcccccTGG‑CCCCC‑GAGAAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTATGaaaa                   <  4:122143/145‑1 (MQ=255)
                                                                                                                                           ccccggcccccc‑GAGAAATAATCTTTA‑TAACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGc              <  4:1449462/141‑1 (MQ=255)
                                                                                                                                            ccTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCAt            >  4:78209/1‑151 (MQ=255)
                                                                                                                                                  ccTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGt       >  4:1126309/1‑151 (MQ=255)
                                                                                                                                                   cTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGt       >  2:1023288/1‑150 (MQ=255)
                                                                                                                                                         aGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCcac  >  1:141310/1‑150 (MQ=255)
                                                                                                                                                           |                                                                                                                                                     
ACGAAATCGACCGCCGCACGCATTTTGTTGACGCTGGATATGCAAAAGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGG‑TTATT‑AA‑GGCGTGAATACCTACAGCATCACATTACCCTGG‑CCTCC‑GAGCAATAATCGCTA‑TTACCGCCATAATCGCGGGCGCA‑CGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCAC  >  NC_000913/573259‑573556

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: ‑