Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 289,257 T→A 9.2% intergenic (‑95/+44) 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,2570TA 9.2% 178.6 / 7.1 73intergenic (‑95/+44)yagI/argFCP4‑6 prophage; putative DNA‑binding transcriptional regulator/ornithine carbamoyltransferase 2, chain F; CP4‑6 prophage
Reads supporting (aligned to +/- strand):  ref base T (51/13);  new base A (4/3);  total (55/18)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTC  >  NC_000913/289109‑289403
                                                                                                                                                     |                                                                                                                                                  
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                                                                                                                                  >  2:134831/1‑150 (MQ=255)
                 tCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGa                                                                                                                                  >  1:264005/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)
                                                                                                                                                     |                                                                                                                                                  
GTTGAACAGGTCGAGGATCTGCAACGCACGTTCAACAGACTGAATAATCGGCATAATGCTGGCATGTCCACGCTGGAATTAAGGCGAAAACGTACCTTTTTCGGG‑ATGAAAAATCAATGAAATGGAGCCGGTGTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCACCAACGCGCAGCGGACCCAATCACTCCCCAAGCGTTGCCATCATCACCGCCTTAATCGTATGCATCCGGTTTTCCGCCTGGTCGAACACGATGCTCGCCGCCGACTCAAACACCTCGTCCGTC  >  NC_000913/289109‑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: ‑