Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,596,501 T→A 6.1% pseudogene (1487/3861 nt) yneO ← pseudogene, AidA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,596,5010TA 6.1% 199.1 / 2.0 64pseudogene (1487/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base A (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

CCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAACGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCGTGGTGATGCTGACATTCCCGACGCCGG  >  NC_000913/1596204‑1596801
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                            
ccACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAACGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACCGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTtc                                                                                                                                                                                                                                                                                                           >  1:320225/1‑301 (MQ=255)
       ggtcATTCAGCGTACCGTTATTCGTCAGGGTTAACGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGAt                                                                                                                                                                                                                                                                                                    <  1:204619/298‑1 (MQ=255)
        gACATTCAGCGTACCGTTATTCATCAGGGTTAACGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCAGTTGAAATATTAACTATCCCATTACCGAAGGCATTtcacccaaaa                                                                                                                                                                                                                                                                                                   >  1:304375/1‑293 (MQ=255)
          cATTCAGCGTACCGTTATTCGTCAGGGTTAACGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTATCGGAGGAAATGCCTCCGGTaag                                                                                                                                                                                                                                                                                                 <  2:88385/301‑2 (MQ=255)
             tctgCGTCCCGTTATCTGTCAGGTTTATCGTACCCGTACCTGATGGTCCTACTGTCATATTGAATGGTTCACGAAGGGAGTTCTGCCCAACCACCATCACGTCGCTCTTACTGTTGTCATTAAGAGCTATCTGGGTATCACGCGCTTTAACGATACCTCCGGTGGTAATATTCAGTGCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATGCCAGAGACTGTCCGTTGAAATATTAACTACCCACTTACCGGAGGCATTTCCTCCGATAAATcc                                                                                                                                                                                                                                                                                              <  1:365080/298‑1 (MQ=255)
                     ccGTTATTGGTCAGGGTTAACGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCg                                                                                                                                                                                                                                                                                      <  1:77197/301‑1 (MQ=255)
                                   ggTTAACGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGCCATTTCCTCCGATAAATCCATAGCCCGTTCACTTAAGGTTg                                                                                                                                                                                                                                                                        >  2:374631/1‑301 (MQ=255)
                                      tAACGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGCCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGAAAACTCAATAGCCAGTTAAATTAAGGGTGGCa                                                                                                                                                                                                                                                                     >  1:617709/1‑301 (MQ=255)
                                        aCGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCAca                                                                                                                                                                                                                                                                   >  2:109651/1‑300 (MQ=255)
                                           tacccgttccAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACCGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCAtt                                                                                                                                                                                                                                                                <  2:320629/301‑1 (MQ=255)
                                                        tGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTTGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACg                                                                                                                                                                                                                                                   <  1:458691/301‑1 (MQ=255)
                                                        ctgccgtccttactttttGAATGTTTGTATAAGAGAGTTCTGTCCATCTACCCTTACGTCGCGCTTACTTTTGTCAGTGAGAGCTATCTGGGGATTACGCGCTTTAACTATACCGCCGATGGGAATATGAAGTGCACCCGTGCATAATGCACCGACATGGAGCAATTGTGCGTTGGTAGATGGTGTCTTTAAATGCCAGAGACTGTCCGTTGAAATATTAACTAGCCCCTTAGCGGAGGCATTTCCGCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACg                                                                                                                                                                                                                                                   <  1:371452/285‑1 (MQ=255)
                                                           cccTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACAATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATa                                                                                                                                                                                                                                                >  2:187563/1‑301 (MQ=255)
                                                                         tGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTATGATACCGCTACCAACAtt                                                                                                                                                                                                                                  >  1:392122/1‑301 (MQ=255)
                                                                                     aaatagagTTCTGCCCATCCACCCTCACGTCTCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGAt                                                                                                                                                                                                                      <  1:232649/297‑1 (MQ=255)
                                                                                          gagTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCca                                                                                                                                                                                                                 <  2:404953/301‑1 (MQ=255)
                                                                                           ttcttccgatctCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcac                                                                                                                                                                                                                <  2:201643/290‑1 (MQ=255)
                                                                                              tCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCAACAga                                                                                                                                                                                                             >  2:183046/1‑300 (MQ=255)
                                                                                                tGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGAAATTGATACCCACAgaaa                                                                                                                                                                                                           >  2:349894/1‑298 (MQ=255)
                                                                                                    cATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGaa                                                                                                                                                                                                       >  2:638925/1‑301 (MQ=255)
                                                                                                       ccACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGAAATTGATACCCACAgatcggaagag                                                                                                                                                                                                    >  1:201400/1‑291 (MQ=255)
                                                                                                         ccccTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGAAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatg                                                                                                                                                                                                  <  2:304779/300‑1 (MQ=255)
                                                                                                            cTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGatgatg                                                                                                                                                                                               <  1:271432/301‑1 (MQ=255)
                                                                                                             tCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGt                                                                                                                                                                                              >  1:84528/1‑301 (MQ=255)
                                                                                                                                                tgtgTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATACCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAACTTACCGACATTGATACCCACAGTATCGAACGATGATGTCGAATCCATCAGGTTAAGTGTGCAGATACCgattt                                                                                                                                                           >  2:424576/1‑300 (MQ=255)
                                                                                                                                                   ttATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCACAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTATGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgcttgctt                                                                                                                                                        <  2:392636/300‑1 (MQ=255)
                                                                                                                                                               ttAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCaaaa                                                                                                                                            >  2:514686/1‑301 (MQ=255)
                                                                                                                                                                  acaataCCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCTACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAg                                                                                                                                         <  1:187321/298‑1 (MQ=255)
                                                                                                                                                                                ggtAATATTCAGTTCACCCGTGCCTAATACACCGCCCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGAATCCATCAGGTTCAGTGTGCCGATACCACTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGaa                                                                                                                           >  2:362656/1‑300 (MQ=255)
                                                                                                                                                                                   aaTCTTCAGTTCACCCGTGCCTAATACACCGACCTGTATCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACg                                                                                                                        <  1:349431/301‑1 (MQ=255)
                                                                                                                                                                                   aaTATTCAGTTCACCCGTGCCTAATTCACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACg                                                                                                                        <  1:109527/301‑1 (MQ=255)
                                                                                                                                                                                      atTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGaaa                                                                                                                     >  2:242989/1‑300 (MQ=255)
                                                                                                                                                                                                           atatccGACCTGTGTCAATGTTGCGTTATTATATGATGTCTTTAAATTCCGGATTCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTg                                                                                                <  1:182849/297‑1 (MQ=255)
                                                                                                                                                                                                               ccGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGt                                                                                            >  2:389652/1‑301 (MQ=255)
                                                                                                                                                                                                                       tAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTacgacg                                                                                    >  2:525071/1‑301 (MQ=255)
                                                                                                                                                                                                                       atgtaagtgttCGTTTGTAGGTGATGTTTTTAAATTCATGAGAAGGTCCGTTGAAATATTAACTAACACCGTACCGGAGGCAATGCCTCCGATAAATCCGTAGGCCGTTGAATTAAGGGTGGTACCATTAATGACAGTTACGATACCGCGACCAAAATTACCGACATTGATACCCACAGTATCGACCGAGTATGGCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACg                                                                                    <  1:374135/294‑1 (MQ=255)
                                                                                                                                                                                                                        aGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGCACTTACCGCCATCGCTAATATTGAGGTTTCCGACGc                                                                                   >  1:103225/1‑301 (MQ=255)
                                                                                                                                                                                                                          aaaTTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGccc                                                                                 <  1:524399/300‑1 (MQ=255)
                                                                                                                                                                                                                                     tagtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAAt                                                                      <  2:103349/301‑1 (MQ=255)
                                                                                                                                                                                                                                      agtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCAGACGCCCGCAACACCAAtg                                                                     >  2:321864/1‑301 (MQ=255)
                                                                                                                                                                                                                                       tcttccgatcTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAAtgg                                                                    <  1:84752/295‑1 (MQ=255)
                                                                                                                                                                                                                                          gatgatGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTTACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAAtggtgg                                                                 <  1:28425/301‑1 (MQ=255)
                                                                                                                                                                                                                                           atgatgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATggtggt                                                                <  1:514014/301‑1 (MQ=255)
                                                                                                                                                                                                                                                  cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAAtggagatcggaaga                                                         >  2:84876/1‑290 (MQ=255)
                                                                                                                                                                                                                                                  cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAAATTACCGCCATCGCTAAAATTGAGGTTTCCGACGCCCGAAACACCAATGGTGGAATAGACa                                                         >  1:105762/1‑301 (MQ=255)
                                                                                                                                                                                                                                                  cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCACAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACa                                                         >  2:72474/1‑301 (MQ=255)
                                                                                                                                                                                                                                                         ttCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGCCATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATCTTGAAGTTTCCGACGCCCGCAACACCAATGGTGGAATAAACATTCTTAt                                                  >  1:388964/1‑301 (MQ=255)
                                                                                                                                                                                                                                                            cAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGAAACGCACTTCCCGCCCTCGCACATATTAAGGTTTCCGACGCCCCCACCACCAATGAGGTTTTAGACCTTCTtttttt                                               >  1:322683/1‑299 (MQ=255)
                                                                                                                                                                                                                                                               agaCTTTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCt                                            <  2:84652/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATGAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTc                                           >  1:259153/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTc                                           >  2:569630/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                 aCTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATGAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTcc                                          <  2:259421/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ttGTCCGGTGAATTATTAACTTTCCCCTTACCGGAGGCATGTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATGGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTccc                                         <  1:638048/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                  cTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCCTAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTccc                                         <  1:389138/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                      tgtgctcttccgatcttATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACAt                                     <  2:70662/285‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      tATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCcacacaaaacggaagaacgct                     >  1:70572/1‑284 (MQ=255)
                                                                                                                                                                                                                                                                                       ctCCTCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACGGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCGTGGTGATGc                    <  1:242722/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           atttATCGGAGGAAATGCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCGTGGTGATGCTGaa                >  2:293335/3‑300 (MQ=255)
                                                                                                                                                                                                                                                                                             ttATCGGAGGAAATGCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACTCCCGCAACACCAATAGAGGTATAGACATTCTTATTAACCTCCAAAATTCCCCCCGTGATGACACTCACAt              >  1:233186/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                  gggtgCATTTCCTCCGATTAATCCATAGCCCTTTGAAGTAATGTTGGCACCATTACTGACATTTGCGATACCGCTACCCAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTGAGTGTGCCGATAGCGCTTGCTTTATCGCCCAAAAAAGTAATATTGTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACAGCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCGTGGTGATGCTGACATTCCCg         <  1:362172/297‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                  ggAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTACTATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGAACTTCTTATTGACCTCCCACATTCCCCCCGTGGTGATGCTAACATTcccc         >  1:308177/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                   gAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCGTGGTGATGCTGACATTCCCGa        >  2:42710/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                      gCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCGTGGTGATGCTGACATTCCCGACGc     >  2:301837/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                       aaTTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCGTGGTGATGCTGACATTCCCGACGcc    <  1:301433/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                         tttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCGTGGTGATGCTGACATTCCCGACGCCgg  <  2:618567/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                            
CCACCTTCGACATTCAGCGTACCGTTATTCGTCAGGGTTAACGTACCCGTACCAGATGTCCCTACGTACATATTGAATGTTTCAAGAAGAGAGTTCTGCCCATCCACCCTCACGTCGCCCTTACTTTTGTCATTGAGAGCTATCTGTGTATCACGCGCTTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCCAAAAAAGTAATATTTTGCGATACGAACTTACCGCCATCGCTAATATTGAGGTTTCCGACGCCCGCAACACCAATGGTGGTATAGACATTCTTATTGACCTCCCACATTCCCCCCGTGGTGATGCTGACATTCCCGACGCCGG  >  NC_000913/1596204‑1596801

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