Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,102 T→C G123G (GGT→GGC)  thrC → threonine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,1020TC98.0% 137.8 / ‑4.5 51G123G (GGT→GGC) thrCthreonine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (18/31);  total (19/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

ATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGA  >  NC_000913/3869‑4335
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                         
aTGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCAt                                                                                                                                                                                                                          <  1:283306/251‑1 (MQ=255)
   cTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGGCCCATATTTCGGGGGCTAAGCCAGTGCCCATTCt                                                                                                                                                                                                                       >  2:299390/1‑251 (MQ=255)
     gAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCTGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCt                                                                                                                                                                                                                       <  2:29305/249‑1 (MQ=255)
             aTTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCCATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCCCCAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTATGCCCCCGAcc                                                                                                                                                                                                             >  1:285996/1‑251 (MQ=255)
                        cGCAGTGCGAAGTTCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATAc                                                                                                                                                                                                  <  2:394455/251‑1 (MQ=255)
                            aaccgaagaatctcccggcggTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGa                                                                                                                                                                                              <  1:310581/236‑1 (MQ=255)
                              gCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGATCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAg                                                                                                                                                                                             <  1:274976/250‑1 (MQ=255)
                              gCGAAGAGCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGc                                                                                                                                                                                            <  2:100661/251‑1 (MQ=255)
                                 aaGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCggg                                                                                                                                                                                         >  1:411310/1‑250 (MQ=255)
                                           cGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCAt                                                                                                                                                                               <  1:390089/251‑1 (MQ=255)
                                           cGGCGTTTATTGGTGATGAAATCCCGCAGGAAAGCCTTGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCAt                                                                                                                                                                               <  1:352932/251‑1 (MQ=255)
                                           cGGCGTTTATTGGGGATGAAAGCCCGCAGGAAATCCTGGAAGCGCGCGGGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTGCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCAt                                                                                                                                                                               <  2:202115/251‑1 (MQ=255)
                                               gTTTAGTGGTGATGAAATCCCGCAGGAAAGCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCCCCGGTCGCCAATGTTGAAAGCGATGTCGTTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCtt                                                                                                                                                                           <  2:299071/251‑1 (MQ=255)
                                                 ttATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATGGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCttt                                                                                                                                                                          <  1:370053/250‑1 (MQ=255)
                                                          atgaAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTa                                                                                                                                                                <  1:55559/251‑1 (MQ=255)
                                                                  ccgcaGGAAAGCCTGGAAGAGCGCGTGCGCGCGGCGGTGGCTTTCCCGGCTCCGGTCGCCAATGGTGAAAGCGATGGCGGTTGTCGGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGGCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAAtgt                                                                                                                                                        <  2:33580/248‑1 (MQ=255)
                                                                        gAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGCCGCTTTTTGGCCCAAAAGCTGACCCCAATTGCTGGCGGTAAGCCCGGGACCATTCTGACCGCGGCCTCCCGTGAAAACGGCGCGGCAGTGGCTCATACTTTCTACGGTTTTCCGacagtga                                                                                                                                                      >  2:119756/1‑241 (MQ=255)
                                                                            tCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGCTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGCTAAGCCAGTGACCATTCTGACCGCGACCTCCGGAGATACCGGGGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGACTGTGAAAGTGGtt                                                                                                                                              >  1:222167/1‑251 (MQ=255)
                                                                                ggAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGCATGTGAAAGTGGGTATcc                                                                                                                                          >  1:132463/1‑251 (MQ=255)
                                                                                 gAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGGGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGGGGCTCCTGCTTTCTACGGTTTTCCGCATGTGACAGTGGTTATCct                                                                                                                                         >  1:315484/1‑251 (MQ=255)
                                                                                     agCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTGTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTAt                                                                                                                                     <  1:103009/251‑1 (MQ=255)
                                                                                       ggcgTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAAGGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATAGTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATcc                                                                                                                                   <  2:378342/250‑1 (MQ=255)
                                                                                                   gcgTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAAt                                                                                                                       >  1:476358/1‑251 (MQ=255)
                                                                                                     gTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGCATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCCCAAATGCTTCCCCATCTTGCGGGCGATTAGCCAGTGACCATTCTGGCCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATACTTTCTACGGTTTACCGAATGTGCAAGTGGTTATCCTCTATCCACGGGGCAAAATc                                                                                                                      >  1:341694/1‑250 (MQ=255)
                                                                                                     gTTTGCCTGCCCGGCTCCGGTGGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGGGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGTTTATCCTCTATCCACGAGGCAAAATCa                                                                                                                     <  1:203150/251‑1 (MQ=255)
                                                                                                               ccGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCa                                                                                                           <  1:377908/251‑1 (MQ=255)
                                                                                                                     ccgggtggtaccgatgaacgctaTGTCTGTTGTCGGTAATTGTTCCACGGGCCAACGCGGGCATTTAAAGATTTCGGCGGTGGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGCCCGCCGGTGATGCCGGAGCGGCAGTGGCTCATCCTTTCTACGGTTTACCGAATGTGAAAGTGGTTAACCTCTATCCACGAGGCAAAATCAGTACCCTGCAAGaaa                                                                                                      <  2:94451/228‑1 (MQ=255)
                                                                                                                             cAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTggtat                                                                                             >  1:473795/1‑247 (MQ=255)
                                                                                                                              aaTGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCtt                                                                                            >  1:365324/1‑250 (MQ=255)
                                                                                                                                   tGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAACGATTTCGGCGGGCGCTTTATGGCACAAATGCTGGCCCATATTGCGGGCGATAAGCCAGTGACCCTTCTGACCGCGGCCTCCGGTGATACCGGGGCGGCAGTGGCTCATTCTTTCTCCGGATTACCCATTGTGAAAGGTGTTATCATCTACCCACGCGGCAAAATCAGGCCACAGCAAGAAAAACAGttttttgcct                                                                                       >  2:366352/1‑243 (MQ=255)
                                                                                                                                      agtcgttgggggTTGTCGGGAATTGTTCCACTGGACAATGTTGTCATTTAAATAGTGCGGCGGGCGCGTTATGGCACAAATGCTGACCAATATTGCGGTCGATAAGCCAGTGACCATTCTGACCGCGGCCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGGTTACCGAATGTGAAAGTGGTTATCTTCTATCCTCGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTCCATTgg                                                                                    <  2:85273/241‑1 (MQ=255)
                                                                                                                                        aagATGTCGGTTGTATGGACTGGTTCCTCGGGCCAACGCTGGCTTTTAAAGACTTCGGAGGTCGTGTTATGTCACAAATGCTGACCCATTTTGCGGGCGATAAGCCAGTGACCATTCTGACCTCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCTTGCTTTCTACGGTTTACCGAATGTGAAAGTGTTTATCCTCTATCCACGAGGCAAAACCAGTCCACTGAAAGAACAACTTTTCTGTCCCTTGggc                                                                                  <  2:143275/249‑1 (MQ=255)
                                                                                                                                         cGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGggcg                                                                                 >  1:358785/1‑251 (MQ=255)
                                                                                                                                             gTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCa                                                                              <  1:299390/250‑1 (MQ=255)
                                                                                                                                               cGGTTGTCTGTAATTGTTCCACAGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTTTGGCCCAAATGCTGACCCATATTAAGGGCGATAAGCCACTGGCCATTCTGACCACGACCTACGGGGATAACCGAGCGGGCGTGGCTCATGCTTTCTTCGGGTTACCCAATGTGGAAGAGGTTATACTTCATCTACGAGGCAAAAGCCGGCCACTACAAGACAAACTTTTTTTTACATTgggtgcgact                                                                            >  2:111910/1‑243 (MQ=255)
                                                                                                                                                ggTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCCGTGACCATTCTTACCGCTACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTATACGGTTTACCGAATGGTAAAGGGGTTATACTCTATCCCCGATGCAAAATCAGTCCACTGCCAGAAAAACTGTTCTGTACATTTGGCGGCCatat                                                                          >  2:458782/1‑251 (MQ=255)
                                                                                                                                                       gggtaTTGTCCCCCGGGCCAACCCTGGCATTTAAAGTGTTCGGTGCACGATTTATGGCACCAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCAATCTGACCGCGACCTCCGGTGATACCGTAGCGGCAGTGTCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCAACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAAGATCGaaa                                                                     <  2:389124/245‑1 (MQ=255)
                                                                                                                                                        ggAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTGCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGaaa                                                                     <  1:146264/248‑1 (MQ=255)
                                                                                                                                                                            cGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCTACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTc                                              <  1:302470/251‑1 (MQ=255)
                                                                                                                                                                            cGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCAc                                                                                                               <  1:330454/186‑1 (MQ=255)
                                                                                                                                                                              cTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCCCAAATGCTGACCCATAATGCGGGCGATAAGCCAGTGACCATTCTGAACGCGACCTCCGGTGATaacg                                                                                                                                                                                                >  2:444216/1‑100 (MQ=255)
                                                                                                                                                                              cTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCg                                                                                                                                                                                                <  1:444216/103‑1 (MQ=255)
                                                                                                                                                                                      tAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCctct                                                                                                                                       <  1:220562/152‑1 (MQ=255)
                                                                                                                                                                                       aaagagtgcggcgtTCGCGTTCTGCCACAAATGCTGACCCAGATTTCGGGCGATAAGCCCGTGACCATTCCGACCGCGCCCGCCGGTGATCCCGGAGCGGTAGTGTCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGGTCTGTACATTGGGCGGCAAGATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCGGTc                                    <  2:481812/242‑1 (MQ=255)
                                                                                                                                                                                                ggcggTCGTTTAATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCGCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGt                          <  2:120096/251‑1 (MQ=255)
                                                                                                                                                                                                 gcggTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTg                         >  1:285555/1‑251 (MQ=255)
                                                                                                                                                                                                          ttATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAgg                  <  1:466537/249‑1 (MQ=255)
                                                                                                                                                                                                          ttATGGCACAAATGCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGc                 <  2:365324/250‑1 (MQ=255)
                                                                                                                                                                                                                   gaaTGCTGACCCATATTGCGGGCGATAAGACAGTGCCCATTCTGACAGCGACCTCCGGTGATCCCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGCGGCAAGCTCAGTCCACTGCAAGACAAACTGTTCTGTACATTGGGCGGCAATATCGCAACTGTTGCCATCGACGGCGATTTCTCTGCCTGGCAGGCCCTGGTGGCGCAGGCGTTtgatgat       >  2:206459/2‑251 (MQ=255)
                                                                                                                                                                                                                       gCTGACCCATATTGCGGGCGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCGGGCGTTtgatgatga     >  2:368602/1‑249 (MQ=255)
                                                                                                                                                                                                                          gccccATATTGCGTGCGATAAGCCAGTGACCATTCTGACCGCGCCCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATgaaga  <  2:358785/247‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                         
ATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGATGTCGGTTGTCTGGAATTGTTCCACGGGCCAACGCTGGCATTTAAAGATTTCGGCGGTCGCTTTATGGCACAAATGCTGACCCATATTGCGGGTGATAAGCCAGTGACCATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGAAGA  >  NC_000913/3869‑4335

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