Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 6,349 C→T E37E (GAG→GAA yaaA ← DNA binding and peroxide stress response protein YaaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009136,3490CT96.8% 176.0 / ‑2.5 63E37E (GAG→GAAyaaADNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (30/31);  total (31/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 = 6.59e-01

GCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAATGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCG  >  NC_000913/6112‑6584
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            
cctcaacccccctcttccattctgCAAATGTTGTTGGGCAAAACCAAAACCGTCTTCGCTGAAGGTTTCGGCTGCCAAGCCGTTGTAGACATCACCTTAAACCGCCAGAATCGCCTGGCGGGAATTCCCCGGCGTAAAATTTGGCTGCCAGTCATAAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAAct                                                                                                                                                                                                                                 <  2:386061/228‑1 (MQ=255)
  cATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGAAAAATCGAAATCGTCTCCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAActgc                                                                                                                                                                                                                               <  1:659388/251‑1 (MQ=255)
   aTACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAActgct                                                                                                                                                                                                                              >  1:153083/1‑251 (MQ=255)
           cGGAAAGCATTCGCAAATGATGTTGGGCAAAATCGAAATCGTCTTAGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCCGTCATGAAGGAGAACGGCGTTGATACCCGCCAATTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAAGTTCCGCGCTTCATGGATAAACGTCTGGGAATTg                                                                                                                                                                                                                      >  2:676507/1‑251 (MQ=255)
              aaaGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTGTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCt                                                                                                                                                                                                                   <  1:166533/251‑1 (MQ=255)
               aaGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCt                                                                                                                                                                                                                   >  1:615531/1‑250 (MQ=255)
                aGCATTCGCAACTGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCAGTGTAGACATCACCTTTAAACGCAAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGGTTACGCGCTTCATGGATCAACTGCTGGGAATTGTCTa                                                                                                                                                                                                                  >  1:372901/1‑250 (MQ=255)
                   aTTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAATCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAAc                                                                                                                                                                                                                >  1:631213/1‑249 (MQ=255)
                      agCAAATGCTGTTAGGCAAAATGGAAATAGTCTTCGCTGAAGGTCTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGACAAAATCGCATGGCGGGCATTCTACGGCGTGACATATGGCTGCAAGTCATGAAAGAGAGGGACGTTGATACCCACAATATTGACGCTGATGCGGAACAACGTGCATATATGAGGAGGAGTCGGGTTACTCGCTTCATGGCACAACGGCTGGGAATGTTCTAACAGCTc                                                                                                                                                                                                           >  2:652653/2‑251 (MQ=255)
                       gCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAAct                                                                                                                                                                                                                                 >  1:386061/1‑228 (MQ=255)
                       gCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTcc                                                                                                                                                                                                          >  1:396470/1‑251 (MQ=255)
                            tGCTGTTGGGCAAAATCGAACTCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGACAGCGAGCGGCGTTGATACCCGCCAGTTTTTCGCTGATGCGCATCCGAGTGCTACACTGCGGAGGCGTAAGTTTCAGCGCTTCCTGGATCAACTGATGGGAATTGTTTACCAGCTCCGGccg                                                                                                                                                                                                     >  2:634000/1‑249 (MQ=255)
                                      cAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAgcgc                                                                                                                                                                                           >  1:466424/1‑251 (MQ=255)
                                       aaaaTCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGcgcg                                                                                                                                                                                          >  1:494217/1‑251 (MQ=255)
                                           ccgaaatcgTCTCAGCTGAAGGTTCCGGCCTGAAAGCCGGTGTAGACATCACTTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCgtggt                                                                                                                                                                                      <  1:664721/250‑1 (MQ=255)
                                              aaatcgTCGTCACTGAAGGTTTCGGACTGCAAGCCGGTGTAGACATCACCTTTAACCGCCAGAATCGCCTGGCGGGGATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGGGCGGCGTTGATACCCGCCAGTTTGGCGCTGATGCGCATCAGCGTGCTAATATGCGGAGGCGTGAGTTTCCGAGCTTCAAGGATCAAATGAGGGGCATTATCTGACAGCTCGGGCAGCGTAAAGCGCGtggtggg                                                                                                                                                                                   >  2:333223/1‑250 (MQ=255)
                                                       tCGCGGAAGGTTTCGGCCTGCAAGCCGGTGTAGCCCTCACCTTTAAACGCCAGAATCGCCTGGCGGCCATTCTCCGGCGTGAAATCGGGCTGCAAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACgg                                                                                                                                                                            <  1:373925/248‑1 (MQ=255)
                                                         gCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCttt                                                                                                                                                                       >  1:434204/1‑251 (MQ=255)
                                                             atggTTGCGGCCTGCAAGCCGTCGGAGAAACCCCCTTAAAACGCCAGAATTGCGTGCCGGGCACTCGCCGGCGATAAGTATGGCTGCATGTCATGAAAGCGAGTGTGGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGGGGTGGTCAACGTGCTTTGGTa                                                                                                                                                                   <  2:631213/249‑1 (MQ=255)
                                                                 tttCGGCCTGCAGGCCGGTGGAGACATCACCTTTAAACGCAAGAACCGCCTGGCGGGCATGCTCCGGCGTGAAATCGGGCTGCCGGTCATGAAAGCGAGGGGGGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATGGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATc                                                                                                                                                                <  1:305372/250‑1 (MQ=255)
                                                                     cggctggaaaGCCGGTGTAGACATCACCTTAAAACGCCATAAGCCTCTGGCGGGCATTCTCCGGTGTGAAATCTGGCTGCCCGTCATGAAAGCGGGCGGGGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGc                                                                                                                                                            <  1:126040/243‑1 (MQ=255)
                                                                                ccGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAgg                                                                                                                                                 >  1:109737/1‑251 (MQ=255)
                                                                                                                  ccgggCGGGCATTCTCTGGCGGTAAATCTGGCTGCCAGTCAGGAAAGCGAGCGGTGTTGATACCCGCCAGTTTGGCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAgg                                                                                                               <  1:634000/248‑1 (MQ=255)
                                                                                                                     ggCGGGCATTCGCCGGCGTGAAGTCGGCTCGCCAGTCATGAAAGCGGGGGGCGTTTACACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCCGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGa                                                                                                              <  1:263979/249‑1 (MQ=255)
                                                                                                                      gggtgctttccccGGCGTGAAATCTGGCTGCCAGCCAGAAAACCGAGCGGGTTTAATCCCCGCCAGGTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGaaaa                                                                                                           <  2:294486/240‑2 (MQ=255)
                                                                                                                        gggCATTCCCCGGCGGAAAATCTGGCTGCCAGTTATAAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAAttt                                                                                                         <  1:518608/251‑1 (MQ=255)
                                                                                                                             ttctccGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGTTCAACGGGATTTGGTAATCAAGCGATTTCGCAGGTGAATTAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGcc                                                                                                    >  2:652747/5‑251 (MQ=255)
                                                                                                                                      gggaatattgggtggcaggTCAGGAAGGCGAGAGGTGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAg                                                                                            <  2:153083/233‑1 (MQ=255)
                                                                                                                                      gTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGaa                                                                                                                                      >  1:298576/1‑208 (MQ=255)
                                                                                                                                      gTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGc                                                                                           >  1:452925/1‑251 (MQ=255)
                                                                                                                                       tGAAATCAGGCTGCCAGACATGAAAGCGAGCGGCGTTGATACCCGCCAGATTGTCGCTGCTGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGATGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGGTTTCGCAGGTGAAATAAGAATCAGCCTATCCAGTCCTTGCAGGAAATTAATGCGGACTTTAGc                                                                                           >  2:585904/1‑250 (MQ=255)
                                                                                                                                          aatgctggtggccagctaTGAAAGCGAGGGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGGCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCaaa                                                                                        <  1:674339/234‑1 (MQ=255)
                                                                                                                                                  tGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCAGATGCGCATAAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGGCAACGGGCCTTGGTAATCAAGCGTTTGAGCAGGTGAAATAAGAATCAGCATATCCAGACCTTGCAGGAAATTTATGCCGACTTTAGCAGAAaaagaga                                                                                >  2:297440/1‑250 (MQ=255)
                                                                                                                                                          aaTGAAAGCGAGCGGCGTGGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGtt                                                                         <  1:469112/249‑1 (MQ=255)
                                                                                                                                                              aaaGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGAGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCg                                                                    >  1:334178/1‑251 (MQ=255)
                                                                                                                                                                agcgagcgGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGAt                                                                  <  1:144626/251‑1 (MQ=255)
                                                                                                                                                                   gagcgGCGTTGATAAACGCCAGTTTGTAGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCAAGTCCTTGCAGGAAATTTATGCCGCCTTTAGCAAAAAAAGAGAATGAGTTGATCGaacgt                                                               >  2:609270/1‑247 (MQ=255)
                                                                                                                                                                     gcgggttttatccccGCCAGGTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTTCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTtg                                                             <  2:434204/239‑1 (MQ=255)
                                                                                                                                                                       ggcgTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTtgtg                                                           <  1:454508/251‑1 (MQ=255)
                                                                                                                                                                       ggcgTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTtgtg                                                           >  1:395495/1‑251 (MQ=255)
                                                                                                                                                                         ggttgacccccGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGAt                                                         <  1:609999/243‑1 (MQ=255)
                                                                                                                                                                            tGATCCCCGCCAGGTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTGCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTa                                                       <  1:465499/250‑1 (MQ=255)
                                                                                                                                                                                   cGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGTAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACgg                                                                                                                                                                            >  2:15309/1‑126 (MQ=255)
                                                                                                                                                                                   cGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACgg                                                                                                                                                                            <  1:15309/126‑1 (MQ=255)
                                                                                                                                                                                   cGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCg                                               <  1:385562/251‑1 (MQ=255)
                                                                                                                                                                                        gTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAAc                                           >  1:569042/1‑250 (MQ=255)
                                                                                                                                                                                         tttGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCTTATAGCGCGTGGTGGTCAACGGGCTTTGTTAATCACGCGTTTTCGAAGGTGAAATAAGAATCAGCATGTCCAGTCCTTGCAGGAAATTTATGCCGACTTTAACAAAACAAGAGAATGAGTTGATCGATAATTGTGATTACTCCTGCGAAAca                                          >  2:26317/1‑250 (MQ=255)
                                                                                                                                                                                          gtgttgcttaTGCTCATCAGCGTGGTAATCTGCGGAGGCGTCAGTCTCCGCGCTTTATGGAACAACGGGTGTGAATTGTCGAACAGCTCCGGCAGCGTATAGCGCGGGTGGGTCAACGGGCTTTGGTAATCAAGCGTTTGCGCATGTTAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGGCTTTAGCAAAAAAAGAGAATGAGTTTTTCGATAGTTGTTATTACTCCTgttgaacat                                         <  2:620178/241‑9 (MQ=255)
                                                                                                                                                                                            gTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAAcatcat                                      >  2:214282/1‑251 (MQ=255)
                                                                                                                                                                                            gTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAAcatcat                                      <  1:345334/251‑1 (MQ=255)
                                                                                                                                                                                               gCTGAGGCGCATCAGCGTGCTAAGCTGTGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTAGAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATTAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATccc                                   <  2:580282/251‑1 (MQ=255)
                                                                                                                                                                                                   aTGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATACACAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCAcgc                               >  1:583284/1‑251 (MQ=255)
                                                                                                                                                                                                        cATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGTGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCGTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCg                          <  1:609270/251‑1 (MQ=255)
                                                                                                                                                                                                         aTCAGCGTGCCAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCgg                         >  1:109994/1‑251 (MQ=255)
                                                                                                                                                                                                           cAGCGTGCTAATCATCGGAGTCGTCAGTTTCCGCGTTTCAGGGATCAACTGCTGGGATTTGTCTAACAGCTCCGGCCTCTTATAGCGCGTGGTGGGCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATTACCAGTCCTTGCATGAAATTTATGCCGACTTTTGCAAAAAAAGATAATTAGTTGCTCGATAGTTGTTATTATTTCGGTGAAACATCATCCCACGCGTCCGgag                       <  2:445188/251‑1 (MQ=255)
                                                                                                                                                                                                                ggCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCAAGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGc                  <  1:26317/250‑1 (MQ=255)
                                                                                                                                                                                                                            gAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTACTCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGat       >  1:281687/1‑250 (MQ=255)
                                                                                                                                                                                                                             aGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGTGCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATAGCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCTAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGatat     <  2:494217/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATc    <  1:652747/251‑1 (MQ=255)
                                                                                                                                                                                                                               gCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGCTAGTTGTGATTACTCCTGCGAAACCTCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATcc   >  2:596861/1‑251 (MQ=255)
                                                                                                                                                                                                                               gCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGAGCGATAGTTGTGATTACTCCTTCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATcc   <  2:566070/251‑1 (MQ=255)
                                                                                                                                                                                                                                cGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCg  >  1:80741/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGc                                                                                           <  1:139123/160‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            
GCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAATGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCG  >  NC_000913/6112‑6584

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