Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,151,689 G→T 33.0% intergenic (+43/+22) yegJ → / ← yegK uncharacterized protein/ser/thr phosphatase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,151,6890GT33.0% 41.5 / 13.7 56intergenic (+43/+22)yegJ/yegKuncharacterized protein/ser/thr phosphatase‑related protein
Reads supporting (aligned to +/- strand):  ref base G (21/16);  new base T (10/9);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

TATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTA  >  NC_000913/2151475‑2151921
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           
tATTGTTAAAAACGGAGAATAGGGCAGGGAGTGAAAAAAAAAGAAAAAAGATATTTCAGTCTGGATTTATACGGGGGGGAAAAAAATTTGCGGTGTTTATACCATTGATCTTTTATTGGTCACTTATCAAAAAAAAGAAGCCGACGAGAAAAGAGAAAAATTGGTTCGCTAACCCGCCTTTATATTTATTGTTAAAACGGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAactc                                                                                                                                                                                                                    <  2:939130/238‑1 (MQ=255)
 aTTGTTACAAACGTAGAATATGGCGAGGAGGGAAAAAAAAAAAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGAATAAAATTTACGGCGGTTATCCCTTTGATCTTTTATTGGTCACCTATCAAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGTTAACCAGCCGTTATATTTATTGTTTAAACGGGCGGTTGAATTGTTTTTTTCAACCGCCATGAGAAactcac                                                                                                                                                                                                                  <  1:842872/239‑1 (MQ=255)
                  atatGGCCGGGGGTGAAAAAAAAAAAAAGAAGATATTCCAGACTGGATTTATACGCGCGGGGATAAAATTTACGGGGGTTATACAATTAATCTTTTATTGGTCACCTACCAAAAAGAAGAAGCCGACGAGCTAAGAGAAAAATTAGTTCCTTAATCAGCCTTTATATTTATTGTTAAACCGGGCGGTTAAATTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGcc                                                                                                                                                                                                   <  1:1370609/237‑1 (MQ=255)
                     tGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTTATCCTTTTTTTGTCACCTATTCAAAAAAAGAAGACGACGAGCTAAGAGCAAAAATAGTTTGCTTATTAGACGTTATATTTATTTTTTAAAATAGCGGGTGTATTTTTTTTTTTTAACACCCTTGAAAAAACACACTTTCCCCACAGCCCCAg                                                                                                                                                                                                 >  1:404175/1‑238 (MQ=255)
                       gcagggaggaaaaaaaaaaaaaaaaaGATTTTACAATCGGGATGTAACCGGGCGGGGAAAAATTTTCGGGGGTTTATCCCATTGATCTTTTATTGGTCCCCTAACAAAAAAAAGAAGCCGACGAGAAAAGAGAAAAATTAGTTCGCAAACCAGCCGTTATATTTATTTTTAAAACGGGGGGTTGATTTTTTTTCTACAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGcc                                                                                                                                                                                             <  2:1272789/218‑1 (MQ=255)
                         cAGGAGTGGAAAAAAAAAAAGAAAGTTTTTTCAGCCGGGATGTATCCGCGGGGGAAAAAATTTACCGGGTGTTATCCCATTGTTCCTTTATGGGTCCCCTACCAAAAAGAAGAACCCAACGAGTTAAGAGAAAAATTGTTCCGCTAATAAGCGTTTATTTTTATTGTTTAAACGGGGGGTTAAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCaa                                                                                                                                                                                           <  1:501181/238‑1 (MQ=255)
                             aGTGGAAAAAAAAGAAAAAAGATTTTACACTCTGTATGTATACGCGGTGGAAAAAATTTTACGGGGTTTATCCCATTGATCTTTTTTGGTCCCCCTATCAAAAAAAAGAGGCCCGGGAGCAAAGAAAAAAATTTTTTCGTAAACCACCCTTTATATTTTTTTTTTAAACGGGGGGTGAAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGtt                                                                                                                                                                                       <  2:1073912/238‑1 (MQ=255)
                                                      ttCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCtt                                                                                                                                                                                                                                         <  1:797680/165‑1 (MQ=255)
                                                               gttGTATACGCGGGGGGAAAAAATTTCCGGGGTTTACACTGTTGATCTTTTTTGGGTAACATATCAAAAAGAAGAAGCCGACGAGTAAAAATAAAAATTATTTCCATAAACATCTTTTATATTTATTTTTAAAACTGGGGGTTTAATTCTTTCCCTAATCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGc                                                                                                                                                     <  1:1296637/236‑1 (MQ=255)
                                                                       cgcgggGGGATAAAATTTACGGGGTTTATCCCATTGATCTTTTATTGGTCACTTATCAAAAAGAAGAAGCCGACGATCTAAGAGAAAAATTATTTCGCAAATCAGCCGTTATATTTTTTGTTAAACCGGGGGGTTAATTTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACaaaa                                                                                                                                            <  1:507052/239‑1 (MQ=255)
                                                                       cgcgggGGAAAAAAATTTACGGCGGTTATACCTTTTATCCTTTTTGGCTCCCCTCCCAAAAAAAAAAAGCCGCCGATCAAAGAGAAAAATTGGTTCGCTAATCAGCCTTTATATTTTTTTTTTAACCGGGGGGTTGATTTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACaaaa                                                                                                                                            <  1:1475065/239‑1 (MQ=255)
                                                                                       ttACGGGGTTTATACATTTGATCTTTTATTGGTCCCCTATCAAAAAGAAAAAGCCGAAGAGCTAAGAGAAAAATTATTCCGCTAATCAGCGTTTATTTTTTTTTTTTAAACGGGCGTTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTc                                                                                                                           <  1:1417096/239‑1 (MQ=255)
                                                                                       ttACGGGGGTTAAACAATTGATCTTTTTTTGGTCCCCTACCAAAAAAAAGAAGCCGACGAGAAAAGAAAAAAATTATTCCGTTAACCAGCTGTTATTTTTATTTTTAAAACGGGGGTTTGATTTTTTGTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTc                                                                                                                           <  1:762124/239‑1 (MQ=255)
                                                                                        tacgggaggTTATAACTTTGATCTTTTATTGGTCCCCTAACAAAAAAAAAAACCCGACGAGTAAAGAAAAAAATTAGTTCGCTAACCAGCCTTTATTTTTTTTGTTAAACCGGGCGGTTGAATTGTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTc                                                                                                                           <  1:30365/232‑1 (MQ=255)
                                                                                                 tttATACCTTTTATCCTTTTTTGGTCCCCTATCAAAAAAAAAAACCGAACGAGCTAAGAAAAAAATTTGTTCGTTAATCACCCTTTTTTTTTTTTTTTAAAACGGGGGTTTGAATTTTTGTTTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGtaat                                                                                                                    <  1:1409538/237‑1 (MQ=255)
                                                                                                         aTTTATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGGTTTAATTTTGTTTTTCAACCCCCTTAGGAAAACAACCTTCCCCCCACGCCCCGGGCCACCGTTTTTTCATATACCTTACCCCCGTTTCCCCCGCGGCTGGGCAAAAACTCTTTTAAAAATTCCgggataaaaat                                                                                                               >  2:254185/1‑227 (MQ=255)
                                                                                                           tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGGATTTTTTTTTTCCACCCGCATTAAAAAATCACTCTTCCCACAACGCCAGTGCCAACGGTTTATCATCATCCCAACCCCCGCTTACCGCCGGGCTGGGCAAAAACTTCTTTTACCATTTCAGGAATTAAATAGAgtt                                                                                                          >  2:1125891/1‑239 (MQ=255)
                                                                                                           tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTTAACCGCCATGAGAAAATAACTCTGCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCCCCTTTACCCCCGGGCTGGGCAAAAACTGCTTTAACAATTCAGGGAATAAATCGAgtta                                                                                                         >  2:87208/1‑238 (MQ=255)
                                                                                                           tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCTACCGCCAAGAGAAACTCACTCTCCCCACCGCGCCACTGCCCACGTTTTATCATCATCCCTACGCTCGTTTAACCCCGCGCTGGGCCACAACCTCTTTTACAATTAAGGTAATAAAACCAgttt                                                                                                         >  2:1402626/1‑238 (MQ=255)
                                                                                                            gATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAATTGTTTTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACCTTTTTTCATCATCCGTATGCCCCCTTACCGCCCGGCTGGACCAAAACTGCTTTTACAATTTAGGTCATAAATCGtgttgt                                                                                                        >  1:842855/1‑239 (MQ=255)
                                                                                                            gATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAACGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttgt                                                                                                        >  1:1218018/1‑239 (MQ=255)
                                                                                                            gATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCACGGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGGCAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgttgt                                                                                                        >  1:620733/1‑239 (MQ=255)
                                                                                                             aTCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCCACCCCCCTGAGAAAATCACTCTGCCCCCAGCGCCAGTGCCAACGTTTTATCATCATCCCTACGCTCGTTTACCGCCGGGCTGGACAAAAAATGCTTTAACAATTCCGGTAATAAATCGAgttggt                                                                                                       >  2:993269/1‑237 (MQ=255)
                                                                                                               ccTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTGTGTTTTCAACCGCCATGAGAAAACCCCTCTGCCCACAGCGCCCGTGCCAACGGTTTATCATCATCCGGACGCTCGTTTACCCCCGGGCTGGGCCAAAACTGCTTTAACAAATCAAGGAATAAATCGAggtttt                                                                                                       >  2:586337/1‑232 (MQ=255)
                                                                                                                  ttATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGGATTGTTGTCTTCAACCCCCCTGAGAAACTCACTCTGCCCCCAGCGCCCGGGCCAACGTTTTATCATCATCCCTACCCCCGTTTTCCGCCGGGGCGGGCAAAAACCGCTTTTACAATTTAGGTAATaaaaccggcttttcctt                                                                                                   >  2:1103948/1‑224 (MQ=255)
                                                                                                                   tATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTGATTGTTTTCTTCCACCGCCATGAGAAACTCCCTCTGCCCCCAGCGCCCGTGCCCACGCTTTTTCATCATCCGTCCCCTCGTTTTCCCCCCGGCTGGCCAAAAACTGCTTTTACAATTTCGGTAATaaaaccggttttttcttc                                                                                                  >  2:490281/1‑223 (MQ=255)
                                                                                                                    aTTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGCATTTTTTTCTTCAACCGCCCTGAGAAACTCACTCTGCCCCCAGCGCCAGCGCCCACGCTTTATCATCCTCCGCACCCCCGTTTACCGCCCGGCTGGCCAAAAACTGCTCTAACAATTCACGCAATAAATCGCGTTCTTCCTGCgg                                                                                                >  2:710382/1‑238 (MQ=255)
                                                                                                                        gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGGATTGTTGTCTTCAACCGCCCTTAGAAACTCACTCTGCCCAACGCGCCAGTGCCAACGTTTTATCACCCTCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTTACAATTCTGGGAATAAATCGAGTTGTTCCTGCGtt                                                                                               >  2:360203/1‑236 (MQ=255)
                                                                                                                        gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTTTTTCACCGCCGTAAGGAACTCACTCTTCCCACAGCGCCAGTGCCAACGTTTTTTCAACCACCGTACCCTCGTTTTCCGCCGGGCCGGGCAAAAACCGCTTTATCAATTTAGGTCATAAAACGAGCTGTTCCTGCgattca                                                                                            >  2:64607/1‑234 (MQ=255)
                                                                                                                          caataaataaaaaaaaaaaaGCCGAGAAGAAAAGAAAAAAATTTGTTCGATAAGCATCTGTTATATTTATTGTATAAACGGGGGGTTAGATTTTTGTCAAAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgct                                                                                            <  2:620733/228‑1 (MQ=255)
                                                                                                                            aCCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCTACCGCCATTAGAAACTCACTCTGCCCACAACGCCCATGCCAACGCTTTTTCATCATCCGTAAGCCCGTTTTCCGCCCGGCTGGGCCAAAACTTCTTTAACCATTTAGGTAATAAATCGAGTTGTTTCCTCGTtgctgc                                                                                          >  2:488447/1‑238 (MQ=255)
                                                                                                                            aCCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCCTGAGAAACTCACTCTGCCCCCCGCGCCAGTGCCCACGTTTTATCATCATCCGTACGCTCGTTTTCCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgctgc                                                                                          >  1:1334267/1‑238 (MQ=255)
                                                                                                                            aCCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCCCAGCGCCAGTGCCAACGGTTTTTCATTATCCGTAAGCCCGTTTACCCCCGGGCCGGACAAAAACTGCTTTTACAATTCAGGGAATAAATCCAGTTTTTCCTGCGtttgtgc                                                                                          >  1:274662/1‑233 (MQ=255)
                                                                                                                              cTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTTTCTTCAACCCCCCTTAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCTCGTTTATCGCCCGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGGTGTTCCTGCGTTGCTgccg                                                                                        >  1:1198077/1‑238 (MQ=255)
                                                                                                                              cTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAAATGTTGTCTTCAAACCCCATGAGAAAATCAATCTTCCCACAGCGCCAGTGGCAACGTTTTTTTATCATCAGTAACCCCCGTTAACGCCAGGGGGGGCAAAAACTGCTTTAACAATTTAGGAACTAAATCAGGTTGTTCCTGCGTTGCTgccg                                                                                        >  2:617755/1‑238 (MQ=255)
                                                                                                                                aTCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAAATGGTGTCTTCAACCGCCCTGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCCCGTTTACCGCCGGGCTGGGGAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCa                                                                                     >  1:1260394/1‑239 (MQ=255)
                                                                                                                                             agCCGACGAGCAAGGGTAAAAATTAGTTCGAAAATCAGCCTTTATTTTTTTTGTTAAACCGGGCGTTTAAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGaaaaa                                                                        <  1:1021560/239‑1 (MQ=255)
                                                                                                                                                   cgaGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTTACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCCTTGAAAAACGGGGt                                                                  >  1:494350/1‑239 (MQ=255)
                                                                                                                                                   agattaaaaaaaaaatTTAGTTCTATAATCAGCGTTTATTTTTTTTTTAAAAACGGGGGGTTAATTTTTTGTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACgggg                                                                   <  1:993269/226‑1 (MQ=255)
                                                                                                                                                        aaagaaaaaaaTTATTTCGCTAATCAGCCTTTTTATTTATTTTTAAAACTGGGGGTTAAATTTTTTTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTaaaga                                                             <  1:1101029/232‑1 (MQ=255)
                                                                                                                                                            agCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGCGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAgg                                                         >  1:568244/1‑239 (MQ=255)
                                                                                                                                                                    tattgcttaattgatccGTTATTTTTTTTTTTAAACTGGGAGTTGGATTTTTTGTCTCAACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                 <  1:1402626/224‑1 (MQ=255)
                                                                                                                                                                    aggTTCGTTACACACCCGTTATTTTATTTGTTAAAACGGGGGGTGAATTGTTGGCTTCAATCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTGATCATCATCCGTACGCGCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTGGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                                 <  1:1101241/237‑1 (MQ=255)
                                                                                                                                                                      gTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCCGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCCGGCTGGACAAAAACTGCTTTAACAATTCCGGTAATAAATCGGGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTTAAAAACGGGGTAAAGAAAAGCACATgaggg                                               >  1:723251/1‑238 (MQ=255)
                                                                                                                                                                         cgaaaacgggcgtatataTTTTTGTTAAAAAGGGGGGGTGAAATTGTTGTCTACACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAAGCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGaa                                            <  1:490281/225‑1 (MQ=255)
                                                                                                                                                                         aaaagaaggccctttttattttttTTTAAACCGGGCGGTTAATTTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGaa                                            <  1:93361/223‑1 (MQ=255)
                                                                                                                                                                            tAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTTTTTTTTTCAACCGCCATGAGAAAACCACTCTTCCCACAGCGCCCGTGCCAACGTTTTATCTTTATTCCGACGCCCGGTTATCACCGGGGTGGACAAAAAATGCCTTTACAATTCAGGTAATAAATCGGATTGTTTCCTCATCGCTGCCGCCCGTCCATTGAAAAAAGGGGTAAAGCAAGGCACATAAGGGGAATTa                                         >  2:409412/1‑239 (MQ=255)
                                                                                                                                                                                  gCCTTTATATTTTTTTTATACACGGGGGGTTGATTGTTTTTCCTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                                   <  1:488447/239‑1 (MQ=255)
                                                                                                                                                                                             tattGTTTAAACTGGCGGTTTAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCCGTGCCAACGTTTTATCATCATCCCTACGCTCGTTTACCCCCCGGCTGGGCAAAAAATTCTTTTACAATTCAGGGAATAAATCGAGTTGTTTCTGCGTTGCCGCCCCCAGGCCATTTAAAAAAGGGGGAAAAAAAAGCACCAGAAGGGGAATATTCAGCATATTTaaag                         >  2:1265149/1‑236 (MQ=255)
                                                                                                                                                                                             tattGTTTAAACTGGCGGTTGAATTGTTGTCTTTAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCCGGGCCAACGTTTTTTCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTTACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCCGTCCATTGAAAAACGGGGGAAAGAAAGGCACATGAGGAGGATTATCCAGCATATTTAACg                         >  1:1478506/1‑238 (MQ=255)
                                                                                                                                                                                                tGTTTAAACTGGCGGTTGAATTTTTTTCTTCAAACGCCATGAGAAACCCACTCTGCCCACAGCGCCAGTGCCAAAGTTTTTTCCTCATTCCTACGCTCGTTTACCCCCCGGCTGGACAAAAAATGCTTTAATAATTCCGGTAATAAAATGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCAATGAAAAAAAGGGTAAAGAAAAGCACATTAGGAGAATTATCCAGGATATTTTACGACaa                     >  2:1177603/1‑238 (MQ=255)
                                                                                                                                                                                                tGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGGTTTATCATCATCCGGACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTTACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCCTTGAAAAACCGGGGAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTactttc                       >  1:1369070/1‑232 (MQ=255)
                                                                                                                                                                                                 tttttttaacttggggttatatttgtttTTTTCATCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGcc                       <  1:1125891/217‑1 (MQ=255)
                                                                                                                                                                                                    tAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTTACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCg                  >  1:1260263/1‑238 (MQ=255)
                                                                                                                                                                                                         gggCGGTTAAATTGTTTTCTACAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGAt            <  1:87208/238‑1 (MQ=255)
                                                                                                                                                                                                           gCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCCGTGCCAACGTTTTTTCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGGACATtagga                                               >  1:1372737/1‑202 (MQ=255)
                                                                                                                                                                                                                      ttttttgtTTCCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTa  <  1:818027/235‑1 (MQ=255)
                                                                                                                                                                                                                      aTTGTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTa  <  1:321336/239‑1 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                                           
TATTGTTACAAACGTAGAATATGGCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATCCGTA  >  NC_000913/2151475‑2151921

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