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,1060GC98.3% 174.6 / ‑3.5 58L118L (CTC→CTGyaaADNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (32/25);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

CCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGC  >  NC_000913/5864‑6348
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
ccGTTCTTCACATCGAGGAACACCGGCTTGATAATCGCGGCATTCAATTTCTGCGGCTTCACCGATTTAAAATCCTCATCTGACGCCAGATTAAGCACCACATTATCGCCTTGGGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGcc                                                                                                                                                                                                                                            <  2:14356/251‑1 (MQ=255)
     cttctCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAGATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATAGAGAGGGAGCAGTACGCCAtgc                                                                                                                                                                                                                                        >  1:223002/1‑248 (MQ=255)
          cATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCGTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCg                                                                                                                                                                                                                                  <  2:104413/251‑1 (MQ=255)
                            tGATAATCTCGGCATTCAATTTCGTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTTGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCATTACGCCATACAAGCCGGAAAGCATTCGCAAATg                                                                                                                                                                                                                 <  1:412905/250‑1 (MQ=255)
                                      ggCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACACGCCGGAAAGCATTCGCAAATGCTGTTGGGCa                                                                                                                                                                                                       >  1:67982/1‑250 (MQ=255)
                                         aTTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAaaat                                                                                                                                                                                                   >  1:227587/1‑251 (MQ=255)
                                               tttCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAaaatcgaaat                                                                                                                                                                                             >  1:140159/1‑251 (MQ=255)
                                                     cGGCTTCATCGATTTAAAATACTCAACTGCCTCCATATTAATCACCAGATTACCGCCTTGTGCTGCGAGCTCCTCGTTCAGCTTGTTGTTGATGAGATCTCCCCAGATTTGATACAGATCATTCCTTCGGGCATGCTCAAGACGGACCCCCATTTCCTGACGATAAGGCTGCATTAAATCGAGCTGGCGCAGTACGCCATACAAGTCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCtt                                                                                                                                                                                       <  2:375654/251‑1 (MQ=255)
                                                     cGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCt                                                                                                                                                                                        <  1:299578/250‑1 (MQ=255)
                                                         ttCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCt                                                                                                                                                                                   >  1:108757/1‑251 (MQ=255)
                                                            cccGATTTAAAATACTCATCTGACAGCAGATTAATCCCCACATTAGCGCCTTGAACTGCGACCACCGCGTTCAGCTTGTTGGTGATCATATCTCCCCAGAAATGTGACAGAACTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGACTGCATTAAATCGAGCGGGCGCATTACGCCATACAAGACGGAAAAAATACGCAAATGCTGTTGAGCTAAATCGCAAGCGTCTTCGCTGat                                                                                                                                                                                >  2:196885/2‑250 (MQ=255)
                                                              cGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTCAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAgg                                                                                                                                                                              >  1:266782/1‑251 (MQ=255)
                                                                aTTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGtt                                                                                                                                                                            <  2:370771/251‑1 (MQ=255)
                                                                  ttaaaaatcccATCTGACGCCAGAGTAATCACTACAGTATCGCTGTGGGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGTTGATATCTCCCCAGAATTGATACATATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAAACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAAGGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGttt                                                                                                                                                                           <  2:169352/240‑1 (MQ=255)
                                                                     aaaTACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGc                                                                                                                                                                       >  1:102995/1‑251 (MQ=255)
                                                                      aaTACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAACATCGAAATCGTCTTCGCTGAAGGTTTCGGcc                                                                                                                                                                      >  1:467607/1‑251 (MQ=255)
                                                                        aaCTCATCTGTAGCCAGATTAATCACCACATTATAGCCTCGTGCTGCGTGCGCTTCGTCCATCTTGTTGGTGATGATATCTCCCCAGAATTTATACAGATCTTTCCCTCGGACATTCTCAAGACGGATCCCCATTTCCAGACTATAATGCTGCATTAAAACGAGCGGGCGCAGTACGCCATACAAGCCGGAATGCATTCGCAAATGCTGTTGGGCAAAATCGAAACCGTCTTCTCTGAAGGTTTCTGCCTg                                                                                                                                                                    <  2:156467/250‑1 (MQ=255)
                                                                            cATCTGACGCCAGATTAAACACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTGTCGGCCTGCAAg                                                                                                                                                                <  2:181019/251‑1 (MQ=255)
                                                                                    gCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGgtgt                                                                                                                                                         >  1:381009/1‑250 (MQ=255)
                                                                                    gCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTa                                                                                                                                                        <  1:300378/251‑1 (MQ=255)
                                                                                     ccAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGgtgt                                                                                                                                                         >  1:93578/1‑249 (MQ=255)
                                                                                       aGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGAc                                                                                                                                                     >  1:479206/1‑251 (MQ=255)
                                                                                       aGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGAc                                                                                                                                                     >  2:24666/1‑251 (MQ=255)
                                                                                            aaTCACCACATTATCGCCTTACGCTGCGAGCGCCTCGTTCAGCTTGTTGGGGATGATAGCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCAc                                                                                                                                                <  2:74851/251‑1 (MQ=255)
                                                                                               caccacATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCtt                                                                                                                                             >  1:300866/1‑251 (MQ=255)
                                                                                                     aTTATCGCCTTGTTCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATAACTCCTCAGAAGTGATACAGATCTTTCCCTCGGGCATTCTCAATACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGTAAAGCATTCGCAAAGGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGGCATCACCTTTAAACg                                                                                                                                       <  2:67982/251‑1 (MQ=255)
                                                                                                          cGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCGCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGCCGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGa                                                                                                                                  <  1:274453/251‑1 (MQ=255)
                                                                                                                       gAGCGCCTCGGTCAGCTTGTTGGTGATGATATCTCCCCATAAGTGAGACAGAGCTTGCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCTGACTATAAGGCTTCATGAAATCGAGCGGGCTCAGTACGCCATACAAGCCGGAAAGCATGCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAGGGTTTCGGCCTGCAAGACTTTTTAGACATCACCTGTAAACGCCTGAATCGCATGGCggg                                                                                                                     <  2:300758/251‑1 (MQ=255)
                                                                                                                                gTTCAGCTTGTTGGTGAGGATCTCTCCCCAGAATTGAAACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGcc                                                                                                             <  2:18648/250‑1 (MQ=255)
                                                                                                                                   cAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGc                                                                                                          <  1:24666/250‑1 (MQ=255)
                                                                                                                                     gCTTGTTGGTGATGATATCCCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGCCGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTg                                                                                                       <  2:243869/251‑1 (MQ=255)
                                                                                                                                                 tgatATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCCCCGGGGTGAAATCTGGCTGc                                                                                           >  2:409484/1‑251 (MQ=255)
                                                                                                                                                                aTTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAg                                                                                        >  1:323797/1‑239 (MQ=255)
                                                                                                                                                                aTTGATACAGATCTTCCCCTCGGGCATTCCCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAg                                                                                        <  2:323797/239‑1 (MQ=255)
                                                                                                                                                                  tGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAagcgagc                                                                          >  1:58918/1‑251 (MQ=255)
                                                                                                                                                                    aTACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAACGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAgcgg                                                                        >  1:374535/1‑251 (MQ=255)
                                                                                                                                                                      aCAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACTCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGATTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAAAGCCAGAATCGCCTGGAGGGCATTCTCCGACATGAAATCTGGGTGACAGTCATGaaggagacctgcg                                                                      >  2:373642/1‑240 (MQ=255)
                                                                                                                                                                       cAGATCTTTCCATCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGGGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGt                                                                     <  1:373642/251‑1 (MQ=255)
                                                                                                                                                                           tCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAACCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGa                                                                  <  2:8798/250‑1 (MQ=255)
                                                                                                                                                                              ttCCCTCGGGCATTCTCAAGACGGATCCCAATATCCAGACGAATAGGCTGCATTAAAGCGAGGGGTCGCAGTACGCCAGACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTGTCGGCCTGCAAGCTGGTGTAGACATCACCTTTAAACTCCAGAATCGCCTGGCGGGCAGTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATAc                                                               <  2:482593/250‑1 (MQ=255)
                                                                                                                                                                              ttCCCCCGGGCTTTCTCAAGACGGATCCCCTTTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATa                                                                <  2:140159/249‑1 (MQ=255)
                                                                                                                                                                              ttCCCCCGGGCATTCCCAAGGCGGACCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAACGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATa                                                                <  2:374535/249‑1 (MQ=255)
                                                                                                                                                                               tCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATAccc                                                             >  1:170476/1‑251 (MQ=255)
                                                                                                                                                                                 ccTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGACGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTTGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGc                                                           >  1:241193/1‑251 (MQ=255)
                                                                                                                                                                                  cTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGcc                                                          >  1:457179/1‑251 (MQ=255)
                                                                                                                                                                                   ccGGGCATCCTCAAGCCGGTTCCCAATTTCCAGGCGAAAAGGCGGCATTAAATCGGGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTGGGGCAAAATCGAAACCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAGTCTCCGGCGTGAAATCTGGCGGCCAGTCAGGAAAGCGAGCGGCGTTGATACCCGCCa                                                         <  2:479206/250‑1 (MQ=255)
                                                                                                                                                                                        cATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTg                                                    >  1:376440/1‑251 (MQ=255)
                                                                                                                                                                                         aTTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTATATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGAAACCCGCCCGTTTGt                                                   >  2:326153/1‑251 (MQ=255)
                                                                                                                                                                                           tctcAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAACTCGAAATCGTCTTCGCTGAAGGTTTCGGCATGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCACCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGACCGGCGTTGATACCCGCCAGTTTGTCg                                                 >  1:227567/1‑251 (MQ=255)
                                                                                                                                                                                                 gACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCCGTTTGTCGCTGATg                                           >  1:167076/1‑251 (MQ=255)
                                                                                                                                                                                                           cATTCCCAGACGATAAGTCTGCATTAAACCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGTTGTTGGGCAAAATCGAAATCTTCTCCGCTGAAGTTTTCGGCCTGCAAGCCGGTGTAGACATCCCCTTTAAACGCCAGAAGCGCCTGGCGGGCATTCTCCGGCGTAAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCg                                 <  1:353416/251‑1 (MQ=255)
                                                                                                                                                                                                                  ataCGATAAGGCTCAATTAAAGCGATCGGGCGCATTCCGCCTTACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCAAAATCGTCTTCGCTGAAGTTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAACCTGGCTGCCAGTCATGAAAGCGAGCGGGGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAAt                          <  1:119191/249‑1 (MQ=255)
                                                                                                                                                                                                                      gATAAGGCTGAATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTCAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCACTGGTGCGCATCAGAGTGCTAATCTGc                      >  2:308388/1‑251 (MQ=255)
                                                                                                                                                                                                                             cTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCg               >  1:383187/1‑251 (MQ=255)
                                                                                                                                                                                                                                cATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCCTGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAggagtac            >  2:462138/1‑246 (MQ=255)
                                                                                                                                                                                                                                     aaTCGAGCGGGCGCGGTACGCCATACAAGCCGGAAAGCATTAGCAAATGCTGTTGGGCAAAATCGACATAGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGGGTAGACCTCACATTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCAGCGTTGATACCAGCCAGTTTGTCGCTGGTGCGCATCAGCGTGCTAATcggcgggggagtaaggttc       >  2:218812/1‑238 (MQ=255)
                                                                                                                                                                                                                                     aaTCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTc       >  1:28122/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aGCGGGCGCAGTACGCCATACAAGCCGGCAAGCATTCGCAACTGCTGTTGGGCAAAATCGCAATCGTCTTCGCTGCCGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGACAGCGACCGGCGTTGATACCCGCCAGTTTGTCGCGGATGCGCATCAGCGTGCTCATCTGCGGAGGCGTCAGATTgcgcgc  >  2:21704/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
CCGTTCTTCTCATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCACATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGC  >  NC_000913/5864‑6348

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