Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,214,955 C→A 6.8% intergenic (+67/+12) rpoD → / ← mug RNA polymerase, sigma 70 (sigma D) factor/G/U mismatch‑specific DNA glycosylase; xanthine DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,214,9550CA 6.8% 237.9 / 6.0 73intergenic (+67/+12)rpoD/mugRNA polymerase, sigma 70 (sigma D) factor/G/U mismatch‑specific DNA glycosylase; xanthine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base C (35/33);  new base A (2/3);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.56e-01

CGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTG  >  NC_000913/3214808‑3215098
                                                                                                                                                   |                                                                                                                                               
cGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTc                                                                                                                                                 <  3:218696/148‑1 (MQ=255)
cGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTGTCg                                                                                                                                                >  1:53077/1‑149 (MQ=255)
   cAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTg                                                                                                                                             <  2:237757/149‑1 (MQ=255)
     gATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGc                                                                                                                                           >  6:311607/1‑149 (MQ=255)
       tCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCg                                                                                                                                                >  5:357014/1‑142 (MQ=255)
       tCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCg                                                                                                                                                <  6:357014/142‑1 (MQ=255)
         gaagcgaagGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGttt                                                                                                                                       <  6:3679/149‑1 (MQ=255)
          cagcgaagGCGCTGCGCAAACTGCGGCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGAGCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGtttt                                                                                                                                      <  6:282788/148‑1 (MQ=255)
          aagcgaagGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGtttt                                                                                                                                      <  4:119939/149‑1 (MQ=255)
           agcgaagGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCGAAAAGGTGGCGtttgt                                                                                                                                     <  4:118822/149‑3 (MQ=255)
            gcgaagGCGCTGCGCAAACTGCGGCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGtttttt                                                                                                                                    <  1:128677/149‑1 (MQ=255)
             cgaagGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTa                                                                                                                                   <  4:121962/149‑1 (MQ=255)
              gaagGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTAt                                                                                                                                  >  7:2978/1‑149 (MQ=255)
              gaagGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTAt                                                                                                                                  >  5:88321/1‑149 (MQ=255)
               aagGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATc                                                                                                                                 >  3:155085/1‑149 (MQ=255)
               aagGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCGAAAAGGTGGCGttt                                                                                                                                       <  6:337192/143‑1 (MQ=255)
               aagGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCGAAAAGGTGGCGttt                                                                                                                                       >  5:337192/1‑143 (MQ=255)
                    gcTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCa                                                                                                                            >  3:81401/1‑149 (MQ=255)
                    gcTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCGAAAAGGTGGCGtttgt                                                                                                                                     >  6:225625/1‑138 (MQ=255)
                    gcTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCGAAAAGGTGGCGtttgt                                                                                                                                     <  5:225625/140‑3 (MQ=255)
                           aaaCTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCCCGCACTa                                                                                                                     >  2:335689/1‑149 (MQ=255)
                            aaCTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTAc                                                                                                                    >  4:142799/1‑149 (MQ=255)
                                  gTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAgcgc                                                                                                              >  1:72983/1‑149 (MQ=255)
                                  gTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGACGCAAACGGAACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCTTTTTTTATCGACCACTCACTACCAgcgc                                                                                                              >  4:265478/1‑149 (MQ=255)
                                     aCCCTAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTg                                                                                                           <  5:65555/149‑1 (MQ=255)
                                      cccGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTgg                                                                                                          <  4:81401/149‑1 (MQ=255)
                                          aGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCa                                                                                                      <  8:205544/149‑1 (MQ=255)
                                          aGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCa                                                                                                      <  2:271990/149‑1 (MQ=255)
                                                 cTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTcgc                                                                                                <  5:211204/148‑1 (MQ=255)
                                                   gAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGa                                                                                              >  5:339188/1‑148 (MQ=255)
                                                       tgctgcGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGc                                                                                         >  7:127105/1‑149 (MQ=255)
                                                         ctgcGTAGCTGCCTGGACGATTAATCGGTAGGCAGGAGCAGGAGTTACGCCGCACCCGGCACTAGGCCCGCTGCACAAACGCCACCTTTTCGGTGGAGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCtt                                                                                       <  3:127077/149‑1 (MQ=255)
                                                           gcGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCa                                                                                     >  8:203975/1‑149 (MQ=255)
                                                                      tGGCCGATTAATCGGTAGGCCGGATCAGGCGTGACGCCGCACCCGGCCCTAGGCCCTCGGCACAAACGCCACCTTTTCGGTGGCGTGGTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTtctc                                                                          >  4:303933/1‑149 (MQ=255)
                                                                      tGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTtctc                                                                          >  4:240158/1‑149 (MQ=255)
                                                                        gACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCa                                                                        <  8:2978/149‑1 (MQ=255)
                                                                        gACGATTAATCGGTAGGCAGGATCAGGCGTTACGCCGCACCCGTCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGTCGTTTTTTATCGCCCACGCACTACCATCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCa                                                                        <  3:265478/149‑1 (MQ=255)
                                                                         aCGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAg                                                                       <  5:311607/149‑1 (MQ=255)
                                                                          cGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGt                                                                      <  2:72983/149‑1 (MQ=255)
                                                                             ttAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGtt                                                                              >  7:190804/1‑138 (MQ=255)
                                                                             ttAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGtt                                                                              <  8:190804/138‑1 (MQ=255)
                                                                               aaTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGaaa                                                                  <  5:158008/148‑1 (MQ=255)
                                                                                aTCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGATCGCGATACGCTTCAACCAGTTTCTCCAGTGAAAcg                                                                >  8:202787/1‑149 (MQ=255)
                                                                                   ggTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGa                                                             >  4:197045/1‑149 (MQ=255)
                                                                                    gTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGCTACGCTACAACCAGTTTCTCCAGTGAAACGCGAc                                                            >  6:364414/1‑149 (MQ=255)
                                                                                         ccGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTaaa                                                       <  1:13508/149‑1 (MQ=255)
                                                                                         ccGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTaaa                                                       >  4:332837/1‑149 (MQ=255)
                                                                                                aGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTgg                                                >  8:347489/1‑149 (MQ=255)
                                                                                                   cGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGAt                                             <  2:92771/149‑1 (MQ=255)
                                                                                                     ttACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGAttt                                           >  1:326383/1‑149 (MQ=255)
                                                                                                         gccgcACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTgg                                         >  2:207394/1‑147 (MQ=255)
                                                                                                         gccgcACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGgca                                       <  7:203975/149‑1 (MQ=255)
                                                                                                           cgcACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGgcag                                      >  4:61927/1‑148 (MQ=255)
                                                                                                              aCCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTgggtgctgc                                  <  6:249997/149‑8 (MQ=255)
                                                                                                                         ggCCCTCTGCACAAACGCCACGTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGc                         <  3:332837/147‑1 (MQ=255)
                                                                                                                         ggCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGTAGCACCCAAATCTGCGt                       >  7:177687/1‑149 (MQ=255)
                                                                                                                         ggCCCTCTGCACAAACGCCACCTTTTCGGGGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATCCGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGt                       <  6:272500/149‑1 (MQ=255)
                                                                                                                           cccTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCg                     <  3:197045/149‑1 (MQ=255)
                                                                                                                             ctctGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTg                                                                     >  4:94683/1‑99 (MQ=255)
                                                                                                                             ctctGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTg                                                                     <  3:94683/99‑1 (MQ=255)
                                                                                                                              tctGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAAc                  >  6:53920/1‑149 (MQ=255)
                                                                                                                                tGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAAcc                 >  7:223811/1‑148 (MQ=255)
                                                                                                                                  cacaAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCaa               <  7:347489/148‑1 (MQ=255)
                                                                                                                                  cacaAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAAt              >  7:202388/1‑149 (MQ=255)
                                                                                                                                  cacaAACGCCACCTTTTCGGTGGCCTTTTTTATCGCCCACGCACTACCAGCTCCTGGTCCAGCTCGCGATAAGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTCAGTCGACCCAAt              >  6:99192/1‑149 (MQ=255)
                                                                                                                                   acaAACGCCACCTTTTCGGGGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATg             >  3:322046/1‑149 (MQ=255)
                                                                                                                                     aaaCGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGt           >  5:72140/1‑149 (MQ=255)
                                                                                                                                     aaaCGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGt           <  3:142799/149‑1 (MQ=255)
                                                                                                                                       aCGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGa         <  6:88321/149‑1 (MQ=255)
                                                                                                                                       aCGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGa         <  1:207394/149‑1 (MQ=255)
                                                                                                                                          ccACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCg      >  7:128922/1‑149 (MQ=255)
                                                                                                                                             ccTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGTAGCACCCAAATCTGCGTCGAACCAATGGTGAGCgttt   <  8:177687/149‑1 (MQ=255)
                                                                                                                                             ccTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCgttt   <  8:223811/149‑1 (MQ=255)
                                                                                                                                              cTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCgtttg  >  2:4076/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CGTCAGATCGAAGCGAAGGCGCTGCGCAAACTGCGTCACCCGAGCCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCACTACCAGCGCCTGGTCCAGCTCGCGATACGCTTCAACCAGTTTCTCCAGTGAAACGCGACTTAAACCGCTGGGATTTGGCAGCACCCAAATCTGCGTCGAACCAATGGTGAGCGTTTG  >  NC_000913/3214808‑3215098

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