Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,162 C→T H255Y (CAC→TAC)  spoT → bifunctional (p)ppGpp synthase/hydrolase SpoT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,823,1620CT98.3% 228.7 / NA 60H255Y (CAC→TAC) spoTbifunctional (p)ppGpp synthase/hydrolase SpoT
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (37/22);  minor base G (0/1);  total (37/23)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATG  >  NC_000913/3822918‑3823398
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
aGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATGCAGAAGATTCTGTCTGAAATCGAAGGGCGTTTGCAGGAATCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAATCATCTTTATTCGATTGACTGCAAATTTGTGCTCAAAGATCAGAGTTTTGACTCg                                                                                                                                                                                                                                          <  2:299833/250‑1 (MQ=255)
                            gAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGAt                                                                                                                                                                                                             >  1:67182/1‑251 (MQ=255)
                                 gCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCa                                                                                                                                                                                                        <  2:352586/251‑1 (MQ=255)
                                          tGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGAGCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCt                                                                                                                                                                                                <  2:84363/250‑1 (MQ=255)
                                               cgcTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATTCCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATTGACATCTACGCTGTCCGCGTGATCGTCAATGATTCTGACAcc                                                                                                                                                                                          <  1:360878/251‑1 (MQ=255)
                                                              gTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTg                                                                                                                                                                          >  2:20322/1‑251 (MQ=255)
                                                                 aTCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCaa                                                                                                                                                                                                       >  2:359623/1‑220 (MQ=255)
                                                                 aTCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCaa                                                                                                                                                                                                       <  1:359623/220‑1 (MQ=255)
                                                                      gTAAGCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAATGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACGCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGTCCAGAt                                                                                                                                                                  <  2:290393/251‑1 (MQ=255)
                                                                         aTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGc                                                                                                                                                                >  1:268584/1‑250 (MQ=255)
                                                                               gaagTGGTGAACGCCGCCAGAGGCAACCGAAAAGCGATGATCCAGCAGATTCTTTCTGAAATCGACGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGCGCAGCATCTTTATTCGCTTTACTGCAAAATGGCGCTCCAAGAGCAGCGTTTTTACTAGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCCGATGCACAGCCt                                                                                                                                                         >  2:428387/1‑251 (MQ=255)
                                                                                       gAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGc                                                                                                                                                 <  1:227918/251‑1 (MQ=255)
                                                                                         aaGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCCGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGTTTTCCTGCAAAATGGGGCTCAACGAGAAGCGTTTTTATTCGATCATGGACATCTACGCTTTCCGAGTGATCGTCAATGATTTTGACACC‑GTTTACGCCGTGTTGGGCCAGATgacaacccggtacaaccgc                                                                                                                                               >  1:410310/1‑233 (MQ=255)
                                                                                            ccGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGt                                                                                                                                            >  1:113264/1‑251 (MQ=255)
                                                                                                 cgcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCggg                                                                                                                                       >  2:197078/1‑251 (MQ=255)
                                                                                                   cgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGATGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGcc                                                                                                                                     >  1:90492/1‑251 (MQ=255)
                                                                                                             aaaGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaa                                                                                                                              >  1:344534/1‑248 (MQ=255)
                                                                                                               agagATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACt                                                                                                                         >  2:335816/1‑251 (MQ=255)
                                                                                                               agagATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACt                                                                                                                         >  1:343067/1‑251 (MQ=255)
                                                                                                                 agaTGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGCAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACtat                                                                                                                       >  1:403237/1‑251 (MQ=255)
                                                                                                                 agaTGATCCAGAAGATGCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACtat                                                                                                                       <  2:209385/251‑1 (MQ=255)
                                                                                                                   atgatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGCCTatat                                                                                                                     >  2:39337/1‑251 (MQ=255)
                                                                                                                   atgatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTatat                                                                                                                     <  1:239356/251‑1 (MQ=255)
                                                                                                                       tcaagaagaTTCTTACTGAAATCGAAGGGGGTTTGCAGGAAGCGGGAATAGCGGGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATCAGACTTACTGCAAAATGGTGCTCAAATGGCAGCGTTTTTACTCGATCATGGACATCTCCGCTTTCCGCTTGATCGTAAATGATTCTGACTCCTGGTATCGC‑TTGCGGGGCCAGAAGCACAGCGTGTACAATCCGCTTCCGGGCCGCGTGAAAGACTATAACGc                                                                                                                  <  2:243517/247‑1 (MQ=255)
                                                                                                                       tCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGcc                                                                                                                 >  1:185159/1‑251 (MQ=255)
                                                                                                                        ccAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACATGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGcc                                                                                                                 <  1:158237/250‑1 (MQ=255)
                                                                                                                            aagaTTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTcc                                                                                                            <  2:343067/251‑1 (MQ=255)
                                                                                                                             aaaTTCTTGCTGAAAGCGAGGGGCGTTTGCAGGAATAGTGAATCCCGTGCCGCGTCAGGGGTCGCGAGAAGCAGCTTTATTCGATGTACAGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGGGATCGTCAATGATTCTGACCCCTGTTACCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGTCGCATGAAAGACTATATCGCCATTCCa                                                                                                           <  2:364634/249‑1 (MQ=255)
                                                                                                                               aTTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCACAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACAATATCACCATTCCaa                                                                                                          >  2:259982/1‑250 (MQ=255)
                                                                                                                                ttctttctGAACTCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGAGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCACTGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCAGGGCCGCGTGCAAGCCTATATCGCCATTCCaaa                                                                                                         >  2:173856/1‑250 (MQ=255)
                                                                                                                                    ttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCg                                                                                                     >  1:106955/1‑250 (MQ=255)
                                                                                                                                        gAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACgg                                                                                                >  1:425044/1‑251 (MQ=255)
                                                                                                                                         aaaTCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGc                                                                                               >  1:379490/1‑251 (MQ=255)
                                                                                                                                                  gtcGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCGCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTAGCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCTCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCt                                                                                      <  1:286531/249‑1 (MQ=255)
                                                                                                                                                         gCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGcac                                                                                >  1:162846/1‑250 (MQ=255)
                                                                                                                                                            ggAATCGGGGATACCGTGCCGCGACAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTTACACCTGTGATCGC‑GTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCATGAAAGACTATATCGCCATTCCCAAAGCGAAAGCCTATCAGTCTATGCACAccg                                                                            >  2:241672/1‑250 (MQ=255)
                                                                                                                                                                  gggAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCCAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATAAGTCTTTGCACACCTCgatg                                                                       >  2:356105/1‑250 (MQ=255)
                                                                                                                                                                  gggAATACCGTGCCGCGTCAGTGGTCGAGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgatg                                                                       >  2:233526/1‑250 (MQ=255)
                                                                                                                                                                        aCCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCGCGATGATCGGcc                                                                <  1:363516/251‑1 (MQ=255)
                                                                                                                                                                                    cAGTGGGCGCGAGAAGCATCTTTATTCGAGTTACTGCAAAAGGGTGCTCCAAGAGCAGTGTGTTTATCCGATCATGGACATCTACGCGTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGtgtg                                                     <  2:18465/250‑1 (MQ=255)
                                                                                                                                                                                     aGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTTCTCAGAGATCAGCGTTTTTACGCGAGCATGGACTTCTCCGCTCTCCGCGTTGTCGTCAATGATTCGGACCCCTGTTATCGC‑GTGCTGGGCCATATGCACAGCTTGTACAAGCCTCGTCCGGGCCGCGTGAAAGACTATAGCGCCATTCCAAAAGCGAACGGCTAGCAGTCTTTGCACCCCTCGATGATCGGCCCGCACGGTGTGct                                                   <  2:131607/251‑2 (MQ=255)
                                                                                                                                                                                        ggTCGCGAGAAGCATCTTTAGTCGATTAACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACAGTGTGCCGGt                                                <  2:3354/251‑1 (MQ=255)
                                                                                                                                                                                              gagaAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGt                                          >  1:219229/1‑251 (MQ=255)
                                                                                                                                                                                                 aaGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTcc                                        >  2:301408/1‑250 (MQ=255)
                                                                                                                                                                                                  aGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTg                                                                                                                                                                                               >  1:16067/1‑99 (MQ=255)
                                                                                                                                                                                                  aGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTg                                                                                                                                                                                               <  2:16067/99‑1 (MQ=255)
                                                                                                                                                                                                      tCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATcc                                  >  2:601/1‑251 (MQ=255)
                                                                                                                                                                                                         ttATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATAGCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCGCGATGAGCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGt                                <  2:62821/250‑1 (MQ=255)
                                                                                                                                                                                                           aTTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGTGGTCCAGCTCCGTAcc                             >  2:322203/1‑251 (MQ=255)
                                                                                                                                                                                                             tCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGcac                                                                                                                                                              <  1:32275/120‑1 (MQ=255)
                                                                                                                                                                                                             tCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGcac                                                                                                                                                              >  2:32275/1‑120 (MQ=255)
                                                                                                                                                                                                                aTTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGCTCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGa                           >  2:358895/1‑248 (MQ=255)
                                                                                                                                                                                                                     cTGCAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATgg                   >  2:291341/1‑251 (MQ=255)
                                                                                                                                                                                                                        cAAAATGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGAcc                <  1:417585/251‑1 (MQ=255)
                                                                                                                                                                                                                           aaTGGTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGCTCCGTACCGAAGATATGGAcc                >  1:284789/1‑248 (MQ=255)
                                                                                                                                                                                                                              ggTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATCCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATgg          <  1:335816/251‑1 (MQ=255)
                                                                                                                                                                                                                               gTGCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGc         >  2:206343/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gCTCAAAGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGAGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCgg       >  1:178551/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cTCAAAGAGAAGCGTTCTTACTCGATCATGGACATCTACGCTTTCCGCGCGATCGTCAATGATTCTGAAACCTGTTATCGC‑GTGCTGGGCAAGATGCACAGCCTGTACACGCCGCGTCCGGGCCGCGTGAAAGACTATCTAGCCATTCCCAAAGCGAACGGCTATCCGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACAAGAGGGCggg      >  2:422564/1‑250 (MQ=255)
                                                                                                                                                                                                                                      aaGAGCAGCGTTTTTACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATg  >  1:23109/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
AGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGC‑GTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGACCAGATGGCGGAGATG  >  NC_000913/3822918‑3823398

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