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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009136,3490CT96.8% 176.0 / ‑2.5 63E37E (GAG→GAA) yaaADNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (30/31);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

GCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAATGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCG  >  NC_000913/6112‑6584
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            
cctcaacccccctcttccattctgCAAATGTTGTTGGGCAAAACCAAAACCGTCTTCGCTGAAGGTTTCGGC‑TGCCAAGCCGTTGTAGACATCACCTTAAACCGCCAGAATCGCCTGGCGGGAATTCCCCGGCGTAAAATTTGGCT‑GCCAGTCATAAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAAct                                                                                                                                                                                                                                 <  2:386061/228‑1 (MQ=255)
  cATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGAAAAATCGAAATCGTCTCCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAActgc                                                                                                                                                                                                                               <  1:659388/251‑1 (MQ=255)
   aTACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAActgct                                                                                                                                                                                                                              >  1:153083/1‑251 (MQ=255)
           cGGAAAGCATTCGCAAATGATGTTGGGCAAAATCGAAATCGTCTTAGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCCGTCATGAAGGAGAACGGCGTTGATACCCGCCAATTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAAGTTCCGCGCTTCATGGATAAACGTCTGGGAATTg                                                                                                                                                                                                                      >  2:676507/1‑251 (MQ=255)
              aaaGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTGTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCt                                                                                                                                                                                                                   <  1:166533/251‑1 (MQ=255)
               aaGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCt                                                                                                                                                                                                                   >  1:615531/1‑250 (MQ=255)
                aGCATTCGCAACTGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCAGTGTAGACATCACCTTTAAACGCAAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGGTTACGCGCTTCATGGATCAACTGCTGGGAATTGTCTa                                                                                                                                                                                                                  >  1:372901/1‑250 (MQ=255)
                   aTTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAATCGCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAAc                                                                                                                                                                                                                >  1:631213/1‑249 (MQ=255)
                      agCAAATGCTGTTAGGCAAAATGGAAATAGTCTTCGCTGAAGGTCTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGACAAAATCGCATGGCGGGCATTCTACGGCGTGACATATGGCT‑GCAAGTCATGAAAGAGAGGGACGTTGATACCCACAATATTGACGCTGATGCGGAACAACGTGCATATATGAGGAGGAGTCGGGTTACTCGCTTCATGGCACAACGGCTGGGAATGTTCTAACAGCTc                                                                                                                                                                                                           >  2:652653/2‑251 (MQ=255)
                       gCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAAct                                                                                                                                                                                                                                 >  1:386061/1‑228 (MQ=255)
                       gCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTcc                                                                                                                                                                                                          >  1:396470/1‑251 (MQ=255)
                            tGCTGTTGGGCAAAATCGAACTCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGACAGCGAGCGGCGTTGATACCCGCCAGTTTTTCGCTGATGCGCATCCGAGTGCTACACTGCGGAGGCGTAAGTTTCAGCGCTTCCTGGATCAACTGATGGGAATTGTTTACCAGCTCCGGccg                                                                                                                                                                                                     >  2:634000/1‑249 (MQ=255)
                                      cAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAgcgc                                                                                                                                                                                           >  1:466424/1‑251 (MQ=255)
                                       aaaaTCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGcgcg                                                                                                                                                                                          >  1:494217/1‑251 (MQ=255)
                                           ccgaaatcgTCTCAGCTGAAGGTTCCGGCCTGA‑AAGCCGGTGTAGACATCACTTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCgtggt                                                                                                                                                                                      <  1:664721/250‑1 (MQ=255)
                                              aaatcgTCGTCACTGAAGGTTTCGGACTGC‑AAGCCGGTGTAGACATCACCTTTAACCGCCAGAATCGCCTGGCGGGGATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGGGCGGCGTTGATACCCGCCAGTTTGGCGCTGATGCGCATCAGCGTGCTAATATGCGGAGGCGTGAGTTTCCGAGCTTCAAGGATCAAATGAGGGGCATTATCTGACAGCTCGGGCAGCGTAAAGCGCGtggtggg                                                                                                                                                                                   >  2:333223/1‑250 (MQ=255)
                                                       tCGCGGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGCCCTCACCTTTAAACGCCAGAATCGCCTGGCGGCCATTCTCCGGCGTGAAATCGGGCT‑GCAAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACgg                                                                                                                                                                            <  1:373925/248‑1 (MQ=255)
                                                         gCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCttt                                                                                                                                                                       >  1:434204/1‑251 (MQ=255)
                                                             atggTTGCGGCCTGC‑AAGCCGTCGGAGAAACCCCCTTAAAACGCCAGAATTGCGTGCCGGGCACTCGCCGGCGATAAGTATGGCT‑GCATGTCATGAAAGCGAGTGTGGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGGGGTGG‑TCAACGTGCTTTGGTa                                                                                                                                                                   <  2:631213/249‑1 (MQ=255)
                                                                 tttCGGCCTGC‑AGGCCGGTGGAGACATCACCTTTAAACGCAAGAACCGCCTGGCGGGCATGCTCCGGCGTGAAATCGGGCT‑GCCGGTCATGAAAGCGAGGGGGGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATGGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATc                                                                                                                                                                <  1:305372/250‑1 (MQ=255)
                                                                     cggctggaaaGCCGGTGTAGACATCACCTTAAAACGCCATAAGCCTCTGGCGGGCATTCTCCGGTGTGAAATCTGGCT‑GCCCGTCATGAAAGCGGGCGGGGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGc                                                                                                                                                            <  1:126040/243‑1 (MQ=255)
                                                                                ccGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAgg                                                                                                                                                 >  1:109737/1‑251 (MQ=255)
                                                                                                                  ccgggCGGGCATTCTCTGGCGGTAAATCTGGCT‑GCCAGTCAGGAAAGCGAGCGGTGTTGATACCCGCCAGTTTGGCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAgg                                                                                                               <  1:634000/248‑1 (MQ=255)
                                                                                                                     ggCGGGCATTCGCCGGCGTGAAGTC‑GGCTCGCCAGTCATGAAAGCGGGGGGCGTTTACACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCCGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGa                                                                                                              <  1:263979/249‑1 (MQ=255)
                                                                                                                      gggtgctttccccGGCGTGAAATCTGGCT‑GCCAGCCAGAAAACCGAGCGGGTTTAATCCCCGCCAGGTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGaaaa                                                                                                           <  2:294486/240‑2 (MQ=255)
                                                                                                                        gggCATTCCCCGGCGGAAAATCTGGCT‑GCCAGTTATAAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAAttt                                                                                                         <  1:518608/251‑1 (MQ=255)
                                                                                                                             ttctccGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGT‑TCAACGGGATTTGGTAATCAAGCGATTTCGCAGGTGAATTAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGcc                                                                                                    >  2:652747/5‑251 (MQ=255)
                                                                                                                                      gggaatattgggtggcaggTCAGGAAGGCGAGAGGTGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAg                                                                                            <  2:153083/233‑1 (MQ=255)
                                                                                                                                      gTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGaa                                                                                                                                      >  1:298576/1‑208 (MQ=255)
                                                                                                                                      gTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGc                                                                                           >  1:452925/1‑251 (MQ=255)
                                                                                                                                       tGAAATCAGGCT‑GCCAGACATGAAAGCGAGCGGCGTTGATACCCGCCAGATTGTCGCTGCTGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGATGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGGTTTCGCAGGTGAAATAAGAATCAGCCTATCCAGTCCTTGCAGGAAATTAATGCGGACTTTAGc                                                                                           >  2:585904/1‑250 (MQ=255)
                                                                                                                                          aatgctggtggccagctaTGAAAGCGAGGGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGGCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCaaa                                                                                        <  1:674339/234‑1 (MQ=255)
                                                                                                                                                  t‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCAGATGCGCATAAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑GCAACGGGCCTTGGTAATCAAGCGTTTGAGCAGGTGAAATAAGAATCAGCATATCCAGACCTTGCAGGAAATTTATGCCGACTTTAGCAGAAaaagaga                                                                                >  2:297440/1‑250 (MQ=255)
                                                                                                                                                          aaTGAAAGCGAGCGGCGTGGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGtt                                                                         <  1:469112/249‑1 (MQ=255)
                                                                                                                                                              aaaGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGAGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCg                                                                    >  1:334178/1‑251 (MQ=255)
                                                                                                                                                                agcgagcgGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGAt                                                                  <  1:144626/251‑1 (MQ=255)
                                                                                                                                                                   gagcgGCGTTGATAAACGCCAGTTTGTAGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCAAGTCCTTGCAGGAAATTTATGCCGCCTTTAGCAAAAAAAGAGAATGAGTTGATCGaacgt                                                               >  2:609270/1‑247 (MQ=255)
                                                                                                                                                                     gcgggttttatccccGCCAGGTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTTCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTtg                                                             <  2:434204/239‑1 (MQ=255)
                                                                                                                                                                       ggcgTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTtgtg                                                           <  1:454508/251‑1 (MQ=255)
                                                                                                                                                                       ggcgTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTtgtg                                                           >  1:395495/1‑251 (MQ=255)
                                                                                                                                                                         ggttgacccccGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGAt                                                         <  1:609999/243‑1 (MQ=255)
                                                                                                                                                                            tGATCCCCGCCAGGTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTGCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTa                                                       <  1:465499/250‑1 (MQ=255)
                                                                                                                                                                                   cGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGTAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACgg                                                                                                                                                                            >  2:15309/1‑126 (MQ=255)
                                                                                                                                                                                   cGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACgg                                                                                                                                                                            <  1:15309/126‑1 (MQ=255)
                                                                                                                                                                                   cGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCg                                               <  1:385562/251‑1 (MQ=255)
                                                                                                                                                                                        gTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAAc                                           >  1:569042/1‑250 (MQ=255)
                                                                                                                                                                                         tttGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCTTATAGCGCGTGGTGG‑TCAACGGGCTTTGTTAATCACGCGTTTTCGAAGGTGAAATAAGAATCAGCATGTCCAGTCCTTGCAGGAAATTTATGCCGACTTTAACAAAACAAGAGAATGAGTTGATCGATAATTGTGATTACTCCTGCGAAAca                                          >  2:26317/1‑250 (MQ=255)
                                                                                                                                                                                          gtgttgcttaTGCTCATCAGCGTGGTAATCTGCGGAGGCGTCAGTCTCCGCGCTTTATGGAACAACGGGTGTGAATTGTCGAACAGCTCCGGCAGCGTATAGCGCG‑GGTGGGTCAACGGGCTTTGGTAATCAAGCGTTTGCGCATGTTAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGGCTTTAGCAAAAAAAGAGAATGAGTTTTTCGATAGTTGTTATTACTCCTgttgaacat                                         <  2:620178/241‑9 (MQ=255)
                                                                                                                                                                                            gTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAAcatcat                                      >  2:214282/1‑251 (MQ=255)
                                                                                                                                                                                            gTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAAcatcat                                      <  1:345334/251‑1 (MQ=255)
                                                                                                                                                                                               gCTGAGGCGCATCAGCGTGCTAAGCTGTGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTAGAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATTAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATccc                                   <  2:580282/251‑1 (MQ=255)
                                                                                                                                                                                                   aTGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATACACAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCAcgc                               >  1:583284/1‑251 (MQ=255)
                                                                                                                                                                                                        cATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGTGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCGTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCg                          <  1:609270/251‑1 (MQ=255)
                                                                                                                                                                                                         aTCAGCGTGCCAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCgg                         >  1:109994/1‑251 (MQ=255)
                                                                                                                                                                                                           cAGCGTGCTAATCATCGGAGTCGTCAGTTTCCGCGTTTCAGGGATCAACTGCTGGGATTTGTCTAACAGCTCCGGCCTCTTATAGCGCGTGGTGG‑GCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATTACCAGTCCTTGCATGAAATTTATGCCGACTTTTGCAAAAAAAGATAATTAGTTGCTCGATAGTTGTTATTATTTCGGTGAAACATCATCCCACGCGTCCGgag                       <  2:445188/251‑1 (MQ=255)
                                                                                                                                                                                                                ggCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCAAGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGc                  <  1:26317/250‑1 (MQ=255)
                                                                                                                                                                                                                            gAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTACTCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGat       >  1:281687/1‑250 (MQ=255)
                                                                                                                                                                                                                             aGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGTGCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATAGCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCTAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGatat     <  2:494217/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATc    <  1:652747/251‑1 (MQ=255)
                                                                                                                                                                                                                               gCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGCTAGTTGTGATTACTCCTGCGAAACCTCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATcc   >  2:596861/1‑251 (MQ=255)
                                                                                                                                                                                                                               gCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGAGCGATAGTTGTGATTACTCCTTCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATcc   <  2:566070/251‑1 (MQ=255)
                                                                                                                                                                                                                                cGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCg  >  1:80741/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGc                                                                                           <  1:139123/160‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                            
GCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGC‑AAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCT‑GCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGG‑TCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAATGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCG  >  NC_000913/6112‑6584

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