Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,647,187 T→C 5.7% pseudogene (170/183 nt) intK ← pseudogene, integrase fragment, Qin prophage;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,647,1870TC 5.7% 225.0 / 2.2 72pseudogene (170/183 nt)intKpseudogene, integrase fragment, Qin prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (37/30);  new base C (2/2);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

TCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGC  >  NC_000913/1646908‑1647455
                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                            
tCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATc                                                                                                                                                                                                                                                                             <  2:1175084/281‑1 (MQ=255)
   aTTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGTTATGCACTCAACCGAAACGATCTTGGTTTAAATCTTTTTTTTCggg                                                                                                                                                                                                                                                                          >  2:1057092/1‑281 (MQ=255)
   aTTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCggg                                                                                                                                                                                                                                                                          >  1:136791/1‑281 (MQ=255)
    ttATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCggg                                                                                                                                                                                                                                                                          >  2:319215/1‑280 (MQ=255)
            tGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGctt                                                                                                                                                                                                                                                                 >  2:679274/1‑281 (MQ=255)
              cAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGcttct                                                                                                                                                                                                                                                               >  1:324089/1‑281 (MQ=255)
                            gTCACTTTTGTTGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGGGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCAATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTc                                                                                                                                                                                                                                                 <  1:1185542/281‑1 (MQ=255)
                            gTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGGACCGCCCAGACTGACTCAGTGCGTTCTTGGCATATCTCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACAAGAAACGCGAATTAAGAACCTAACCTATGAATGGCTACACACTCAACAGAATAGATCTTGGTTTCAATCTTTTTCAACGGGATCAGGTTTATTTTTAGGTAACTTc                                                                                                                                                                                                                                                 >  2:266452/1‑281 (MQ=255)
                            gTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTc                                                                                                                                                                                                                                                 >  2:1209229/1‑281 (MQ=255)
                            gTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTc                                                                                                                                                                                                                                                 >  1:660387/1‑281 (MQ=255)
                            gTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCCCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATCTCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATAGTCACTCAACCGTATCTATCTTGGTTTCAATCTTTTTTATCGGGATCAGGATTTTTTTTAGGTAAATTc                                                                                                                                                                                                                                                 >  2:1185542/1‑281 (MQ=255)
                            gTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCCCTTAGCGGACATTAGCATAGGCTACTTACCATAACGCCTCATTACGCGCCCCCCCAAGAATGAATCAGCGCGTTTACGGCATATCCCCGGTAAAACCAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGAATTAAGAATCAACCATATGAAGGGTTATGCACTCAACCGAATCGATCGTGGCTTCAATTTTTCTTATCGGGACCAGGCTTCTTGCTAGGTAACTTc                                                                                                                                                                                                                                                 >  2:448945/1‑281 (MQ=255)
                                  ttttTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGc                                                                                                                                                                                                                                           <  2:136791/281‑1 (MQ=255)
                                           aCAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTCAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGCATCGATCTTGGTTTCAATCTTTTTTATCGGGATCa                                                                                                                                                                                                                                                                      >  1:305118/1‑245 (MQ=255)
                                               tCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTACCGGGATCAGGCTTCTTTTTAGGTAACTTCg                                                                                                                                                                                                                                                >  1:316892/1‑263 (MQ=255)
                                               tCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTACCGGGATCAGGCTTCTTTTTAGGTAACTTCg                                                                                                                                                                                                                                                <  2:316892/263‑1 (MQ=255)
                                                      aaCTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTTAATTTTTTTTATCGGGATCAGGCTCCTTTTTAGGTAACTTGGGGGGATTAACTTGCTAATGACTTTg                                                                                                                                                                                                                       >  2:923027/1‑281 (MQ=255)
                                                         tGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGt                                                                                                                                                                                                                    >  1:244911/1‑281 (MQ=255)
                                                          gCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGtt                                                                                                                                                                                                                   <  1:1209229/281‑1 (MQ=255)
                                                           ccAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTACCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTc                                                                                                                                                                                                                  >  1:495924/1‑281 (MQ=255)
                                                            cAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCg                                                                                                                                                                                                                 <  2:1206493/281‑1 (MQ=255)
                                                                 tCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGcgcg                                                                                                                                                                                                            <  1:131034/281‑1 (MQ=255)
                                                                 tCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGcgcg                                                                                                                                                                                                            <  2:418943/281‑1 (MQ=255)
                                                                 tCCGCTTCATGCTCTTATCAGACGTTATACTGGGCTATTTACCATAACTCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCATTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGcgcg                                                                                                                                                                                                            <  1:266452/281‑1 (MQ=255)
                                                                   cGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGACACGCGACTTAACAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGTTGACTACTTTGCGTTCGGGGCGTa                                                                                                                                                                                                          >  2:404235/1‑281 (MQ=255)
                                                                                tAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGt                                                                                                                                                                                             <  1:653805/281‑1 (MQ=255)
                                                                                       cATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCttt                                                                                                                                                                                      <  1:1139754/281‑1 (MQ=255)
                                                                                            gCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGtt                                                                                                                                                                                 <  1:1057092/281‑1 (MQ=255)
                                                                                             cATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTACCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGtt                                                                                                                                                                                 <  2:495924/280‑1 (MQ=255)
                                                                                                aGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAAc                                                                                                                                                                             <  2:324089/281‑1 (MQ=255)
                                                                                                           cacatagCCTCATTGGGCGCTCCGACCGGACTGACCCAGTGCGTTTCTAGCATATCCCCGGTAAAACAAGGAACAAACCGCCCGAGAATCAACGCCTTATACCCGACTTCAGAATCTTAACTATGGTTGGTTATGCACTCAACCGAATCGATCTTGGTTTCAATCTGTTTTATCGGGATCAGGCTTCTTTTTAGGTAAGGTGGGGGGCTTTACTTGCTGATGACTTTGCGTTCGGGGCGTAAGCCAGGGATGGACAGCGTTAGGTCTAACATAGGAtgtt                                                                                                                                                                   <  2:313616/274‑1 (MQ=255)
                                                                                                                     cATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGtt                                                                                                                                                                                 >  2:971807/1‑256 (MQ=255)
                                                                                                                     cATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGtt                                                                                                                                                                                 <  1:971807/256‑1 (MQ=255)
                                                                                                                     cATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTAGGCGCGTAAGCCAGGGTTGGTCAGCTTTCGGTTTAACATAGTCTTTTGAGCGTAAATc                                                                                                                                                        >  1:1056388/1‑281 (MQ=255)
                                                                                                                                        actgactCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACt                                                                                                                                     <  1:319215/281‑1 (MQ=255)
                                                                                                                                             ctCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAAc                                                                                                                                                                             >  2:56286/1‑236 (MQ=255)
                                                                                                                                             ctCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAAc                                                                                                                                                                             <  1:56286/236‑1 (MQ=255)
                                                                                                                                             ctCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGtt                                                                                                                                 >  2:335557/1‑280 (MQ=255)
                                                                                                                                                  cgcgTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAggg                                                                                                                                                                                                  <  2:1138835/210‑1 (MQ=255)
                                                                                                                                                  cgcgTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAggg                                                                                                                                                                                                  >  1:1138835/1‑210 (MQ=255)
                                                                                                                                                  cgcgTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTAGTTGGAGCGTAAATCAATACGGGCATTATCCACTAGTTCATGGa                                                                                                                           >  2:1035288/1‑281 (MQ=255)
                                                                                                                                                    cgTTTCTGGCATATCCCCGGTAAAGCAAGTAACAAACCACCCGAAGATGGACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTAACTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGac                                                                                                                          <  1:303677/280‑1 (MQ=255)
                                                                                                                                                      tttCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCtga                                                                                                                                                                                                                               >  1:169773/1‑177 (MQ=255)
                                                                                                                                                      tttCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCtga                                                                                                                                                                                                                               <  2:169773/177‑1 (MQ=255)
                                                                                                                                                               tatCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGAGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTCGCTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTGAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCGTGAGGATTAATATAGACTTTTGAGCGTAAACCAAAAAGGGCATTATCCACTCGTTCATGGATACTCTTTtaatc                                                                                                              >  2:777105/1‑281 (MQ=255)
                                                                                                                                                                 tCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGt                                                                                                                                                                                             <  2:726167/200‑1 (MQ=255)
                                                                                                                                                                 tCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGt                                                                                                                                                                                             >  1:726167/1‑200 (MQ=255)
                                                                                                                                                                      ggTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCCTCATCCAg                                                                                                        >  1:309509/1‑280 (MQ=255)
                                                                                                                                                                      ggTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGATCCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGATTTTGCGTTCGGCGCGTAAGCAAGGGCTGGTCAGTCATAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGAGCATTATCCACTAGTTCATGGACAATATTTTCATATTCCAg                                                                                                        >  2:364047/1‑280 (MQ=255)
                                                                                                                                                                        aaaaaCAATAACCAACCACCCGAAAATTCACACCAGAACCGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGGTCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGCCTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGt                                                                                                    <  2:244911/280‑1 (MQ=255)
                                                                                                                                                                            aCAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGc                                                                                                 <  1:923027/281‑1 (MQ=255)
                                                                                                                                                                                         caccCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCa                                                                                    <  1:679274/281‑1 (MQ=255)
                                                                                                                                                                                               aaaaTGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTCCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACt                                                                              >  2:651456/1‑281 (MQ=255)
                                                                                                                                                                                                         ccAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTGCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATTGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGa                                                                    <  2:660387/281‑1 (MQ=255)
                                                                                                                                                                                                                 gcgACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCa                                                            <  1:335557/281‑1 (MQ=255)
                                                                                                                                                                                                                    aCTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCa                                                            >  2:706123/1‑278 (MQ=255)
                                                                                                                                                                                                                      ttAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGTTCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTTATGACTTTGCGTTCAGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACg                                                          <  1:448945/278‑1 (MQ=255)
                                                                                                                                                                                                                               tACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCGTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGTCTCCATAATTGCAGGCACTAGCGCCGTGAAATTTTCACGCTTACCCCTgg                                               >  1:296736/1‑280 (MQ=255)
                                                                                                                                                                                                                                  ccTATGAATGGATATGCACTCAACCGAATCGATCTTGGATTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACCTTCCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATGATCTACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTCCAGGCACTCGCGCGGTGAACTTTGCACGCTTACCCCTGGTGTc                                           >  1:80413/1‑281 (MQ=255)
                                                                                                                                                                                                                                      tGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCAc                                                                               <  2:935411/241‑1 (MQ=255)
                                                                                                                                                                                                                                      tGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCAc                                                                               >  1:935411/1‑241 (MQ=255)
                                                                                                                                                                                                                                         aTGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTAAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCAATAGAGCCGTGAACTTTTCACGCTTATCCCTGGTTTCGAAAGcc                                    >  1:811088/1‑281 (MQ=255)
                                                                                                                                                                                                                                           ggATATGCACTCAACCGAATCGATCTTGGTTCCAATCTTTTTTATCGTGATCAGGCTTCTTTTTCGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATTAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCt                                   >  1:900323/1‑280 (MQ=255)
                                                                                                                                                                                                                                            gATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGATTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACAATCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGGTTCGATAGCTTTc                                 >  2:322621/1‑281 (MQ=255)
                                                                                                                                                                                                                                                tGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTCCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCTCTCGTTCATGGACACTCTTTTCATCATCCAGTGGTACCCTACATAATTCCAGTCACTAGCGCCGTGAACTCTTCACGCTCATCCCTGGTTTCGATGGCCTTCCAGc                             >  2:1211867/1‑281 (MQ=255)
                                                                                                                                                                                                                                                       aaCCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGTGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCAACCGTTCaa                       >  1:552471/1‑280 (MQ=255)
                                                                                                                                                                                                                                                       aaCCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGcgcg                                                                                                                                                                                                            <  2:357539/99‑1 (MQ=255)
                                                                                                                                                                                                                                                       aaCCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGcgcg                                                                                                                                                                                                            >  1:357539/1‑99 (MQ=255)
                                                                                                                                                                                                                                                       aaCCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCaa                       >  1:795333/1‑280 (MQ=255)
                                                                                                                                                                                                                                                                           cAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAAGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGc  <  1:651456/281‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aaTCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCt                                                                                                                                                                                        >  2:682308/1‑98 (MQ=255)
                                                                                                                                                                                                                                                                            aaTCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCt                                                                                                                                                                                        <  1:682308/98‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                            
TCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATCCCCGGTAAAACAAGTAACAAACCACCCGAAAATGAACACCAGAAACGCGACTTAAGAATCTACCCTATGAATGGATATGCACTCAACCGAATCGATCTTGGTTTCAATCTTTTTTATCGGGATCAGGCTTCTTTTTAGGTAACTTCGGGGGCTTAACTTGCTGATGACTTTGCGTTCGGCGCGTAAGCCAGGGATGGTCAGCTTTAGGTTTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGC  >  NC_000913/1646908‑1647455

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: