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

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