Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,375,737 G→A A259T (GCC→ACC)  ycjQ → D‑guloside 3‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,375,7370GA98.3% 164.0 / NA 60A259T (GCC→ACC) ycjQD‑guloside 3‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (22/35);  minor base C (1/0);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CGCTGGCGCTTCTGTGGTGATTGGCGTCGATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTGCCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGGTGGTGACCTTTGCCAATAGCCCGGAAAGCTATATGCAGTATGTCGACCAGCATCCG  >  NC_000913/1375499‑1375982
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
cgcTGGCGCTTCTGTGGTGATTGGCGTCGATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTg                                                                                                                                                                                                                                            <  1:188609/251‑1 (MQ=255)
       gcTTCTGTGGTGATTGGCGTCGATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAgg                                                                                                                                                                                                                                    <  1:675295/251‑1 (MQ=255)
          tCTGTGGTGATTGGCGTCGATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGtt                                                                                                                                                                                                                                  <  1:244094/250‑1 (MQ=255)
             gtggttATTGGCGTCGATCCTATCGCCCAGCGTTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCATCATCTCCTATGTCGCGTTTACCAAGCCGGTTGCCGAAGGttt                                                                                                                                                                                                                                 <  2:635482/248‑1 (MQ=255)
                gtgATTGGCGTCGATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGGTTGCCGCAGGGTTTAACCt                                                                                                                                                                                                                           >  1:55852/1‑251 (MQ=255)
                 tgATTGGCGTCGATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAGGCAGGGGGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTGGCGTTTACCAAGCCGTTTGCCGATGGTTTTAACCTc                                                                                                                                                                                                                          >  1:151486/1‑251 (MQ=255)
                  gATTGGCGTCGATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTGTACCAAGCCGTTTGCCGAAGGTTTTAACCTc                                                                                                                                                                                                                          <  1:186992/250‑1 (MQ=255)
                   aTTGGCGTCGATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATATCCTATGTTGCGTTTACCAAGCCGTTTGGCGCAGGTTTTAAACTCgg                                                                                                                                                                                                                        >  1:386745/1‑251 (MQ=255)
                   aTTGGAGTCGATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGCTGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGTTGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGGCTGGCTTAAGGCGGCACCCACTCCCATGTGGCCGTTTACAAGCCGTTTGCCGGAAGGtttgacccccg                                                                                                                                                                                                                         >  2:554011/1‑243 (MQ=255)
                         gTCGATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGCTTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGAcgcg                                                                                                                                                                                                                   >  1:283955/1‑250 (MQ=255)
                            gATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGAc                                                                                                                                                                                                                      <  1:168509/244‑1 (MQ=255)
                            gATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTAGCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCGTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAgc                                                                                                                                                                                                               <  1:84414/251‑1 (MQ=255)
                                   ccgcgcatccctttttgatTGCCCGTCGCTACGGCGCGGATTTTTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGTGATCAAAACGCTGACCGTCAAGCAGGGTGCCGCTGTGATTATCGAAACCATCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTAGGGCGGCACCATCTCCTATGTCGCGGTTACCAAGACGGTTGCCGAAGGTTTTAACCTCGGTCGCGAAGCGCATTTc                                                                                                                                                                                                        <  2:458724/234‑1 (MQ=255)
                                           cgccgttatatttccccacGCCACTGGGCGGATTACTGCCTTTCCCCCATCGGCACTGTGGATGGTAAAGAGATCAAAACGCGGACCGGCACGCAGTGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGGCTGGCTGACGGCGGCACCATCCCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGTGAAGCGCATTTCAATAACGc                                                                                                                                                                                                <  2:406298/233‑1 (MQ=255)
                                           cGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGCCCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCCTCTCCTATATAGCGATTACCAAGCCGTTTGCCGGAGGTTTTAACCTCTGACGGGAAGCGCATTTCAATAAcac                                                                                                                                                                                                >  2:447022/1‑249 (MQ=255)
                                                      ttcccgtggctactgcgcGGATTGCTGCCTTAACCCCAGCGGCCCTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAAGAACGCCAAAATTGtc                                                                                                                                                                                      <  2:86716/236‑1 (MQ=255)
                                                      tGCCCGTCGCCACGCCGCGGATTTCTGCCTCAACCCCATCGGCTCTGATGTAGGTAAAGAGATCAAACCGCTGACCGGCAAGCAGGGTGCCGATGTGATTCTCCAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCGCCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGtc                                                                                                                                                                                      <  1:452246/250‑1 (MQ=255)
                                                       gCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGAGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGGCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGtctt                                                                                                                                                                                    <  2:120167/251‑1 (MQ=255)
                                                              gCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGAAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCCTCTCCTATGTCGCGTTTAACAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTtctctagc                                                                                                                                                                             >  1:372/1‑248 (MQ=255)
                                                                        ggATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGTCTTTCAGGGTGCCGATGTGATTAGCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGATGGTTTTTTCCTCgg                                                                                                                                                                                                                        <  2:421280/198‑1 (MQ=255)
                                                                        ggATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGcg                                                                                                                                                                   <  1:480289/251‑1 (MQ=255)
                                                                        ggATTTCTGCCTTAACCCCATCGGCACTGAGGAAGGGAAAGAGAGAAAAAAGCTGAACGGCAAGAAGGGAGACGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCtagg                                                                                                                                                                                                                        >  1:421280/1‑195 (MQ=255)
                                                                            ttCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGcgaac                                                                                                                                                                <  1:619793/250‑1 (MQ=255)
                                                                                gCCTTAACCCCATCGGCACTGATGGAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGGTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGcgaaccg                                                                                                                                                              <  1:44083/248‑1 (MQ=255)
                                                                                           aTCGGCACTGATGTAGGTAAAGATATCAAAACGCTGTCCGGCGAGACGGGTGCGGATGTGATTATCGGTCCCAGCGGCTACGCCGACGCTCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTATCAAGCCGTCTGCCGAAGGTTTTAACCTCGGACGCGCAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCATCGAACCGAACCCGGATTTTcc                                                                                                                                                <  2:103878/251‑1 (MQ=255)
                                                                                             gggCACTGATGGAGGTAAAGAGATCAAAACGCTTACCGGCAAGCAGGGTGCCGATGTGAATATCGAAACCAGCGGCTACGCCGACGCGCTGCAAGCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGATTGCCGAAGGTTTTAACCGCGGACGCGAAGCGCATTTCAATAACGCCAAAAGTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCcgc                                                                                                                                              <  2:642221/250‑1 (MQ=255)
                                                                                              ggCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCcgcg                                                                                                                                             <  1:274006/251‑1 (MQ=255)
                                                                                              ggCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACCCCGACGCGCTGCAATCGGACCCCCGCGGTCTGGCTTATGGCGGCAACTTCTCCTATGTCGCATTTACCACGCCGTTTGCCGAAGGTTTTAACCCCGGACGCGAAGCGCATTTCCATAACGCCAAAATTGTCTTCCCTCGCGCGAGCAGCGAACCGCCCCCGGATTATCcgcc                                                                                                                                             >  1:642221/1‑250 (MQ=255)
                                                                                               gCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCgcgc                                                                                                                                            >  1:369507/1‑251 (MQ=255)
                                                                                                   tGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGttt                                                                                                                                                                                                                                             >  1:191301/1‑151 (MQ=255)
                                                                                                   tGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGttt                                                                                                                                                                                                                                             <  2:191301/151‑1 (MQ=255)
                                                                                                    gATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTgg                                                                                                                                         <  1:401387/249‑1 (MQ=255)
                                                                                                            gTAAAGAGATCAAAACGCTGACAGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCCACCCCTATTTCGCGTTTCACCAGCCGTTTGCCGAGGGGTTTAACCTCGGACGCGAAGCTCATTTCAATAACGCCAAAATTGTCTTCTCTCCCGCGGGCAGCGCAACGCACCCGGATTAtaccagctggtgccgcc                                                                                                                                  >  2:495474/1‑233 (MQ=255)
                                                                                                              aaaGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCCAGCg                                                                                                                              >  1:124700/1‑251 (MQ=255)
                                                                                                               aaGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGt                                                                                                                             <  1:562775/251‑1 (MQ=255)
                                                                                                                  agaTCAAAACGCTGACCGGCAAGCAGGGTGCAGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATATCCCGTTTACCAAGCCGGTTGACGAAGGTTTTAACCTCGGATGCGAAGCGCATTTCAATAAAGCCATAATTGTTTTCTCTTGCGCGTGCGACGAACCGAAACCGGCATATCCGCGCTGGAGCCGCAAgggattt                                                                                                                          >  2:65957/1‑245 (MQ=255)
                                                                                                                                aCCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGATGCAATCGGCGCTCCGCGATCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTAGGTCGCGAAACAAATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGCACCGAACCCGGATTATCCGCGCTGGAGCCGCAAACATATTGAAGAAACCTGTTg                                                                                                            >  1:386977/1‑251 (MQ=255)
                                                                                                                                    gCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTCTGTCGCGTTTACCAAGCCGTTTGCCGAAGGATTTAACCTCGGCCGCGAAGCGCCTTTCAATAAAGCAAAAATTGTCTTCTCTAGCGCGTGCAGCGAACCGAACCCGGCTTATCCGCCGTGGATCCGCAAGAGTGTTGCAGAACCCTGTTGGGa                                                                                                         >  2:61511/1‑250 (MQ=255)
                                                                                                                                    gCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAg                                                                                                                                >  1:325067/1‑227 (MQ=255)
                                                                                                                                    gCAAGCAGGGTGCCGATGTGATTAGCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGTTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAg                                                                                                                                <  2:325067/227‑1 (MQ=255)
                                                                                                                                     cAAGCAGGGTGCCGATGGGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCGCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCAGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAAc                                                                                                       <  2:282560/251‑1 (MQ=255)
                                                                                                                                          aGGGTGCCGATGTGATTATCGAAAACAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTAGGGCGGCACCATCTCCTATGTCGCGTTTACCAAACCGTTTGCCGAAGGTATTAACCTCGGACGCGAAGCGCATTACACTAACGCCAAAATTGTCTTCTTTCGGGCGTGCAGCGAAACGACCACGGATTATCCGCGCTGGAGACGACAGCGTATTgacgaacccaggttggaaaggctg                                                                                                  >  2:606849/1‑233 (MQ=255)
                                                                                                                                             gcgccgttgcgATTATCGAAACCAGCGGCTACGGCGACGCGTTGCAATCGGCGCTCCGCGGTCTGGCTTAGGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCtgatg                                                                                               <  2:151486/241‑1 (MQ=255)
                                                                                                                                                    tgtgATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGtt                                                                                        >  1:438701/1‑251 (MQ=255)
                                                                                                                                                      tgATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGCAGGTTTTAACCTCGGACGCGAAGAGCATTTCAATAAAGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAAACCGGATTATCCGCGCTTGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGGACTGCTGATGAACGGTTAt                                                                                      >  1:399433/1‑251 (MQ=255)
                                                                                                                                                       gATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCGAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAAc                                                                                                                                                                                                  <  1:458592/142‑1 (MQ=255)
                                                                                                                                                       gATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCGAGCCGTTTGACGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAAc                                                                                                                                                                                                  >  2:458592/1‑142 (MQ=255)
                                                                                                                                                              gAAACCAGCGGCTACGCCGACGCGCTGCAATCGGGGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGGCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCGCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATTTCAATTg                                                                              <  2:283315/251‑1 (MQ=255)
                                                                                                                                                                       ggCTACGCCGACGCGCTGCAATCGGAGCTCCGAGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGAGCATTTCAATAAAGGCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGATGGAGACGCGAGCGTATTGAAGAAACATGTTGGGAACGGCTGCTGAACGGGTATATAAATTGAGAAGAttt                                                                     >  1:662686/1‑251 (MQ=255)
                                                                                                                                                                           aCGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAAt                                                                  <  1:447022/250‑1 (MQ=255)
                                                                                                                                                                                gtctctcTGCAGTCGGCACTCCGCGGTCTGGCTTAGGTCGGCACCATGCCGCATTTCGCGTTTACCAAGCCGTTTGCCGAAGGTGTTATCCTCGGCCGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGTGAACCGTACCTGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAATCTGTTGGGATCTGCTGATGAACGGATATCTCCATTGCGAAGTGTTAATCGAccc                                                            <  2:79551/245‑1 (MQ=255)
                                                                                                                                                                                        gcaattggtgcgccGCGGTCTGGCTTATGGCGGCACCATCTCCTATTTCGCGTTTATCAAACCGTTTTCCGAAGGTTTTACCCTCGGACGTGAATCGCATTTCAATAACGCCAAAATGGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCGGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCGGATGAACGGTTAGCTCAATTGCGAAGATTTAATCGTCCCggtggt                                                      <  2:542775/237‑1 (MQ=255)
                                                                                                                                                                                        gCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTGCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGgtggtg                                                     <  1:554011/250‑1 (MQ=255)
                                                                                                                                                                                                  ccccGCGGTCTGGCTTCTGGCTGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTGCTCTCGCTCGTGCAGCGTACCGATCCCGGATTATCCGCTCTGGAGCCGCAAGCGTAATGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGGTGGTGACCTTTGcc                                            <  2:386745/247‑1 (MQ=255)
                                                                                                                                                                                                  atcctcgGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTACCAAGCCGTTTGCCGGAGGTTTTAACCCCGGACGCGAAGCGCATTTCAATAACACCCCAATTGTCTTCTCTCGCCCGTGCAGCGGACCGAACCCGGCTTATCCGCTCTGGAGCCGCAAGAGTATTGAAGAAACCTGTTGGGCACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGAcccaggggggacccttgcga                                           >  2:477092/6‑233 (MQ=255)
                                                                                                                                                                                                         gTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTGCCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGGATTGAAGAAACCTGTTGGGACCTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGGTGGTGACCTTTGCCAATAGCCCg                                   >  1:225196/1‑251 (MQ=255)
                                                                                                                                                                                                                                 ccTATGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGGTGGTGACCTTTGCCAATAGCCCGGAAAGCTATATGCAGTATGTCGa            <  1:606849/250‑1 (MQ=255)
                                                                                                                                                                                                                                    aTGTCGCGTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGGTGGTGACCTTTGCCAATAGCCCGGAAAGCTATATGCAGTATGTCGACCAg        <  1:495474/251‑1 (MQ=255)
                                                                                                                                                                                                                                        cgcgTTTACCAAGCCTTTTGCCGAAGGGTTTAACCTCTGACGCGAAGCGAATTGCAATAACGCCAAAATTGGCTTCTCCCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACGGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGGTGGTGACCTTTGCCAATAGCCCGGAAAGCTATATGCAGTATGCCGACCAgtat     <  2:3332/250‑4 (MQ=255)
                                                                                                                                                                                                                                          cgTTTACCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCACTAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGCAGATTTAATCGACCCGGTGGTGACCTTTGCCAATAGCCCGGAAAGCTATATGCAGTATGTCGACCAGCATCCg  >  2:92141/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
CGCTGGCGCTTCTGTGGTGATTGGCGTCGATCCTATCGCCCATCGCTGTGATATTGCCCGTCGCCACGGCGCGGATTTCTGCCTTAACCCCATCGGCACTGATGTAGGTAAAGAGATCAAAACGCTGACCGGCAAGCAGGGTGCCGATGTGATTATCGAAACCAGCGGCTACGCCGACGCGCTGCAATCGGCGCTCCGCGGTCTGGCTTATGGCGGCACCATCTCCTATGTCGCGTTTGCCAAGCCGTTTGCCGAAGGTTTTAACCTCGGACGCGAAGCGCATTTCAATAACGCCAAAATTGTCTTCTCTCGCGCGTGCAGCGAACCGAACCCGGATTATCCGCGCTGGAGCCGCAAGCGTATTGAAGAAACCTGTTGGGAACTGCTGATGAACGGTTATCTCAATTGCGAAGATTTAATCGACCCGGTGGTGACCTTTGCCAATAGCCCGGAAAGCTATATGCAGTATGTCGACCAGCATCCG  >  NC_000913/1375499‑1375982

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