Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,202 A→G A189A (GCT→GCC pyrE ← orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,815,2020AG98.3% 214.1 / ‑6.9 59A189A (GCT→GCCpyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (30/27);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACC  >  NC_000913/3814956‑3815436
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          
gAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCaa                                                                                                                                                                                                                                        <  1:339230/251‑1 (MQ=255)
   ccccTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGCCTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCa                                                                                                                                                                                                                                     <  1:170626/250‑1 (MQ=255)
   aCCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATc                                                                                                                                                                                                                                      >  1:370668/1‑250 (MQ=255)
   aCCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATc                                                                                                                                                                                                                                      >  1:346322/1‑250 (MQ=255)
               acccatacccttgcttttatgccaGGTTTGCTCTCATGTTTTGTCCCCAGCGAATTTAACTCCCACCCGTAGGAGTATTTGCGCGCCCGCTGGCCGCTAATCGCGGCACATTTTAAGTATTATGCCGCCGGATAATTTTTCCTCCGCCGAGTTTCTTTAAACGCCAAACTCTTCGCGATATGCTTTAACCGCCGCCTGATGTTCCGCCATTTCCGGCTTCTCTGCCAGATAGGCAATCAGATCTTTCAggg                                                                                                                                                                                                                         <  2:265147/228‑1 (MQ=255)
                     acacTGGCCTGCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATGGCGGCACAGTTGCAGAAAGATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTGTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGtgatga                                                                                                                                                                                                                   <  2:370668/248‑1 (MQ=255)
                          ggCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATag                                                                                                                                                                                                                >  1:265619/1‑249 (MQ=255)
                            ccTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATag                                                                                                                                                                                                                >  1:72601/1‑247 (MQ=255)
                             cTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATag                                                                                                                                                                                                                >  1:37360/1‑246 (MQ=255)
                                    gaaggtgtggtgtaagtttttGTCCGCAGCGAATTTAAATACCGCCGGACGGATGATTTGGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAACCTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGc                                                                                                                                                                                                    <  2:339171/234‑1 (MQ=255)
                                    gAAGGTATGCTGTCACTTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGTATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGc                                                                                                                                                                                                    <  1:316506/251‑1 (MQ=255)
                                          atgcgggcACGTTTTGTCCGCAGCGAATTTAAATCCCGCCCGACGGATGATTGTGACGCCCGGTGGTCGGTAAGTGCGGCACAGTTGCAGTAATATGGCGCCGGATGACTTTTCATCCGGGGAGTTTCTTTAAACGCCAAGCTCTTCGCGATAGGCGTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTTTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTg                                                                                                                                                                                               <  1:266859/243‑1 (MQ=255)
                                                       ttGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAAc                                                                                                                                                                                 >  2:335850/1‑251 (MQ=255)
                                                         gcgcgtaGCGATTTTAAGTACGGGCCGACGAATGATTTTGACGCCCGTTGGCAGTTAATTTCGGCAAAGTTGCAGTATTAGGCCACCGGATGTCTTTTCAACCGGCGAGTTTCTTTAATCGCAAATCTCTTCGCTATAAGCCTTAACCGCCGACTTATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAAACAGAGCTTTCAGGGTGATGATAGATATCACTTTGCAGTCGTAATCTCGCTTAACt                                                                                                                                                                                <  2:313417/244‑1 (MQ=255)
                                                         gTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGAAAGGCAATCAGATCTTTCAGGGAGATGATAGANATCACTTTGCAGTTGTAATCACGCCCAACtt                                                                                                                                                                               >  2:203251/1‑251 (MQ=255)
                                                             gCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTg                                                                                                                                                                           >  1:10542/1‑251 (MQ=255)
                                                              cAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTg                                                                                                                                                                           >  1:389549/1‑250 (MQ=255)
                                                                gCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAAt                                                                                                                                                                        <  2:389549/251‑1 (MQ=255)
                                                                     tttAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCATCTAATTGCGGACCAGTTGCAGCACTATGCCGCAGGATGAATTTTCTTCCTGCGAGTTTCTTTAACTGCCAAACTCTTCGCGATATGCCTTAACAGCCGCCAGATGCTCCACCATTTCCGGCTTCTCTTCCCGATAAGCACTCAGTTTTTTCAAGATGATGATAGAAATCACTTTGAAGTTGTAATCACGCTCAACTTCCTGAATCGCCg                                                                                                                                                                   >  2:107728/1‑251 (MQ=255)
                                                                                              ttgccGCCCGCGGGCCGTTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTAAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCt                                                                                                                                           <  2:335408/246‑1 (MQ=255)
                                                                                                 aCGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACg                                                                                                                                       <  1:334445/251‑1 (MQ=255)
                                                                                                  cGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACg                                                                                                                                       >  1:43009/1‑250 (MQ=255)
                                                                                                       gCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATcaag                                                                                                                                 <  1:352322/251‑4 (MQ=255)
                                                                                                        gtgccctctattggCGGCACAGTGGAAGTAAAATAACGCCGGAGAACTTTCCATCCGGCGAGTTTCTTTAATCGCAAAACTCTTCGCGTGAGTCCTTAACCGCCTCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGTTAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATcaag                                                                                                                                 <  2:194984/237‑4 (MQ=255)
                                                                                                            ccGCTAATTGCGGCACAGTTGCATTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGaaa                                                                                                                            <  2:345093/251‑1 (MQ=255)
                                                                                                            ccGCTAATTGCGGCACAGTTGAAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGACGCGCCCCCGTTCCTGACGATCAAGCGaa                                                                                                                             >  2:423298/1‑250 (MQ=255)
                                                                                                                       ggCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGcc                                                                                                                  >  1:432589/1‑250 (MQ=255)
                                                                                                                                     attaggACGCCGGATGACTTTCCAGCCGGCGAGTTTCTTTAAACGCCAATCTCTTCGGGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATATGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCGGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAGGCGTCGCGcc                                                                                                   <  2:422798/246‑1 (MQ=255)
                                                                                                                                     aaTATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGCCGATCAAGCGAAATCAATAAGCCAGCAAGCGTCGCGcc                                                                                                   >  2:151919/1‑251 (MQ=255)
                                                                                                                                          gACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCAtt                                                                                                >  2:285109/1‑249 (MQ=255)
                                                                                                                                                ggATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTcgcc                                                                                                                                                         >  1:184438/1‑186 (MQ=255)
                                                                                                                                                ggATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTAAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTcgcc                                                                                                                                                         <  2:184438/186‑1 (MQ=255)
                                                                                                                                                 ggtGACTTTTCATCCGGGGAGTTTCTTTAAANCCCAAACTCTTCGCGATCGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCTCCATTAGCCTGaat                                                                                       <  2:212996/249‑1 (MQ=255)
                                                                                                                                                  aTGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGaat                                                                                       >  1:346154/1‑250 (MQ=255)
                                                                                                                                                    gACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGaataat                                                                                    <  1:285109/251‑1 (MQ=255)
                                                                                                                                                         ttCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTcc                                                                                <  1:380916/250‑1 (MQ=255)
                                                                                                                                                              ccGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCACGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGat                                                                          >  2:51858/1‑250 (MQ=255)
                                                                                                                                                                ggCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGattc                                                                        >  1:26589/1‑248 (MQ=255)
                                                                                                                                                                  cGAGTTTCTTTAACCGCCAATCTCTTCGCGATAGGCCTTACCCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTGCCAGATAGGCAATCAGATCTTTCAGGGTGACGATAGAGATCACTTTGCAGTTTTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTcgc                                                                      <  2:423434/251‑1 (MQ=255)
                                                                                                                                                                           ttAAACGCCAACCTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACTATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCAGCGATTCGCGAATCg                                                                <  1:353073/248‑1 (MQ=255)
                                                                                                                                                                             aaaCGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACg                                                                                                                                       >  2:253456/1‑175 (MQ=255)
                                                                                                                                                                             aaaCGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACg                                                                                                                                       <  1:253456/175‑1 (MQ=255)
                                                                                                                                                                                 gCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCcgg                                                       >  1:298604/1‑251 (MQ=255)
                                                                                                                                                                                    aaaCTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCcggcgg                                                    >  2:322680/1‑251 (MQ=255)
                                                                                                                                                                                     aaCTCTTCGCGATAGGCCTTAACCGCCGCCCGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCcggcggg                                                   >  1:138404/1‑250 (MQ=255)
                                                                                                                                                                                                  aGGCTTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCt                                      <  2:346322/251‑1 (MQ=255)
                                                                                                                                                                                                     ccTTAACCGCCGCCAGATGTTCCGACATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAAGGAGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGACGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGttctaaactcatcttac                                   >  1:196220/1‑235 (MQ=255)
                                                                                                                                                                                                         aaCCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGc                                >  1:46368/1‑250 (MQ=255)
                                                                                                                                                                                                          aCCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTAAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGTCGATCAAGCGAAATCAATACGCCCGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACCTCTTCTACCGGCAt                              >  1:86476/1‑251 (MQ=255)
                                                                                                                                                                                                             gccgccAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTAc                           <  1:51858/251‑1 (MQ=255)
                                                                                                                                                                                                               cgccAGATGTTCCGCCATGTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTGCCGGCGGTGAGCACATCATCTACCAGCATTAcgc                         <  2:180019/251‑1 (MQ=255)
                                                                                                                                                                                                                   agctGTTCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCTGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTcc                     <  2:243854/248‑1 (MQ=255)
                                                                                                                                                                                                                         tCCGCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTa                <  2:72601/250‑1 (MQ=255)
                                                                                                                                                                                                                            gCCATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCACCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCCTCTACCAGCATTCCGCGTCCTTGTAAcgc            >  2:291383/1‑251 (MQ=255)
                                                                                                                                                                                                                              cATTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTTCGCGTCCTTGTAACgcgc          >  2:46016/1‑251 (MQ=255)
                                                                                                                                                                                                                               aTTTCCGGCTTCTCTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCCCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGAACACATCATCTACCCGCATTACGCGTCCTTGTAACGCGCt         >  1:296541/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ccGGCTTCTCTTCCAGATAGGCAATCAGCTCTTTCAGGGTGATGATAGAGATCACTTCGCAGTTGTAATCACGCTCAACTTCCTGACTCGTCGAAATCTCGCCGCGCCCCCGTTCCTGACGAACAAGCTAAATCAAGACGCCAGCAAGCGGCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAAcgctgtacca     <  2:46884/251‑8 (MQ=255)
                                                                                                                                                                                                                                         tctcTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCCTCTACCAGCATTTCGCGTCCTTGTAACGCGCTacccacc  >  2:282872/1‑248 (MQ=255)
                                                                                                                                                                                                                                         tctcTTCCAGATAGGCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTaccaacc  >  1:199728/1‑248 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          
GAAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACC  >  NC_000913/3814956‑3815436

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