Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 289,256 G→C 8.2% intergenic (‑94/+45) yagI ← / ← argF CP4‑6 prophage; putative DNA‑binding transcriptional regulator/ornithine carbamoyltransferase 2, chain F; CP4‑6 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913289,2560GC 8.2% 204.6 / 10.8 73intergenic (‑94/+45)yagI/argFCP4‑6 prophage; putative DNA‑binding transcriptional regulator/ornithine carbamoyltransferase 2, chain F; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base G (52/14);  new base C (3/3);  total (55/18)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

TCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTC  >  NC_000913/289107‑289403
                                                                                                                                                      |                                                                                                                                                   
tCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTTAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCtg                                                                                                                                                     >  2:1203714/1‑150 (MQ=255)
  gTTGAACAGGTCGAGGATCTGCAACGCACGTTCAGCAGACTGAATAATCGGCATAATGCTGGCATGTACACGCTGGAAATAAGGCGAAAACGTACATTATTCGGGATGAAAAATCAATGAAATGGAGCAGGTGTTCTCCCTCTCCCtgtg                                                                                                                                                   >  2:1806422/1‑150 (MQ=255)
  gTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGTCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCtgtg                                                                                                                                                   <  1:34757/150‑1 (MQ=255)
       aCATGTCGAGGAGCTGCAACGCACGTTCAACAGACTGAATAATAGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgaga                                                                                                                                              <  2:899956/150‑1 (MQ=255)
        cAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGagag                                                                                                                                             >  2:1706654/1‑150 (MQ=255)
         tgttCGAGTATTTGAAACCAAAGTTCAAAATATTTAATATTGGGAATATTGCTGCCATTACCAGCATGTAATTAAGGCTAAAAGGTACCTTTTTCGGGATTAAAATTAAATCAATGGAGCCGGTGTTCTCCCTCTCCCACAGGGAGAggg                                                                                                                                           <  2:692661/147‑1 (MQ=255)
         aGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCACAGGGAGAgg                                                                                                                                            <  1:1664049/150‑1 (MQ=255)
          gggcGAGGAGCGGCCACGGACGGTCAACAGACTGAATAAGCGGCATAATGCCGGCTTGTCCACGCTGGAATTATGGCGAAAACGTACCTTTTTCGGGATGCAAAATCAATGAACTGGAGCCGGTGTTCTCCCTCTCCCACAGGGAGAggg                                                                                                                                           <  2:145164/147‑1 (MQ=255)
          ggTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGTATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAggg                                                                                                                                           >  2:896603/1‑150 (MQ=255)
           gTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAggga                                                                                                                                          >  1:1822780/1‑149 (MQ=255)
             cGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                        >  2:1590873/1‑150 (MQ=255)
             cGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                        >  2:661365/1‑150 (MQ=255)
             cGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                        >  1:222438/1‑150 (MQ=255)
               aGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCTTGTCCCCGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAAGGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                      >  1:1222226/1‑150 (MQ=255)
               aGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                      >  2:236518/1‑150 (MQ=255)
                ggATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTGCCTTTTTAGGGATGGTAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCATGTGGGAGAGGGCCgggg                                                                                                                                     >  2:267625/1‑150 (MQ=255)
                ggATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCgggg                                                                                                                                     >  1:203661/1‑150 (MQ=255)
                   tCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGa                                                                                                                                  >  1:264005/1‑150 (MQ=255)
                   tCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGa                                                                                                                                  >  2:134831/1‑150 (MQ=255)
                     tGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAACTCAATGAAATGGAGCCGCTGTTCTCCCTCTCCCTGAGGGAGAGGGCCGGGGTGAgg                                                                                                                                >  2:484370/1‑150 (MQ=255)
                     tGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                >  2:1835927/1‑150 (MQ=255)
                     tGCAACCCACGTTCAACAGACGGAATAATCGGCATACTGCTGGCATGTCCACGCTGTAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAAGGGAGCCGGTGTTCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                <  1:1593188/150‑1 (MQ=255)
                     ctgaaCGCACGTTCAACATACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                <  2:695384/147‑1 (MQ=255)
                       cAACGCTCGTTCAACAGACAGAATAAACGGCAATATCCTGGCAGGTCCACGCGGGAATTAAGGCGAAAACGTTCCTTTTTCGGGATGAAAAATCAATGTAATGGAGCCGTTGTTCTCCCTCTCCCTGTGGGAGAGGGCCTGGGTGAGGGc                                                                                                                              <  2:220061/150‑1 (MQ=255)
                       cAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGc                                                                                                                              <  1:989469/150‑1 (MQ=255)
                        aaCGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                             >  2:1317741/1‑150 (MQ=255)
                         ccGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATCAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTTTGGGAGAGGGCCGGGGTGAGGGCAc                                                                                                                            <  2:1209661/149‑1 (MQ=255)
                          cGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTAGCTTTTTCGGGATGAAAAATCAATTAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAcc                                                                                                                           <  2:1085078/150‑1 (MQ=255)
                                 tCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGTTGAGGGCACCAACgcgc                                                                                                                    >  2:923132/1‑150 (MQ=255)
                                  cAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCACCgcgcc                                                                                                                   >  2:192586/1‑149 (MQ=255)
                                   aaCAGACTGAATAATCGGCATAATGCTGGCATGCCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAATAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACgcgctc                                                                                                                  >  2:19875/1‑148 (MQ=255)
                                          tGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGAccc                                                                                                           >  2:1213674/1‑150 (MQ=255)
                                           gAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCa                                                                                                          >  2:1087377/1‑150 (MQ=255)
                                                atCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCAc                                                                                                     >  2:1426993/1‑150 (MQ=255)
                                                   ggCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTcc                                                                                                  >  2:1185716/1‑150 (MQ=255)
                                                    ccATAAGGCTGCCATGCCAACCCGGCAATTACGGCAAAACCGACCCTTTTCGGGTATAAAAAATCATTCAAATGGCGCCGGTGTGCTCCCCCTCCCTGGGGGTGAGGGCCGGGGGGAGGGCACCAACGCGCAGCGGACCCAATCACTccc                                                                                                 <  2:961134/149‑1 (MQ=255)
                                                          tGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAAAGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCg                                                                                           >  1:1814545/1‑150 (MQ=255)
                                                             tGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTg                                                                                        >  1:59420/1‑150 (MQ=255)
                                                              ggCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGATGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGc                                                                                       >  1:224329/1‑150 (MQ=255)
                                                                cATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCAGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCca                                                                                     >  2:392576/1‑150 (MQ=255)
                                                                     ccACGCTGGAATTAAGGCGAAACCGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCCACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCatcatc                                                                                >  2:635269/1‑150 (MQ=255)
                                                                     ccACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCATCGGACCCAATCACACCCCAAGCGTTGCcatcata                                                                                >  1:1120687/1‑149 (MQ=255)
                                                                      cACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCAtcatca                                                                               >  2:1788625/1‑150 (MQ=255)
                                                                          cTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACGCCCCAAGCGTTGCCATCATCACCGc                                                                           >  2:377745/1‑150 (MQ=255)
                                                                              aaTTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGTGCACCAACGCGCAGCGGACCCAATCACTCCCCACGCGTTGCCATCATCACCGCCTTa                                                                       >  1:175498/1‑150 (MQ=255)
                                                                                   aGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACACCCCAAGCGTTGCCATCATCACCGCCTTAATCGt                                                                  >  2:1660716/1‑150 (MQ=255)
                                                                                    ggCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTa                                                                 >  1:1117499/1‑150 (MQ=255)
                                                                                              tACCTTTTCGGGTATTAAAAATCAATTAATTCGCGCCGGTGTCCCCCCTCTCTCTGGGTGTGTGGGCGGGGTTGCGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGt                                                       <  2:110143/150‑1 (MQ=255)
                                                                                                ccTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGCCCCAATCACTCCCCCAGCGTTGCCTTCATCACCGCCTTACTCGTATGCCTCCGGttt                                                     >  1:891586/1‑150 (MQ=38)
                                                                                                ccTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGttt                                                     >  1:1690564/1‑150 (MQ=255)
                                                                                                        gggATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTg                                             >  1:786681/1‑150 (MQ=255)
                                                                                                         ggATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGATTTCCGCCTgg                                            >  2:174746/1‑150 (MQ=255)
                                                                                                            aTGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTc                                          >  2:1610736/1‑150 (MQ=255)
                                                                                                            aTGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTc                                          >  1:1052731/1‑150 (MQ=255)
                                                                                                             tGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCg                                         >  1:582522/1‑150 (MQ=255)
                                                                                                             tGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCg                                         >  2:1100483/1‑150 (MQ=255)
                                                                                                              caaaaaTCAATGAAAGGGCGCCGGTGTCCCCCCTCTCCCTGTGGGCGTGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGa                                        <  1:5108/149‑1 (MQ=255)
                                                                                                                aaaaTCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCGTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAac                                      >  2:837945/1‑150 (MQ=255)
                                                                                                                aaaaTCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGACGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAac                                      <  2:1822780/150‑1 (MQ=255)
                                                                                                                aaaaTAAATTAAATGGAGCCGGTGTTCTCCCGCTCCCGGTGGGAGTGGGCGGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAac                                      <  2:1058060/150‑1 (MQ=255)
                                                                                                                     cAAAGAAATCGAGCCGGTGATCTCCATCTCCCGGTGGGAGAGGGCCCGGGTGAGGGCACCAACGCGCAGCGGCCCCAATCACTCCCCAAGCGTTGCCATCTTCACCGCCTTAATCGTTTGCTTCCGGGTGTCGCCCTGGTCGAACACgcg                                 >  2:1085184/1‑148 (MQ=255)
                                                                                                                         gAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCCCCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTc                             >  1:1365928/1‑150 (MQ=255)
                                                                                                                           aaTGGAGCAGGTGTTCTCACTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCACCACCGCGCAGCGGACCCCATCACTCCCCAAGCGTTGCCATCATCCCCGCCTTAAGCGTATGCAGCCGGTTTTCCGGCTGGTCGAAAGAGATGCTcgc                           >  2:992734/1‑150 (MQ=37)
                                                                                                                            agggtgCCGGGGTTCTCCCTCTCCCGGTGGCTGCGGGCCGGGGTGCGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACTGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTcgcc                          <  2:264005/145‑1 (MQ=255)
                                                                                                                             tGGCGCCGGGGCTCACCCTCCCCCTGTGGGAGTGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGTTcgccg                         <  2:20957/150‑1 (MQ=255)
                                                                                                                                       gtTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCCAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCaaa               >  1:1060123/1‑150 (MQ=255)
                                                                                                                                        ccggccctctcccaCAGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGCCCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATACGGTTTTCCGCCTGGTAGAACACGATGCTCACCGCCGACTCAAac              >  2:1503211/5‑150 (MQ=255)
                                                                                                                                         tctccctctcccTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAaca             >  2:707356/1‑150 (MQ=255)
                                                                                                                                          ggccctctcccaCAGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCCCTCCCCAAGCGTTGGCATCATCGCCGCCTTAATCGTATGCAGGCGGGTTTCCGCGTGGTCGATCACGATGCTCGCCGGCGAGTCAAaccc            >  1:242965/3‑148 (MQ=255)
                                                                                                                                          ggccctctcccaCAGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGCCCACGCTGCTCGCCGCCGCCTCAAacac            >  1:1823869/3‑150 (MQ=255)
                                                                                                                                          ctccctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTCTGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGACGACTCaacccc            >  1:813914/1‑146 (MQ=255)
                                                                                                                                          ctccctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAacac            >  2:1764818/1‑150 (MQ=255)
                                                                                                                                             cctctccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAAGCACCCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAAGACCTc         >  2:1207784/1‑150 (MQ=255)
                                                                                                                                                 tccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTcgtcc     >  2:245422/1‑150 (MQ=255)
                                                                                                                                                  ccctGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTAGGCATCCGGTGTTCCGCCTGGTCGACCACGATGCTCGCCGCCGACTCAAACCCATcgtcgg    >  1:916325/1‑148 (MQ=255)
                                                                                                                                                    ctGTGGGAGAGGGCGGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGTTCGCCGCCGACTCAAACACCTcgtccgtc  <  2:1066538/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
TCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGGATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTC  >  NC_000913/289107‑289403

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