Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 6,106 G→C L118L (CTC→CTG 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,1060GC96.9% 185.1 / ‑0.4 64L118L (CTC→CTGyaaADNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base C (37/25);  total (38/26)
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.68e-01

CATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCT  >  NC_000913/5874‑6350
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       
cATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCg                                                                                                                                                                                                                                       >  1:549015/1‑251 (MQ=255)
       gAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTAc                                                                                                                                                                                                                                                    >  2:150882/1‑231 (MQ=255)
       gAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTAc                                                                                                                                                                                                                                                    <  1:150882/231‑1 (MQ=255)
          cacCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTc                                                                                                                                                                                                                            <  1:46480/251‑1 (MQ=255)
                 ttGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATg                                                                                                                                                                                                                     >  1:489368/1‑251 (MQ=255)
                      aaTCTCGGGATTCAATTGCTTCGGCTTACCCGATTAAAAATACCAATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCGTGTTGGTGATGATATAGCCCCAGAATAGATACAGATCTTTCCCGCGGGCAGTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAAACGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGtt                                                                                                                                                                                                                <  2:350470/251‑1 (MQ=255)
                           cGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCa                                                                                                                                                                                                           <  2:286857/251‑1 (MQ=255)
                              cATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGCAGCCAGATTAATCACCACATTATAGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGCGCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGATGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGATAGCAGTCGCAAATGATGTTGGGCaaa                                                                                                                                                                                                         >  2:469112/1‑250 (MQ=255)
                                ttCATTTTCTTCTGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAACCACCACATTATCGCCTTGTGCTTCGATCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATCCGCAAATGCTGTTGGGCAaaat                                                                                                                                                                                                       <  2:602221/250‑1 (MQ=255)
                                       tcttaGGCCTCACCGATTTCCAATACTCATCTGACGCCAGATTACTCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGCTCCCCATTTCCAGACGATCAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCAGGCACGCATTCGCAAATGCTGTTGGGCAaaatcgaaaccg                                                                                                                                                                                               >  2:674339/1‑248 (MQ=255)
                                          tcGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTc                                                                                                                                                                                             <  1:659312/250‑1 (MQ=255)
                                          ccggtttCACCGATTTAAAACACCCACCTGACGCCAGATTAATCACCTCATTATCGCCTGGTGCTGCGCGCGCCTCGTTCAGCTTGTTTGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTTTTGTGCAAAATCGAAATCGTc                                                                                                                                                                                             <  2:320074/245‑1 (MQ=255)
                                                    cGATTTAAAATTCTCATCTGTCGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGATCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGTCTAAATCGAAATCGTTTGCGCTGAAgg                                                                                                                                                                                  <  2:483204/251‑1 (MQ=255)
                                                    cGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCAGGAAAGCATTCGCCAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAgg                                                                                                                                                                                  >  1:620178/1‑251 (MQ=255)
                                                          aaaaTACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTTATCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGTCAAAATCGAAATAGTCTTCGCTGAAGGTTTCgg                                                                                                                                                                            >  2:139157/1‑251 (MQ=255)
                                                             aTACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGATGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCAATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAAAAGTCTTCGCTGAAGGTTTAGGCCt                                                                                                                                                                         >  2:78756/1‑251 (MQ=255)
                                                               ccTCATCTTACGCCATATTAATACTATCATTATCGTCTTGTGCTGCGAGCGCCTCGTTCTGCTTGGTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTGCAAGACGAGAAGGCTGCATTAAATCGATCGGGCGCAGTACGCCATACAAGCCGGAAATCATTCGCAAATGCTGTTGGTCAAAATCGAAATCGTCGTCGCTGAAGGTTTCGGCCTGc                                                                                                                                                                       <  2:361146/250‑1 (MQ=255)
                                                                       gACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTCCAAGCCGgt                                                                                                                                                               <  1:139157/251‑1 (MQ=255)
                                                                              gATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGAc                                                                                                                                                         >  1:566070/1‑250 (MQ=255)
                                                                               aTTAATCACCCCATTATCGCCTTGTGCTGCGAGCGCCCCGTTCAGCTTGTTGGTTATGATATCTCCCCAGAATTGATACAGATCTTTCCCCCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACAt                                                                                                                                                       <  2:549015/251‑1 (MQ=255)
                                                                                  aaTCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCAc                                                                                                                                                    >  2:166533/1‑251 (MQ=255)
                                                                                  aaTCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTACAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTTTAGACATCAc                                                                                                                                                    >  1:250342/1‑251 (MQ=255)
                                                                                     caccacATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTATTCGCTAAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCtt                                                                                                                                                 >  2:518608/1‑251 (MQ=255)
                                                                                      accacATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGACTTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTATACATCACCttt                                                                                                                                                >  2:659388/1‑251 (MQ=255)
                                                                                        cacaTTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCACTTTCCAGACGATAAGGCTGCCTTAAATCGCGCGGGCGCAGTACGCCATACAAGCCGGAAAGCCTTCGCAAACGCTGTTGGGCACACTGGACATCGTCTTCCCTGCAGCTTTCGCCCTGCAAGCCGGTGAAGaactcaccatta                                                                                                                                               >  2:184134/1‑239 (MQ=255)
                                                                                           caTTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAgg                                                                                                                                                                                  >  2:355389/1‑213 (MQ=255)
                                                                                           caTTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAgg                                                                                                                                                                                  <  1:355389/213‑1 (MQ=255)
                                                                                            aTTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGCTCTTTACCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAACTCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAAcc                                                                                                                                           >  1:445188/1‑250 (MQ=255)
                                                                                                 cGCCTTTTGCCGCGAGCGCCTCTTTCATCTTGGTGGTGATGATATCGCCCCAGAATTGATTCAGATATTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCATACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAATTCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTATACGCCAGa                                                                                                                                      <  2:469103/251‑1 (MQ=255)
                                                                                                        ggctgctAGCGCCTCGTTCAGCTTTGTGGGGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCTATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCAGTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCt                                                                                                                               <  2:598658/250‑1 (MQ=255)
                                                                                                             cGAGGGCCTCGGTCAGATTGTTGGTGATGATATCTCCCCAGAATTGATACAAATCTGTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCgg                                                                                                                          <  2:32262/251‑1 (MQ=255)
                                                                                                                 cgcCTGGTTCAGCTTGTTGGTGATGATCTCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCTCTGAATGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCTGAATCGCCTGGCTGGCAt                                                                                                                      <  2:467893/251‑1 (MQ=255)
                                                                                                                 cgcCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGcct                                                                                                                      >  2:30001/1‑249 (MQ=255)
                                                                                                                 cgcCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGATGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTACACTCCAGCATCGCCTGGCGGTCAt                                                                                                                      >  2:144626/1‑251 (MQ=255)
                                                                                                                  gcCTCGTTCAGCTTGTTGGTCATGCTATCTCCCCAGAATTGATACAGATCTTTACCTCGGGCATTCTCAAGACGGCTCCCCATTTCCAGACGCTAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAACCGGAAAGCATTCCCAAATGCTGTTGGGCAAAAACGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCACGCCGGAGTAAGCATCACCTTTAAACGCCCGAATCGCCTGGCGGGcctt                                                                                                                     >  2:347417/1‑248 (MQ=255)
                                                                                                                     tCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAACCGCCAGAATCGCCTGGCGGGCATTctc                                                                                                                  >  2:454508/1‑249 (MQ=255)
                                                                                                                       gTTCAGTTTGTTGGTGATGATCTCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCTTTCACAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAAACGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTctcc                                                                                                                 <  2:319862/250‑4 (MQ=255)
                                                                                                                        ttCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGACGGAAAGCATTCGCACATGCTGTTGGGCAACATCGAAATCGTCTTCGCTTAAGGTTTCGGCCTGCACGACGGTGTAGCCATCCCCTTTCAACGCCAGAATCGCCTGGCGGGCAttatccgg                                                                                                               >  2:174764/1‑245 (MQ=255)
                                                                                                                                  tgGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGGGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCCTTCTCCGGCTTGGAATCt                                                                                                     >  2:609999/1‑251 (MQ=255)
                                                                                                                                    gTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTgg                                                                                                   >  1:580282/1‑251 (MQ=255)
                                                                                                                                      gctgatATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTTGCt                                                                                                 >  2:45848/3‑251 (MQ=255)
                                                                                                                                          atatCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCa                                                                                             <  1:78756/251‑1 (MQ=255)
                                                                                                                                             tctcCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTc                                                                                          >  1:443623/1‑251 (MQ=255)
                                                                                                                                             tctcCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTc                                                                                          >  1:83315/1‑251 (MQ=255)
                                                                                                                                             tctcCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGACGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCAGGCTGCCCGTc                                                                                          >  2:385562/1‑251 (MQ=255)
                                                                                                                                               cccccAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCAt                                                                                        <  1:497737/250‑1 (MQ=255)
                                                                                                                                                ccccAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCa                                                                                         <  1:30001/249‑1 (MQ=255)
                                                                                                                                                       aTTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAagcg                                                                                 >  1:463931/1‑250 (MQ=255)
                                                                                                                                                         tGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAagcgagc                                                                              >  1:145433/1‑251 (MQ=255)
                                                                                                                                                           aTACAGATCTTTCCCTCGGGCATTCGTAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCTGGAAAGCATTCGCAAATGCGGTTGGGCAAAATCGAAATCGTCTTCGCGGAAGGTTTCGGTCTGCAAGCCGGTGGAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCGGGCTGCCAGTCATGAAAGCGAgcgg                                                                            <  2:618537/251‑1 (MQ=255)
                                                                                                                                                             aCAGATCTTTCCCTCGGGCATTCTCAAGACGGCTCCCCATTTCCAGACGATATGGCTGCCTTAAATCGAGCGGGCGCAGTCCGCCATCCACGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGCAATCGTCTTCGATGAAGGTTTCGGCCTGCACGCCGGTGTAGCCATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCCGTAAGGAAAGCGAgcggcg                                                                          >  2:305372/1‑251 (MQ=255)
                                                                                                                                                                     ttCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACAt                                                                                                                                                       >  1:636339/1‑166 (MQ=255)
                                                                                                                                                                     ttCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTg                                                                      <  1:49905/246‑1 (MQ=255)
                                                                                                                                                                     ttCCCCCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACAt                                                                                                                                                       <  2:636339/166‑1 (MQ=255)
                                                                                                                                                                      ccccTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGGCATGAAAGCGAGCGGCGTTGATAccc                                                                <  2:114642/250‑1 (MQ=255)
                                                                                                                                                                                             aTCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCAt                                         >  1:294486/1‑251 (MQ=255)
                                                                                                                                                                                              tCCCCATTTCCAGACGATAAGGATGCCTTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGAAAATGCTGTTGGGCACAATCGAAATCGTCTTCGATGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCCTGCAGTCTGGCTGCCAGTCCTGCAAGGGGGCGGCGTTTGAACCCGCCAGtttttggcttatgagcat                                         >  2:436236/1‑235 (MQ=255)
                                                                                                                                                                                                  cattgcaaGACGATAAGGCGGCATAAAATCGAGCGGGCGCAGTACCCCATACAAGCCGGAAAGTATTCGCAAATGCGGTTGGGCAAAAGCGAAATCGTCTGCGCGTAAGTTTTCGGCCTGCAAGCCGGTGTATACATCACCTTTAAACGCCATAATCGCCTGGCGGGCATTCTCCGTCGTGAAATCTGGCGGCCAGTCAGGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCg                                    <  1:436236/244‑1 (MQ=255)
                                                                                                                                                                                                   aTTTCCAGACGATAAGAATGCATTAAATCGAGCGGGCGCAGTACTCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGAAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTACGTAGTGTAATCTGGCTGCCAGTCATTAAAGCGCGCGGCATTGATACCCGCCAGTTTGACGCTGATGCGCATccgatt                                   >  2:44535/1‑246 (MQ=255)
                                                                                                                                                                                                           gcaaTAAGGCAGAATTAAAGCGAGCGGGCGCAGTACCCCAACCAGGCCGGAAGCCTTGCGAAAATGCTGTTGGGAAAAATCGAAATCTTCTTCGCTGAGGGTTCCGGCCGGAAAGCCGGTGTAGACATCACGTTTAACCGCCAGAATCGCCTGGCGGGCATTCCCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGGGAGGGGCGTTAATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCt                          <  2:145433/248‑1 (MQ=255)
                                                                                                                                                                                                                       cATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCAACTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGAAGCGCATCAGCGTGCTCATCTGCGGCGGCGTc               >  1:99905/1‑250 (MQ=255)
                                                                                                                                                                                                                            aaTCGAACGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACACCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGGCATGAAAGAGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAACCTGCGGAGGCGTCAGTTTc         >  1:291768/1‑251 (MQ=255)
                                                                                                                                                                                                                                gagtgggcgcAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCTGCGttt                                                                       <  2:636121/183‑2 (MQ=255)
                                                                                                                                                                                                                                gAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTg                                                                      >  1:636121/1‑187 (MQ=255)
                                                                                                                                                                                                                                aaGCGGACGCAGTACGCCATACAAAACGGAAAGCATTCGCAAATGCTGTTGGGCACAATCGAAATCGTCTTCGCTGAAGGTTTAGGCCTGCAAGCCGGTGTAGACATCACCTTTACACGCCAGAATCGCCTGGCGGGCCTTCTCCGGCGTGAAATCTGGCGGCACGTCCAGAAAGCTCGCGGATTTGATCCCCGCAAGTTTGTGTCTTGTGCGATAAGCATGGCTAATCTGCGGAGGGGTCAGTTTccacg     >  2:501103/2‑248 (MQ=255)
                                                                                                                                                                                                                                   ggggcgcAGTCCGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCgcgctt  <  1:45848/250‑3 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       
CATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCT  >  NC_000913/5874‑6350

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