Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,792,209 G→A V313M (GTG→ATG)  gabD → succinate‑semialdehyde dehydrogenase (NADP(+)) GabD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,792,2090GA98.4% 229.4 / NA 62V313M (GTG→ATG) gabDsuccinate‑semialdehyde dehydrogenase (NADP(+)) GabD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/40);  minor base C (1/0);  total (22/40)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAGTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGC  >  NC_000913/2791967‑2792452
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
ccGGTTCGACCGAAATGGGCCGCCAGTAAATGGAACAGTGCGCGAAAGCCATCAAGAAGGGGCCGCGGGAGCTGGGCGGTACCGCGCGGTTTATCGTCTTTGAGGAGGCCGCCCCCGACAAAGCCGGGAAAGGCGGGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCGGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCaaa                                                                                                                                                                                                                                             <  2:162740/251‑1 (MQ=255)
  ggtacGACCGAAATTGGCCGCCAGTTAATGGACCAGTGCGCGAAAGACATAAAGAAAGTGTCGCTGGAGCTGGGGGGTAACGCGCCGTTTATCTTCTTTGACGATGCCGACCCCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACt                                                                                                                                                                                                                                           <  2:18684/247‑1 (MQ=255)
  ggTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGCCCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAACCCGGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACt                                                                                                                                                                                                                                           <  1:8961/251‑1 (MQ=255)
            aaaTTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGc                                                                                                                                                                                                                                 >  1:284515/1‑251 (MQ=255)
             aaTTGGCCGCCAGTTAATGGAAAAGTGCGCGAAAGACATCAAGAAAGTGTCGCGGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGCCCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGTCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCg                                                                                                                                                                                                                                <  2:181523/251‑1 (MQ=255)
                    cgccAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCt                                                                                                                                                                                                                         >  1:168207/1‑251 (MQ=255)
                    cgccAGTTAATGGAAAAGGGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCGGGGGGGTAACGCGCCGTTTATCGTCTTTGACGAGGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCGGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCt                                                                                                                                                                                                                         <  1:313323/251‑1 (MQ=255)
                           tAATGGAACAGTGCGCAAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAAc                                                                                                                                                                                                                  >  1:71870/1‑251 (MQ=255)
                            aaTGGAACAGTGCGGGAAAGACATCAAGAAAGTGTCGCGGGGGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGGGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAACCCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACg                                                                                                                                                                                                                 <  1:275473/251‑1 (MQ=255)
                            aaTGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACg                                                                                                                                                                                                                 >  1:40714/1‑251 (MQ=255)
                            aaTGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACg                                                                                                                                                                                                                 >  1:258896/1‑251 (MQ=255)
                              tGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGc                                                                                                                                                                                                               >  2:171314/1‑251 (MQ=255)
                                      ggcgcGAAAGACATCAAGAAAGTGTCGCGGGAGCTGGGCGGTAACGCGCGGTTTATCGTCTTGGACGAGGCCGACCCCGACAAGGCCGGGGAAGGCGGGCG‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCAt                                                                                                                                                                                                       <  2:294749/250‑1 (MQ=255)
                                       gcgcGAAAGACATCAAGAAAGTGTCGCTGGAGCGGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCGGGGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATc                                                                                                                                                                                                      <  1:147802/251‑1 (MQ=255)
                                        cgcgAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCg                                                                                                                                                                                                     >  1:6164/1‑251 (MQ=255)
                                            aaaaaCATCAAGAAAGGGTCGCTGGAGCTGGGCGGTACCGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGcc                                                                                                                                                                                                 <  1:6334/247‑1 (MQ=255)
                                               aacataaaaaaaGTGTCGCTGGAGCTGGGCGTTAACGCGCCGTTTATCGTCTTTGACGATGCCGCCCTCGACAAGGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCACCCCCGGGCAAACCTGCTTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCt                                                                                                                                                                                              <  1:106823/242‑1 (MQ=255)
                                                    aaagaagttgtCGCTGGAGCTGGGCGGAACCCCGCCGTTTATCGTCTTTGACGATGCCGCCCTCGACAAGGCCGTGGAAGGCGCGCT‑GGCCTCGAAATCCCGCAACGCCGGGCAACCCTGCGTCGGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCg                                                                                                                                                                                         <  2:175591/243‑1 (MQ=255)
                                                       aaaagggtccctgaggCTGGGTGGAAACCCGCCGTTAATCGTCTTTGCCGAGGCCGCCCTCGACAATGCCGTGGAAGGCCCGCT‑GGCCTCGAAATTCCGCAGCGCGGGGCAACCCTGCGTCGGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATg                                                                                                                                                                                      <  2:12224/236‑1 (MQ=255)
                                                            tgtCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGaaaa                                                                                                                                                                                  >  2:208251/1‑250 (MQ=255)
                                                            ggtggcgggAGCGGGGCGGTACCGCGCGGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCGGTGGAGGGCGCGCG‑GGCCTCGAAATTCCGCAACGCGGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGaaaaa                                                                                                                                                                                 <  2:108752/244‑1 (MQ=255)
                                                             tgcGCTGTAGCTGGGCGTTAACGCTCTGTTGAGCGCCTTTGCCGATGCCGACCCGGACAAAGCCGTGGAAGG‑GCCCTTGTCCTCGAAATTCCCCACCGCCGGGCAACCCTGCGTCTGCGCAAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGaaaaa                                                                                                                                                                                 <  1:347724/248‑1 (MQ=255)
                                                                     aGCTGGGCGGAACCGCGCTGTTTATCTTCTTGGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAACCCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCa                                                                                                                                                                        <  1:75450/251‑1 (MQ=255)
                                                                         gggCGGTAACGCGCCGTTTATCGTTTTTGACGATGCCGCCCTCGAAAAGGCCGGGGAAGGCGCGCT‑GGCCTCAAAATCCCGCACCGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAg                                                                                                                                                                    <  2:156108/251‑1 (MQ=255)
                                                                                 aCGCGCCGTTTTTCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAagag                                                                                                                                                            <  1:356682/251‑1 (MQ=255)
                                                                                    cgcCGTTTATCGTCTTTGACGATGCCGCCCTCGAAAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCat                                                                                                                                                         <  2:64298/251‑1 (MQ=255)
                                                                                       cGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATAtt                                                                                                                                                      >  2:204181/1‑251 (MQ=255)
                                                                                          ttATCGTCTTTGACGATGCCGACCTGGACAAACCGGGGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGc                                                                                                                                                    <  1:216002/250‑1 (MQ=255)
                                                                                            aTCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGACGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCg                                                                                                                                                  >  1:143619/1‑250 (MQ=255)
                                                                                                 cTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCTTCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATgcgc                                                                                                                                            >  1:118935/1‑251 (MQ=255)
                                                                                                        gATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGgaga                                                                                                                                      >  1:345946/1‑250 (MQ=255)
                                                                                                              gcccTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGgcgc                                                                                                                               >  1:3006/3‑251 (MQ=255)
                                                                                                                 cgcGACAAGGCCGTGGAAGGCGCGCG‑GGCCTCGAAATCCCGCAACGCGGGGCAACCCGGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGgcgcgcg                                                                                                                            <  2:285217/249‑1 (MQ=255)
                                                                                                                    gACAAAGCCGTGGAAGGAGGGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGCCGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGAc                                                                                                                                                                                                                              >  2:232311/1‑150 (MQ=255)
                                                                                                                    gACAAAGCCGTAGAAGCCGCGCT‑GGCCTCGAAATCCCGCAACGCCGGGCAACCCTGCGCCTGCGCCATCCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGAc                                                                                                                                                                                                                              <  1:232311/150‑1 (MQ=255)
                                                                                                                        aaGCCGTGGAAGGCGCGCT‑GGCCTCAAAATTCCGCAACGCGGGGCAACCCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGtt                                                                                                                     <  1:346977/251‑1 (MQ=255)
                                                                                                                          gCCGGGGAAGGCGCGCT‑GGCCTCGAAATTCCGCACCGCCGGGCAACCCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTg                                                                                                                   <  2:258896/251‑1 (MQ=255)
                                                                                                                           ccGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTgc                                                                                                                  >  1:279400/1‑251 (MQ=255)
                                                                                                                                     cgcgcT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAg                                                                                                        <  1:309629/251‑1 (MQ=255)
                                                                                                                                            ggCCTCGAAATCCCGCACCGCCGGGAAAACCGGCGTCTGCGCCACCCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACg                                                                                                  <  2:284515/251‑1 (MQ=255)
                                                                                                                                              ccTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCACTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTACAGCGCACgc                                                                                                 >  2:108732/1‑249 (MQ=255)
                                                                                                                                                 cGAAATTCCGCAACGCCGGGCAACCCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAAcg                                                                                              <  1:231077/250‑1 (MQ=255)
                                                                                                                                                    aatcccGCAACGCGGGGCAAACCGGTGTTGGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACgcggc                                                                                          <  2:46401/247‑1 (MQ=255)
                                                                                                                                                            accGCCGGGCAAACCGGCGTCTGCCCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAActt                                                                                  <  2:158895/249‑1 (MQ=255)
                                                                                                                                                            aaCGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGCACGCGGCGGCAAcct                                                                                  >  2:309574/1‑249 (MQ=255)
                                                                                                                                                             aCGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAg                                                                                                                                                                                >  2:236675/1‑156 (MQ=255)
                                                                                                                                                             aCGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAg                                                                                                                                                                                <  1:236675/156‑1 (MQ=255)
                                                                                                                                                                 cGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTcc                                                                             <  1:268148/251‑1 (MQ=255)
                                                                                                                                                                    gCAACCCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGc                                                                          <  1:253974/251‑1 (MQ=255)
                                                                                                                                                                       aaCCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGa                                                                       <  1:97015/251‑1 (MQ=255)
                                                                                                                                                                         ccTGCGTCTGCGCCATCCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGAcc                                                                     <  1:237811/251‑1 (MQ=255)
                                                                                                                                                                              gtccgcgcCACCCGCCTGTATGGGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAATCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGTAGAAAGGCGCGCTCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCt                                                                <  1:333552/247‑1 (MQ=255)
                                                                                                                                                                                     ccagccGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCCCCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCtggtgg                                                           <  2:198723/245‑1 (MQ=255)
                                                                                                                                                                                              tGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGAAGAAGGCAATGAGCACACTGCACATCGGAGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGAGCAGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGAAAAGCGCACGAACGCGGCGGGAACTTCTTCCAGCCGACCCTTATGGTGGAAGTTCCgggc                                                >  1:200931/1‑249 (MQ=255)
                                                                                                                                                                                                    tgCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAACGCACGCAAGCGGCGGGCACTTCTTTCAGCCGACCCTTTTGGTGGACGTTCCCGCCAACGcc                                          >  2:361156/1‑251 (MQ=255)
                                                                                                                                                                                                                   aTGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAaagaa                           <  1:278817/251‑1 (MQ=255)
                                                                                                                                                                                                                   aTGACCGTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAaaga                            <  1:208251/250‑1 (MQ=255)
                                                                                                                                                                                                                         gTTTTGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACg                     <  1:177465/251‑1 (MQ=255)
                                                                                                                                                                                                                            ttGCCGAAAACTTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGGTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTc                  >  1:62221/1‑251 (MQ=255)
                                                                                                                                                                                                                            ttGCCGAAAAATTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCTTCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGtt                   <  2:118935/250‑1 (MQ=255)
                                                                                                                                                                                                                                   aaaaTTGCAGCAGGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGccccc           <  1:92466/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cagcagGCAATGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCCACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGcccc    >  1:270338/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gcaggcaaTGAGAAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGCCCGCTGATCGATGAAAACGCGGTAGCAAAAGTGGACGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGAAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGCCCATTCTGGTGGCCGTTCCGGCCAACGCCAAAGAGACGCAAGAAGAGACGTTCGGGCCCCTCGCCCCGc  >  2:26059/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
CCGGTTCGACCGAAATTGGCCGCCAGTTAATGGAACAGTGCGCGAAAGACATCAAGAAAGTGTCGCTGGAGCTGGGCGGTAACGCGCCGTTTATCGTCTTTGACGATGCCGACCTCGACAAAGCCGTGGAAGGCGCGCT‑GGCCTCGAAATTCCGCAACGCCGGGCAAACCTGCGTCTGCGCCAACCGCCTGTATGTGCAGGACGGCGTGTATGACCGTTTTGCCGAAAAATTGCAGCAGGCAGTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGC  >  NC_000913/2791967‑2792452

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