Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,636,382 C→A 100% intergenic (‑15/‑374) nohQ ← / → ynfO pseudogene, Qin prophage; Phage DNA packaging protein Nu1 family;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin/uncharacterized protein, Qin prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,636,3820CA95.8% 3232.6 / 60.2 1312intergenic (‑15/‑374)nohQ/ynfOpseudogene, Qin prophage; Phage DNA packaging protein Nu1 family;Phage or Prophage Related; DNA packaging protein NU1 homolog from lambdoid prophage Qin/uncharacterized protein, Qin prophage
Reads supporting (aligned to +/- strand):  ref base C (33/30);  new base A (683/558);  total (719/593)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 1313 total aligned reads displayed.

GGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAAT  >  NC_000913/1636236‑1636530
                                                                                                                                                  |                                                                                                                                                    
 gcAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCCGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACa                                                                                                                                                    >  1:291255/1‑148 (MQ=17)
   agagTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACaaa                                                                                                                                                  >  5:307309/1‑148 (MQ=17)
   agagTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACaaa                                                                                                                                                  >  5:336056/1‑148 (MQ=17)
    gagTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACaaaaa                                                                                                                                                <  4:97303/149‑1 (MQ=17)
      gTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACTGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACa                                                                                                                                              <  5:41346/149‑1 (MQ=17)
      gTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACaaaaa                                                                                                                                                >  2:14613/1‑147 (MQ=17)
      gTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACa                                                                                                                                              <  7:235209/149‑1 (MQ=17)
      gTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACa                                                                                                                                              <  8:301661/149‑1 (MQ=17)
      gTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACa                                                                                                                                              <  4:178100/149‑1 (MQ=17)
      gTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCGTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACa                                                                                                                                              <  6:175331/149‑1 (MQ=17)
      gTCATAAAGCACATCATTACCCTTGCCACCGCCTCGAAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACa                                                                                                                                              <  3:84791/149‑1 (MQ=12)
       ttaTAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAATTGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACag                                                                                                                                             <  2:199003/147‑1 (MQ=12)
       tCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACaaaaa                                                                                                                                                >  4:320411/1‑146 (MQ=17)
       tCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACag                                                                                                                                             <  2:1096/149‑1 (MQ=17)
       tCATAAAGCACCTCATTACCCTTGCAACAGCCTCGCAGAACGGGCATTCCCTGTTCCTTCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGACAGCTGCTTTTTGTTGCCTTCAATTGTTCATTCCACGGACaaaaa                                                                                                                                                >  6:45241/1‑146 (MQ=11)
       tCATAAAGCACCTCAGTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGGTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACag                                                                                                                                             <  6:74774/149‑1 (MQ=12)
        cATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGa                                                                                                                                            >  3:326147‑M2/1‑148 (MQ=255)
         aTAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACa                                                                                                                                              >  7:115790/1‑146 (MQ=17)
         aTAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACa                                                                                                                                              >  7:218811/1‑146 (MQ=17)
         aTAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACa                                                                                                                                              >  3:14301/1‑146 (MQ=17)
         aTAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACa                                                                                                                                              >  6:123377/1‑146 (MQ=17)
         aTAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACa                                                                                                                                              >  6:147554/1‑146 (MQ=17)
         aTAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGag                                                                                                                                           >  8:240286‑M2/1‑147 (MQ=255)
         aTAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGa                                                                                                                                            >  7:197076‑M2/1‑147 (MQ=255)
         aTAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGa                                                                                                                                            >  7:349242‑M2/1‑147 (MQ=255)
          tAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGaga                                                                                                                                          >  8:49630‑M2/1‑146 (MQ=255)
          tAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCATGTTCCTGACAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGag                                                                                                                                           <  4:106976‑M2/148‑3 (MQ=255)
           aaaGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGACAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGaga                                                                                                                                          <  1:61440‑M2/148‑4 (MQ=255)
           aaaGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGCTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGaga                                                                                                                                          <  7:256700‑M2/148‑4 (MQ=255)
            aaGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaa                                                                                                                                        >  1:381403‑M2/1‑144 (MQ=255)
            aaGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaa                                                                                                                                        >  8:327161‑M2/1‑144 (MQ=255)
            aaGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaa                                                                                                                                        >  3:352713‑M2/1‑144 (MQ=255)
            aaGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaa                                                                                                                                         >  3:259175‑M2/1‑144 (MQ=255)
             aGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaag                                                                                                                                       >  6:281764‑M2/1‑143 (MQ=255)
             aGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaag                                                                                                                                       <  4:76924‑M2/149‑7 (MQ=255)
             aGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaag                                                                                                                                       >  4:26005‑M2/1‑143 (MQ=255)
             aGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaag                                                                                                                                       >  2:355676‑M2/1‑143 (MQ=255)
              gCACCTCATTACCCTTGCCCCCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaag                                                                                                                                       <  8:220791‑M2/148‑7 (MQ=255)
               cACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaag                                                                                                                                       <  6:356100‑M2/147‑7 (MQ=255)
                aCCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaagga                                                                                                                                     <  1:162336‑M2/148‑9 (MQ=255)
                aCCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaagga                                                                                                                                     <  1:242983‑M2/148‑9 (MQ=255)
                aCCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaagga                                                                                                                                     <  2:377537‑M2/148‑9 (MQ=255)
                 ccTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTCCTCCAGCCAGTTCTGAATGGTACGGAAACTCGCACCGAAAATCTCAGCCACCTGCTCTTTGTTGACTTCCATTGTTCATTCCCCGGACAAAACCACagaaaggaac                                                                                                                                   >  4:343583‑M2/1‑139 (MQ=255)
                  cTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACaaa                                                                                                                                                  <  3:10726/133‑1 (MQ=17)
                  cTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaac                                                                                                                                  >  8:72317‑M2/1‑138 (MQ=255)
                  cTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaac                                                                                                                                  >  1:267820‑M2/1‑138 (MQ=255)
                  cTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaac                                                                                                                                  <  6:89505‑M2/149‑12 (MQ=255)
                  cTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaac                                                                                                                                  <  2:179247‑M2/149‑12 (MQ=255)
                  cTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaac                                                                                                                                  <  3:193711‑M2/149‑12 (MQ=255)
                  cTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGa                                                                                                                                            <  8:135095‑M2/139‑2 (MQ=255)
                  cTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGa                                                                                                                                            >  7:135095‑M2/1‑138 (MQ=255)
                  cTCAATACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTCTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaac                                                                                                                                  <  1:188068‑M2/149‑12 (MQ=255)
                   tCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaac                                                                                                                                  <  7:385675‑M2/148‑12 (MQ=255)
                   tCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaac                                                                                                                                  <  8:390012‑M2/148‑12 (MQ=255)
                   tCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaac                                                                                                                                  <  1:100587‑M2/148‑12 (MQ=255)
                     aTTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaacgac                                                                                                                               >  4:203686‑M2/1‑135 (MQ=255)
                     aTTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaacga                                                                                                                                <  8:277196‑M2/148‑14 (MQ=255)
                      ttACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaacgaca                                                                                                                              >  3:171840‑M2/1‑134 (MQ=255)
                       tACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaacga                                                                                                                                <  3:219267‑M2/146‑14 (MQ=255)
                       tACCATTGCCAGCGCCTCGCAGAACGGGGATTCCCTGTTGCTGTCAGTTCTGATTGGTACGGATAATAGCACCGAAAATGTCAGCGAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACTTACAAAAACAGagaaaggaaacgact                                                                                                                              >  8:105194‑M2/1‑133 (MQ=255)
                       tACACTTGCCACCGCCTCGCAGAACGTGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCAAGCTTCTTTTTGTTGACTTCCATTCTTCATACCACGGACAAAAAAATagaaaggaaacga                                                                                                                                <  2:306673‑M2/146‑14 (MQ=255)
                        aCCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaacgacag                                                                                                                             >  5:159897‑M2/1‑132 (MQ=255)
                        aCCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaacgaca                                                                                                                              <  2:281102‑M2/147‑16 (MQ=255)
                        aCCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCGGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaacgaca                                                                                                                              <  5:69555‑M2/147‑16 (MQ=255)
                         cccTTGCCACCGCCTCGCAGCACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACTGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACagag                                                                                                                           <  1:312599/149‑1 (MQ=12)
                         cccTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaag                                                                                                                                       >  3:300750‑M2/1‑131 (MQ=255)
                         cccTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACagag                                                                                                                           <  8:184220/149‑1 (MQ=17)
                         cccTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACaga                                                                                                                            <  2:331598/148‑1 (MQ=17)
                          ccTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaac                                                                                                                                  <  5:11769‑M2/141‑12 (MQ=255)
                          ccTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAgg                                                                                                                          >  1:218278/1‑149 (MQ=17)
                          ccTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAgg                                                                                                                          >  5:36125/1‑149 (MQ=17)
                          ccTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAgg                                                                                                                          <  5:246175/149‑1 (MQ=17)
                          ccTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAgg                                                                                                                          >  7:73632/1‑149 (MQ=17)
                          ccTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAgg                                                                                                                          >  4:292061/1‑149 (MQ=17)
                          ccTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTGCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAACCGACCGAgg                                                                                                                          >  3:60283/1‑149 (MQ=255)
                          ccTTGCCACCGCCTCGCAGAACGGGCATACCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAgg                                                                                                                          >  5:216343/1‑149 (MQ=17)
                            ttGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTTACTTCCATTGTTCATTCAACGGACAAACACAGAGAAAGGAAACGACAGAGGc                                                                                                                         >  1:85139/1‑148 (MQ=255)
                              gCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGcca                                                                                                                       <  7:49680/148‑2 (MQ=14)
                               ccACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGccaa                                                                                                                      >  5:183828/1‑146 (MQ=11)
                                 aCCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACagag                                                                                                                           >  2:48700/1‑141 (MQ=17)
                                  ccGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaa                                                                                                                                         >  6:229446‑M2/1‑122 (MQ=255)
                                                gggCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaacgacaggtgctgaaagcgagctttttggcct                                                                                                    >  7:161058‑M2/1‑108 (MQ=255)
                                                     ttCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAAc                                                                                                                                               >  8:357445/1‑101 (MQ=17)
                                                            ttacTGCCAGTTCTGAATGTTACGGATACTCGCACCTAAAATGTCAGCCAGCTGCTTTTTGTTGACATCCATGGTTCATTACACGGCCAAAAACAGAGAAAGGAAACGACTGAGGACAAAAAGCTCGCTTTCAGCACCTGTCGTTTc                                                                                          <  2:84423/144‑1 (MQ=11)
                                                            ttCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaacgacaggtgctgaaagcgagctttttggcctctgtcgtttcc                                                                                         >  7:289213‑M2/1‑96 (MQ=255)
                                                            ttCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaacgacaggtgctgaaagcgagctttttggcctctgtcgtttc                                                                                          >  1:364746‑M2/1‑96 (MQ=255)
                                                            ttCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaacgacaggtgctgaaagcgagctttttggcctctgtcgtttc                                                                                          >  4:59278‑M2/1‑96 (MQ=255)
                                                            ttCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaacgacaggtgctgaaagcgagctttttagcctctgtcgtttc                                                                                          >  4:231965‑M2/1‑96 (MQ=255)
                                                               cTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaagg                                                                                                                                      <  7:326146‑M2/100‑8 (MQ=255)
                                                                tGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCATTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGcca                                                                                                                       >  8:160435/1‑113 (MQ=9)
                                                                 gCCAGTGATGAATGGTACGGATACTCGCACCGTAAATGTCAGCCATCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACagag                                                                                                                           <  6:321278/109‑1 (MQ=11)
                                                                      ttCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaacgaca                                                                                                                              <  1:339814‑M2/101‑16 (MQ=255)
                                                                               ggaCGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAACACAGAGAAAGGAAACGACAGAGGCCAAAAAGCACGCTTTCAGCACCTGTCGGTTCCTTTCTATTCAGAGGGTAatt                                                                     >  6:25510/3‑149 (MQ=11)
                                                                               gTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGTGGGTAttt                                                                     >  6:204474/1‑149 (MQ=255)
                                                                               gTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGAGAAAGGAAACGACAGAGGc                                                                                                                         >  5:143236/1‑97 (MQ=18)
                                                                                 aCGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaa                                                                                                                                    <  1:128873‑M2/84‑10 (MQ=255)
                                                                                         tCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaa                                                                                                                                    >  2:228343‑M2/1‑67 (MQ=255)
                                                                                                 aaaaTGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACAGagaaaggaaacgacaggtgctgaaagcgagctttttggcctctgtcgtttc                                                                                          >  4:148871‑M2/1‑59 (MQ=255)
                                                                                                       tCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACAAAAACag                                                                                                                                             >  2:248320/1‑53 (MQ=255)
                                                                                                         agccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGACaa                                                                                                                                                   >  5:273092/1‑45 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGCAGAGTCATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAAT  >  NC_000913/1636236‑1636530

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