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→CTG) yaaADNA 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

CATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCCT  >  NC_000913/5874‑6350
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       
cATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCg                                                                                                                                                                                                                                       >  1:549015/1‑251 (MQ=255)
       gAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTAc                                                                                                                                                                                                                                                    >  2:150882/1‑231 (MQ=255)
       gAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTAc                                                                                                                                                                                                                                                    <  1:150882/231‑1 (MQ=255)
          cacCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑c                                                                                                                                                                                                                            <  1:46480/251‑1 (MQ=255)
                 ttGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATg                                                                                                                                                                                                                     >  1:489368/1‑251 (MQ=255)
                      aaTCTCGGGATTCAATTGCTTCGGCTTACCCGATTAAAAATACCAATCTGACGCCAGATTAATCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCGTGTTGGTGATGATATAGCCCCAGAATAGATACAGATCTTTCCCGCGGGCAGTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAAACGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGtt                                                                                                                                                                                                                <  2:350470/251‑1 (MQ=255)
                           cGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCa                                                                                                                                                                                                           <  2:286857/251‑1 (MQ=255)
                              cATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGCAGCCAGATTAATCACCA‑CATTATAGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGCGCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGATGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGATAGCAGT‑CGCAAATGATGTTGGGCaaa                                                                                                                                                                                                         >  2:469112/1‑250 (MQ=255)
                                ttCATTTTCTTCTGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAACCACCA‑CATTATCGCCTTGTGCTTCGATCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATC‑CGCAAATGCTGTTGGGCAaaat                                                                                                                                                                                                       <  2:602221/250‑1 (MQ=255)
                                       tcttaGGCCTCACCGATTTCCAATACTCATCTGACGCCAGATTACTCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGCTCCCCATTTCCAGACGATCAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCAGGCACGCATT‑CGCAAATGCTGTTGGGCAaaatcgaaaccg                                                                                                                                                                                               >  2:674339/1‑248 (MQ=255)
                                          tcGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTc                                                                                                                                                                                             <  1:659312/250‑1 (MQ=255)
                                          ccggtttCACCGATTTAAAACACCCACCTGACGCCAGATTAATCACCT‑CATTATCGCCTGGTGCTGCGCGCGCCTCGTTCAGCTTGTTTGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTTTTGTGCAAAATCGAAATCGTc                                                                                                                                                                                             <  2:320074/245‑1 (MQ=255)
                                                    cGATTTAAAATTCTCATCTGTCGCCAGATTAATCACCA‑CATTATCGCCTTGTGCTGCGATCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGTCTAAATCGAAATCGTTTGCGCTGAAgg                                                                                                                                                                                  <  2:483204/251‑1 (MQ=255)
                                                    cGATTTAAAATACTCATCTGACGCCAGATTAATCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCAGGAAAGCATT‑CGCCAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAgg                                                                                                                                                                                  >  1:620178/1‑251 (MQ=255)
                                                          aaaaTACTCATCTGACGCCAGATTAATCACCA‑CATTATCGCCTTTATCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGTCAAAATCGAAATAGTCTTCGCTGAAGGTTTCgg                                                                                                                                                                            >  2:139157/1‑251 (MQ=255)
                                                             aTACTCATCTGACGCCAGATTAATCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGATGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCAATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAAAAGTCTTCGCTGAAGGTTTAGGCCt                                                                                                                                                                         >  2:78756/1‑251 (MQ=255)
                                                               ccTCATCTTACGCCATATTAAT‑ACTATCATTATCGTCTTGTGCTGCGAGCGCCTCGTTCTGCTTGGTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTGCAAGACGAGAAGGCTGCATTAAATCGATCGGGCGCAGTACGCCATACAAGCCGGAAATCATT‑CGCAAATGCTGTTGGTCAAAATCGAAATCGTCGTCGCTGAAGGTTTCGGCCTGc                                                                                                                                                                       <  2:361146/250‑1 (MQ=255)
                                                                       gACGCCAGATTAATCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTCCAAGCCGgt                                                                                                                                                               <  1:139157/251‑1 (MQ=255)
                                                                              gATTAATCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGAc                                                                                                                                                         >  1:566070/1‑250 (MQ=255)
                                                                               aTTAATCACCC‑CATTATCGCCTTGTGCTGCGAGCGCCCCGTTCAGCTTGTTGGTTATGATATCTCCCCAGAATTGATACAGATCTTTCCCCCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACAt                                                                                                                                                       <  2:549015/251‑1 (MQ=255)
                                                                                  aaTCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCAc                                                                                                                                                    >  2:166533/1‑251 (MQ=255)
                                                                                  aaTCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTACAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTTTAGACATCAc                                                                                                                                                    >  1:250342/1‑251 (MQ=255)
                                                                                     cacca‑cATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTATTCGCTAAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCtt                                                                                                                                                 >  2:518608/1‑251 (MQ=255)
                                                                                      acca‑cATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGACTTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTATACATCACCttt                                                                                                                                                >  2:659388/1‑251 (MQ=255)
                                                                                        ca‑caTTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCACTTTCCAGACGATAAGGCTGCCTTAAATCGCGCGGGCGCAGTACGCCATACAAGCCGGAAAGCCTT‑CGCAAACGCTGTTGGGCACACTGGACATCGTCTTCCCTGCAGCTTTCGCCCTGCAAGCCGGTGAAGaactcaccatta                                                                                                                                               >  2:184134/1‑239 (MQ=255)
                                                                                           caTTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAgg                                                                                                                                                                                  >  2:355389/1‑213 (MQ=255)
                                                                                           caTTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAgg                                                                                                                                                                                  <  1:355389/213‑1 (MQ=255)
                                                                                            aTTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGCTCTTTACCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAACTCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAAcc                                                                                                                                           >  1:445188/1‑250 (MQ=255)
                                                                                                 cGCCTTTTGCCGCGAGCGCCTCTTTCATCTTGGTGGTGATGATATCGCCCCAGAATTGATTCAGATATTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCATACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAATTCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTATACGCCAGa                                                                                                                                      <  2:469103/251‑1 (MQ=255)
                                                                                                        ggctgctAGCGCCTCGTTCAGCTTTGTGGGGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCTATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCAGT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCt                                                                                                                               <  2:598658/250‑1 (MQ=255)
                                                                                                             cGAGGGCCTCGGTCAGATTGTTGGTGATGATATCTCCCCAGAATTGATACAAATCTGTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCgg                                                                                                                          <  2:32262/251‑1 (MQ=255)
                                                                                                                 cgcCTGGTTCAGCTTGTTGGTGATGATCTCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCTCTGAATGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCTGAATCGCCTGGCTGGCAt                                                                                                                      <  2:467893/251‑1 (MQ=255)
                                                                                                                 cgcCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGcct                                                                                                                      >  2:30001/1‑249 (MQ=255)
                                                                                                                 cgcCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGATGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTACACTCCAGCATCGCCTGGCGGTCAt                                                                                                                      >  2:144626/1‑251 (MQ=255)
                                                                                                                  gcCTCGTTCAGCTTGTTGGTCATGCTATCTCCCCAGAATTGATACAGATCTTTACCTCGGGCATTCTCAAGACGGCTCCCCATTTCCAGACGCTAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAACCGGAAAGCATT‑CCCAAATGCTGTTGGGCAAAAACGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCACGCCGGAGTAAGCATCACCTTTAAACGCCCGAATCGCCTGGCGGGcctt                                                                                                                     >  2:347417/1‑248 (MQ=255)
                                                                                                                     tCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAACCGCCAGAATCGCCTGGCGGGCATTctc                                                                                                                  >  2:454508/1‑249 (MQ=255)
                                                                                                                       gTTCAGTTTGTTGGTGATGATCTCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCTTTCACAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAAACGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTctcc                                                                                                                 <  2:319862/250‑4 (MQ=255)
                                                                                                                        ttCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGACGGAAAGCATT‑CGCACATGCTGTTGGGCAACATCGAAATCGTCTTCGCTTAAGGTTTCGGCCTGCACGACGGTGTAGCCATCCCCTTTCAACGCCAGAATCGCCTGGCGGGCAttatccgg                                                                                                               >  2:174764/1‑245 (MQ=255)
                                                                                                                                  tgGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGGGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCCTTCTCCGGCTTGGAATCt                                                                                                     >  2:609999/1‑251 (MQ=255)
                                                                                                                                    gTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTgg                                                                                                   >  1:580282/1‑251 (MQ=255)
                                                                                                                                      gctgatATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTTGCt                                                                                                 >  2:45848/3‑251 (MQ=255)
                                                                                                                                          atatCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCa                                                                                             <  1:78756/251‑1 (MQ=255)
                                                                                                                                             tctcCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTc                                                                                          >  1:443623/1‑251 (MQ=255)
                                                                                                                                             tctcCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTc                                                                                          >  1:83315/1‑251 (MQ=255)
                                                                                                                                             tctcCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGACGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCAGGCTGCCCGTc                                                                                          >  2:385562/1‑251 (MQ=255)
                                                                                                                                               cccccAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCAt                                                                                        <  1:497737/250‑1 (MQ=255)
                                                                                                                                                ccccAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCa                                                                                         <  1:30001/249‑1 (MQ=255)
                                                                                                                                                       aTTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAagcg                                                                                 >  1:463931/1‑250 (MQ=255)
                                                                                                                                                         tGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAagcgagc                                                                              >  1:145433/1‑251 (MQ=255)
                                                                                                                                                           aTACAGATCTTTCCCTCGGGCATTCGTAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCTGGAAAGCATT‑CGCAAATGCGGTTGGGCAAAATCGAAATCGTCTTCGCGGAAGGTTTCGGTCTGCAAGCCGGTGGAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCGGGCTGCCAGTCATGAAAGCGAgcgg                                                                            <  2:618537/251‑1 (MQ=255)
                                                                                                                                                             aCAGATCTTTCCCTCGGGCATTCTCAAGACGGCTCCCCATTTCCAGACGATATGGCTGCCTTAAATCGAGCGGGCGCAGTCCGCCATCCACGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGCAATCGTCTTCGATGAAGGTTTCGGCCTGCACGCCGGTGTAGCCATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCCGTAAGGAAAGCGAgcggcg                                                                          >  2:305372/1‑251 (MQ=255)
                                                                                                                                                                     ttCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACAt                                                                                                                                                       >  1:636339/1‑166 (MQ=255)
                                                                                                                                                                     ttCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTT‑g                                                                      <  1:49905/246‑1 (MQ=255)
                                                                                                                                                                     ttCCCCCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACAt                                                                                                                                                       <  2:636339/166‑1 (MQ=255)
                                                                                                                                                                      ccccTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGGCATGAAAGCGAGCGGCGTT‑GATAccc                                                                <  2:114642/250‑1 (MQ=255)
                                                                                                                                                                                             aTCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCAt                                         >  1:294486/1‑251 (MQ=255)
                                                                                                                                                                                              tCCCCATTTCCAGACGATAAGGATGCCTTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGAAAATGCTGTTGGGCACAATCGAAATCGTCTTCGATGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCCTGCAGTCTGGCTGCCAGTCCTGCAAGGGGGCGGCGTTTGA‑ACCCGCCAGtttttggcttatgagcat                                         >  2:436236/1‑235 (MQ=255)
                                                                                                                                                                                                  cattgcaaGACGATAAGGCGGCATAAAATCGAGCGGGCGCAGTACCCCATACAAGCCGGAAAGTATT‑CGCAAATGCGGTTGGGCAAAAGCGAAATCGTCTGCGCGTAAGTTTTCGGCCTGCAAGCCGGTGTATACATCACCTTTAAACGCCATAATCGCCTGGCGGGCATTCTCCGTCGTGAAATCTGGCGGCCAGTCAGGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCg                                    <  1:436236/244‑1 (MQ=255)
                                                                                                                                                                                                   aTTTCCAGACGATAAGAATGCATTAAATCGAGCGGGCGCAGTACTCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGAAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTACGTAGTGTAATCTGGCTGCCAGTCATTAAAGCGCGCGGCATT‑GATACCCGCCAGTTTGACGCTGATGCGCATccgatt                                   >  2:44535/1‑246 (MQ=255)
                                                                                                                                                                                                           gcaaTAAGGCAGAATTAAAGCGAGCGGGCGCAGTACCCCAACCAGGCCGGAA‑GCCTTGCGAAAATGCTGTTGGGAAAAATCGAAATCTTCTTCGCTGAGGGTTCCGGCCGGAAAGCCGGTGTAGACATCACGTTTAACCGCCAGAATCGCCTGGCGGGCATTCCCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGGGAGGGGCGTT‑AATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCt                          <  2:145433/248‑1 (MQ=255)
                                                                                                                                                                                                                       cATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCAACTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGAAGCGCATCAGCGTG‑CTCATCTGCGGCGGCGTc               >  1:99905/1‑250 (MQ=255)
                                                                                                                                                                                                                            aaTCGAACGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACACCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGGCATGAAAGAGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAACCTGCGGAGGCGTCAGTTTc         >  1:291768/1‑251 (MQ=255)
                                                                                                                                                                                                                                gagtgggcgcAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCTGCGttt                                                                       <  2:636121/183‑2 (MQ=255)
                                                                                                                                                                                                                                gAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTT‑g                                                                      >  1:636121/1‑187 (MQ=255)
                                                                                                                                                                                                                                aaGCGGACGCAGTACGCCATACAAAACGGAAAGCATT‑CGCAAATGCTGTTGGGCACAATCGAAATCGTCTTCGCTGAAGGTTTAGGCCTGCAAGCCGGTGTAGACATCACCTTTACACGCCAGAATCGCCTGGCGGGCCTTCTCCGGCGTGAAATCTGGCGGCACGTCCAGAAAGCTCGCGGATTT‑GATCCCCGCAAGTTTGTGTCTTGTGCG‑ATAAGCATGGCTAATCTGCGGAGGGGTCAGTTTccacg     >  2:501103/2‑248 (MQ=255)
                                                                                                                                                                                                                                   ggggcgcAGTCCGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCgcgctt  <  1:45848/250‑3 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                       
CATCGAGGAACACCGGCTTGATAATCTCGGCATTCAATTTCTTCGGCTTCACCGATTTAAAATACTCATCTGACGCCAGATTAATCACCA‑CATTATCGCCTTGTGCTGCGAGCGCCTCGTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATT‑CGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTT‑GATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTG‑CTAATCTGCGGAGGCGTCAGTTTCCGCGCCT  >  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: ‑