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

TCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTC  >  NC_000913/289107‑289403
                                                                                                                                                      |                                                                                                                                                   
tCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTTAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCtg                                                                                                                                                     >  2:1203714/1‑150 (MQ=255)
  gTTGAACAGGTCGAGGATCTGCAACGCACGTTCAGCAGACTGAATAATCGGCATAATGCTGGCATGTACACGCTGGAAATAAGGCGAAAACGTACATTATTCGGG‑ATGAAAAATCAATGAAATGGAGCAGGTGTTCTCCCTCTCCCtgtg                                                                                                                                                   >  2:1806422/1‑150 (MQ=255)
  gTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGTCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCtgtg                                                                                                                                                   <  1:34757/150‑1 (MQ=255)
       aCATGTCGAGGAGCTGCAACGCACGTTCAACAGACTGAATAATAGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGgaga                                                                                                                                              <  2:899956/150‑1 (MQ=255)
        cAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGagag                                                                                                                                             >  2:1706654/1‑150 (MQ=255)
         tgttCGAGTATTTGAAACCAAAGTTCAAAATATTTAATATTGGGAATATTGCTGCCATTACCAGCATGTAATTAAGGCTAAAAGGTACCTTTTTCGGG‑ATTAAAATTAAAT‑CAATGGAGCCGGTGTTCTCCCTCTCCCACAGGGAGAggg                                                                                                                                           <  2:692661/147‑1 (MQ=255)
         aGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCACAGGGAGAgg                                                                                                                                            <  1:1664049/150‑1 (MQ=255)
          gggcGAGGAGCGGCCACGGACGGTCAACAGACTGAATAAGCGGCATAATGCCGGCTTGTCCACGCTGGAATTATGGCGAAAACGTACCTTTTTCGGG‑ATGCAAAATCAATGAACTGGAGCCGGTGTTCTCCCTCTCCCACAGGGAGAggg                                                                                                                                           <  2:145164/147‑1 (MQ=255)
          ggTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGT‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAggg                                                                                                                                           >  2:896603/1‑150 (MQ=255)
           gTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAggga                                                                                                                                          >  1:1822780/1‑149 (MQ=255)
             cGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                        >  2:1590873/1‑150 (MQ=255)
             cGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                        >  2:661365/1‑150 (MQ=255)
             cGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCg                                                                                                                                        >  1:222438/1‑150 (MQ=255)
               aGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCTTGTCCCCGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAAGGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                      >  1:1222226/1‑150 (MQ=255)
               aGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                      >  2:236518/1‑150 (MQ=255)
                ggATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTGCCTTTTTAGGG‑ATGGTAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCATGTGGGAGAGGGCCgggg                                                                                                                                     >  2:267625/1‑150 (MQ=255)
                ggATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCgggg                                                                                                                                     >  1:203661/1‑150 (MQ=255)
                   tCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGa                                                                                                                                  >  1:264005/1‑150 (MQ=255)
                   tCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGa                                                                                                                                  >  2:134831/1‑150 (MQ=255)
                     tGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAACTCAATGAAATGGAGCCGCTGTTCTCCCTCTCCCTGAGGGAGAGGGCCGGGGTGAgg                                                                                                                                >  2:484370/1‑150 (MQ=255)
                     tGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                >  2:1835927/1‑150 (MQ=255)
                     tGCAACCCACGTTCAACAGACGGAATAATCGGCATACTGCTGGCATGTCCACGCTGTAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAAGGGAGCCGGTGTTCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                <  1:1593188/150‑1 (MQ=255)
                     ctgaaCGCACGTTCAACATACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAgg                                                                                                                                <  2:695384/147‑1 (MQ=255)
                       cAACGCTCGTTCAACAGACAGAATAAACGGCAATATCCTGGCAGGTCCACGCGGGAATTAAGGCGAAAACGTTCCTTTTTCGGG‑ATGAAAAATCAATGTAATGGAGCCGTTGTTCTCCCTCTCCCTGTGGGAGAGGGCCTGGGTGAGGGc                                                                                                                              <  2:220061/150‑1 (MQ=255)
                       cAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGc                                                                                                                              <  1:989469/150‑1 (MQ=255)
                        aaCGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                             >  2:1317741/1‑150 (MQ=255)
                         ccGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATCAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTTTGGGAGAGGGCCGGGGTGAGGGCAc                                                                                                                            <  2:1209661/149‑1 (MQ=255)
                          cGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTAGCTTTTTCGGG‑ATGAAAAATCAATTAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAcc                                                                                                                           <  2:1085078/150‑1 (MQ=255)
                                 tCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGTTGAGGGCACCAACgcgc                                                                                                                    >  2:923132/1‑150 (MQ=255)
                                  cAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCACCgcgcc                                                                                                                   >  2:192586/1‑149 (MQ=255)
                                   aaCAGACTGAATAATCGGCATAATGCTGGCATGCCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAATAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACgcgctc                                                                                                                  >  2:19875/1‑148 (MQ=255)
                                          tGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGAccc                                                                                                           >  2:1213674/1‑150 (MQ=255)
                                           gAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCa                                                                                                          >  2:1087377/1‑150 (MQ=255)
                                                atCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCAc                                                                                                     >  2:1426993/1‑150 (MQ=255)
                                                   ggCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTcc                                                                                                  >  2:1185716/1‑150 (MQ=255)
                                                    ccATAAGGCTGCCATGCCAACCCGGCAATTACGGC‑AAAACCGACCCTTTTCGGGTATAAAAAATCATTCAAATGGCGCCGGTGTGCTCCCCCTCCCTGGGGGTGAGGGCCGGGGGGAGGGCACCAACGCGCAGCGGACCCAATCACTccc                                                                                                 <  2:961134/149‑1 (MQ=255)
                                                          tGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAAAGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCg                                                                                           >  1:1814545/1‑150 (MQ=255)
                                                             tGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTg                                                                                        >  1:59420/1‑150 (MQ=255)
                                                              ggCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGATGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGc                                                                                       >  1:224329/1‑150 (MQ=255)
                                                                cATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCAGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCca                                                                                     >  2:392576/1‑150 (MQ=255)
                                                                     ccACGCTGGAATTAAGGCGAAACCGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCCACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCatcatc                                                                                >  2:635269/1‑150 (MQ=255)
                                                                     ccACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCATCGGACCCAATCACACCCCAAGCGTTGCcatcata                                                                                >  1:1120687/1‑149 (MQ=255)
                                                                      cACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCAtcatca                                                                               >  2:1788625/1‑150 (MQ=255)
                                                                          cTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACGCCCCAAGCGTTGCCATCATCACCGc                                                                           >  2:377745/1‑150 (MQ=255)
                                                                              aaTTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGTGCACCAACGCGCAGCGGACCCAATCACTCCCCACGCGTTGCCATCATCACCGCCTTa                                                                       >  1:175498/1‑150 (MQ=255)
                                                                                   aGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACACCCCAAGCGTTGCCATCATCACCGCCTTAATCGt                                                                  >  2:1660716/1‑150 (MQ=255)
                                                                                    ggCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTa                                                                 >  1:1117499/1‑150 (MQ=255)
                                                                                              tACCTTTTCGGGT‑ATTAAAAATCAATTAATTCGCGCCGGTGTCCCCCCTCTCTCTGGGTGTGTGGGCGGGGTTGCGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGt                                                       <  2:110143/150‑1 (MQ=255)
                                                                                                ccTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGCCCCAATCACTCCCCCAGCGTTGCCTTCATCACCGCCTTACTCGTATGCCTCCGGttt                                                     >  1:891586/1‑150 (MQ=38)
                                                                                                ccTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGttt                                                     >  1:1690564/1‑150 (MQ=255)
                                                                                                        ggg‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTg                                             >  1:786681/1‑150 (MQ=255)
                                                                                                         gg‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGATTTCCGCCTgg                                            >  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)
                                                                                                                                                      |                                                                                                                                                   
TCGTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTC  >  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: ‑