Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,913,050 G→A L207F (CTT→TTT)  glmS ← L‑glutamine‑‑D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,913,0500GA98.4% 216.1 / ‑4.2 62L207F (CTT→TTT) glmSL‑glutamine‑‑D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (38/23);  total (39/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

GTCGGCGTTCGGTCCCAGCTCGCTTAAATCAACCTGACCGTGGCTGATGCGTCCGGTAAGGGTGTTTTTGATCGCGTTCGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGATAGCACGCAGAACGGCCTCACGCAGAGTCCCGCCTTGTTTCAGCTCCCAGTTCACCAGATGGGCAATCACTT  >  NC_000913/3912808‑3913292
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
gtcgGCGTTCGGTCCCAGCTCGCTTAAATCAACCTGACCGTGGCTGATGCGTCCGGTAAGGGTGTTTTTGATCGCGTTCGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATa                                                                                                                                                                                                                                             >  1:151988/1‑250 (MQ=255)
                 gCTCGTTTAAATCAACCTGACCGTGGCTGATGCGCCCGCTCAGGGTGTTTTTGATCGCGTTCGGCGGTTCGTAGATCGCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCCATATCCTGACGTTTTACTTCCGCGCCGGTTTTATCGACGAGGTTTACCGACCGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGAAGGGTCACCGGc                                                                                                                                                                                                                           >  2:348316/1‑251 (MQ=255)
                     gCTTAAATCAACCTGACCGTGGCTGATGCGTCCGGTAAGGGTGTTTTTGATCGCGTTCGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAAc                                                                                                                                                                                                                        >  2:300801/1‑250 (MQ=255)
                        tAAATCAACCTGACCGTGGCTGATGCGTCCGGTAAGGGTGTTTTTGATCGCGTTAGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAg                                                                                                                                                                                                                      <  1:86688/249‑1 (MQ=255)
                            tCAACCTGACCGTGGCTGATGCGTCAGGCAAGGGCGATTTTGATCGCGTTCGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTACGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGcc                                                                                                                                                                                                                  >  2:175298/1‑249 (MQ=255)
                             cAACCTGACCGTGGCTGATGCGTCCGGCAAGGGTGTTTTTGATCGCGTTCGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGc                                                                                                                                                                                                               >  1:103042/1‑251 (MQ=255)
                               aCCTGACCGTGGCTGATGCGTCCGGTAAGGGTGTTTTTGATCGCGTTCGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAAAAGCGCCAGCTg                                                                                                                                                                                                             >  2:172431/1‑251 (MQ=255)
                                ccttaCCGTGGCTGATGCGTCCGGTAAGGGTGTTTTTGATCGCGTTCGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGATCGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTgg                                                                                                                                                                                                            <  2:352308/247‑1 (MQ=255)
                                   gACCGTGGCTGATGCGTCCGGTAAGGGTGTTTTTGATCGCGTTCGGCTGTTCGTAGATCTCTTTCTGCATGGAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCa                                                                                                                                                                                                         <  2:402203/251‑1 (MQ=255)
                                      cGTGGCTGATGCGTCCGGTAAGGGTGTTTTTGATCGCGTTCGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGa                                                                                                                                                                                                       >  2:200423/1‑250 (MQ=255)
                                              aTGCGTCCGGTAAGGGTGTTTTTGATCGCGTTCGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGACAGATGGTCAGAAGCGATaa                                                                                                                                                                                               >  1:212762/1‑250 (MQ=255)
                                                gCGTCCGGTAAGGGTGTTTTTGATCGCGTTCGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGt                                                                                                                                                                                            <  1:156495/251‑1 (MQ=255)
                                                    ccGGTCAGGGTGTTTTTGATCGCGTTCGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTNTTc                                                                                                                                                                                        >  2:358978/1‑251 (MQ=255)
                                                    ccGGTACGGGTGTGTTTGATCGCGTGCGGCTGTTCGTAGATCGCTTTTTGCATGTAGTGACGGTAAATGCCTTTATCGCCAGCATCATATTGCAGATTGGATTCGAGATCCTGACGTTTTGCTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAAGATCGCGCTCTTCAAAGAAGATAAAGAGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCAATAAAGTTTTc                                                                                                                                                                                        <  2:350737/251‑1 (MQ=255)
                                                          aGGGTGTTTTTGATCGCGTTCGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCAt                                                                                                                                                                                  <  2:337683/251‑1 (MQ=255)
                                                          aGGGTGTTTTTGATCGCGTTCGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGAAGGGTCACCGGCAACAGCGCCAGCTGGACAGAAGCGATAAAGTTTTCGCCCAt                                                                                                                                                                                  >  1:161857/1‑251 (MQ=255)
                                                              tgtTTTTGATCGCGTTCGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGACTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGCAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATCAAGTTTTCGCCCATcccc                                                                                                                                                                              >  2:289003/1‑251 (MQ=255)
                                                                           gTTCGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCAcc                                                                                                                                                                  >  1:227752/1‑250 (MQ=255)
                                                                              cGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGc                                                                                                                                                               <  2:13012/250‑1 (MQ=255)
                                                                                   gTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGGTTTCGCCCATCCCCAGGCCAATCACCAGCGGACt                                                                                                                                                          <  1:346854/250‑1 (MQ=255)
                                                                                          gATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAAc                                                                                                                                                  >  2:156739/1‑251 (MQ=255)
                                                                                               ctTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGCCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTgcc                                                                                                                                             >  1:381292/1‑251 (MQ=255)
                                                                                                     gCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAATAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGaa                                                                                                                                                   <  2:174784/239‑1 (MQ=255)
                                                                                                     gCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGaa                                                                                                                                                   >  1:174784/1‑239 (MQ=255)
                                                                                                           aGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTa                                                                                                                                 <  2:151988/251‑1 (MQ=255)
                                                                                                              gACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATcc                                                                                                                              >  2:167558/1‑251 (MQ=255)
                                                                                                               aCGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCg                                                                                                                             <  1:156739/251‑1 (MQ=255)
                                                                                                                   tAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTg                                                                                                                         >  1:22821/1‑251 (MQ=255)
                                                                                                                   tAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGCATCCGGGTg                                                                                                                         >  2:114162/1‑251 (MQ=255)
                                                                                                                     aaTGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGCATCCGGGTGAc                                                                                                                       >  1:389658/1‑251 (MQ=255)
                                                                                                                             tATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCgg                                                                                                                                                                                                                            <  2:137859/142‑1 (MQ=255)
                                                                                                                             tATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGCTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCgg                                                                                                                                                                                                                            >  1:137859/1‑142 (MQ=255)
                                                                                                                                       gTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTg                                                                                                     >  2:99072/1‑251 (MQ=255)
                                                                                                                                        tCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGCTTTTCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGTGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGAAGGGCATCCGGGTGACGGGAGACCATGATCACTGt                                                                                                    >  2:184545/1‑251 (MQ=255)
                                                                                                                                             ttGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGt                                                                                               >  1:13801/1‑251 (MQ=255)
                                                                                                                                                cAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACg                                                                                            <  1:167558/251‑1 (MQ=255)
                                                                                                                                                         ttCGATATCCTGACGTTTTACTTCAGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACt                                                                                                                                                          >  2:90128/1‑180 (MQ=255)
                                                                                                                                                         ttCGATATCCTGACGTTTTACTTCAGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACt                                                                                                                                                          <  1:90128/180‑1 (MQ=255)
                                                                                                                                                           cGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAg                                                                                  <  2:56481/250‑1 (MQ=255)
                                                                                                                                                             atatCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTg                                                                               <  1:114162/251‑1 (MQ=255)
                                                                                                                                                               atCCTGACGTCTTCCTTCAGCGCAAGATTTATCGAAGATGTTTCCCGAACGACACGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGCTCAGAAGCGCTAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGTACTACCAGAACATGCCGCCAGCAGGGTAACCGGGTGACGGCAGACCAAGATCACTGTACCCTACGCACCACGCCGCTGCg                                                                             >  2:56171/1‑251 (MQ=255)
                                                                                                                                                                  cTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCCCGCAGCTGCggg                                                                           >  2:294531/1‑250 (MQ=255)
                                                                                                                                                                    gACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCTCCACGCAGCTGCGGGATa                                                                        <  1:300801/251‑1 (MQ=255)
                                                                                                                                                                    gACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGAt                                                                         >  2:155133/1‑250 (MQ=255)
                                                                                                                                                                     aCGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAG‑‑GGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGATAGc                                                                      >  1:47855/1‑250 (MQ=255)
                                                                                                                                                                       gTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGACGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAACTGGTCAGAAGCGATAAAGTTTTCGCCCATCCACAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGAATCCGGGTGACGGGAGTCCATGATCATTGTACCGTACGCACCACGCCGCTGCGGGATAGcc                                                                     >  2:390892/1‑250 (MQ=255)
                                                                                                                                                                        ttttACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCGTCAAAGAAGATAAAGCGACGGGTCACCGTCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGATAGCa                                                                     <  2:254900/250‑1 (MQ=255)
                                                                                                                                                                          ttACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGCTAGCACGc                                                                  >  2:120037/1‑251 (MQ=255)
                                                                                                                                                                              ttCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGATAGCACGCAGa                                                               >  2:207611/1‑250 (MQ=255)
                                                                                                                                                                                 cgcgCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGATAGCACGCAGAACgg                                                           >  1:125929/1‑251 (MQ=255)
                                                                                                                                                                                     ccAGTTTTATCGAAGCTGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTTAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGATGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGATAGCACGCAGAACGGCCTc                                                       >  1:231938/1‑251 (MQ=255)
                                                                                                                                                                                                gAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGATAGCACGCAGAACGGCCTCACGCAGAGTcc                                            <  2:284683/251‑1 (MQ=255)
                                                                                                                                                                                                gAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGATAGCACGCAGAACGGCCTCACGCAGAGTcc                                            >  1:285117/1‑251 (MQ=255)
                                                                                                                                                                                                   gATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGACAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTACCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGATAGCACGCAGAACGGCCTAACGCAGAGTCCCGc                                         >  2:123918/1‑251 (MQ=255)
                                                                                                                                                                                                            cGAACGGCGAGTGATGTCCGCAATATCGCCCTCGTCAAAGAAGATAAAGCGACGGGTGACCGGCAACAGCGCCAGCTGGTCAGAAGCGACAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAGCGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTGCCGTACTTACGACGCAGCCGCGGGATAGCACGCAGAACGGCCTCACGCAGAGTCCCGCCTTGTTTCa                                <  1:219023/251‑1 (MQ=255)
                                                                                                                                                                                                              aaCGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCcagc                                                                                                                                       <  2:336001/146‑1 (MQ=255)
                                                                                                                                                                                                              aaCGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCcagc                                                                                                                                       >  1:336001/1‑146 (MQ=255)
                                                                                                                                                                                                                         aTTTCCGCAATATCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGATAGCACGCAGAACGGCCTCACGCAGAGTCCCGCCTTGTTTCAGCTCCCAGTTCAc                   <  1:38903/251‑1 (MQ=255)
                                                                                                                                                                                                                                   tatCGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGATAGCACGCAGAACGGCCTCACGCAGAGTCCCGCCTTGTTTCAGCTCCCAGTTCACCAGATGGGc          >  2:283045/1‑250 (MQ=255)
                                                                                                                                                                                                                                      cGCCCTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGATAGCACGCAGACCGGCCTCACGCAGAGTCCCGCCTTGTTTCAGCTCCCAGTTCACCAGATGGGCAATc      >  2:264012/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cccTCTTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCCGGCCAACCACCAGCGGACTCCCCGAACGTTCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGATAGCACGCAGAACGGCCTCACGCAGAGTCCCGCCTTGTTTCAGCTCCCAGTTCACCAGATGGGCAATCAc    <  2:45638/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ctctTCAAAGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGATAGCACGCAGAACGGCCTCACGCAGAGTCCCGCCTTGTTTCAGCTCCCAGTTCACCAGATGGGCAATCACtt  >  2:304136/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
GTCGGCGTTCGGTCCCAGCTCGCTTAAATCAACCTGACCGTGGCTGATGCGTCCGGTAAGGGTGTTTTTGATCGCGTTCGGCTGTTCGTAGATCTCTTTCTGCATGTAGTGACGGTAAATGCCTTTATCGCCCGCGTCATATTGCAGATTGGATTCGATATCCTGACGTTTTACTTCCGCGCCAGTTTTATCGAAGATGTTTACCGAACGGCGAGTGATTTCCGCAATATCGCCCTCTTCAAGGAAGATAAAGCGACGGGTCACCGGCAACAGCGCCAGCTGGTCAGAAGCGATAAAGTTTTCGCCCATCCCCAGGCCAATCACCAGCGGACTACCAGAACGTGCCGCCAGCAGGGTATCCGGGTGACGGGAGTCCATGATCACTGTACCGTACGCACCACGCAGCTGCGGGATAGCACGCAGAACGGCCTCACGCAGAGTCCCGCCTTGTTTCAGCTCCCAGTTCACCAGATGGGCAATCACTT  >  NC_000913/3912808‑3913292

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: