Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,418,862 C→T G255D (GGT→GAT)  ulaG ← L‑ascorbate‑6‑phosphate lactonase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,418,8620CT95.5% 190.9 / ‑6.6 66G255D (GGT→GAT) ulaGL‑ascorbate‑6‑phosphate lactonase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/29);  minor base G (0/3);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGATAGTGGTACTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATACCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAGCGCATGAATTTCAATGTCTT  >  NC_000913/4418621‑4419108
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      
cgGATAGTGGTACTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGt                                                                                                                                                                                                                                               >  1:300623/1‑251 (MQ=255)
           aCTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCCTATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTAc                                                                                                                                                                                                                                    >  2:39534/1‑251 (MQ=255)
                aagatGTCTTTATCCAGCGGCCAGGGAAATTTGCTACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCCGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATcc                                                                                                                                                                                                                               >  2:360348/5‑251 (MQ=255)
                    tGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGCACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAAc                                                                                                                                                                                                                           >  1:78325/1‑251 (MQ=255)
                      tCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTAAGATCCTAACGc                                                                                                                                                                                                                         >  2:111622/1‑251 (MQ=255)
                       cTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGCAGTTTGACCAGATATCGTGGTGGAACGGGATCCCTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGTTCCTAACGcc                                                                                                                                                                                                                        >  2:54326/1‑251 (MQ=255)
                                gCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACATCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGGCCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCCTATCGGCGCTGGTCATGTTGTCGGTGATATCGCGCGGATTGTCGCCGTACGATCCTAACGCCACGTCGATc                                                                                                                                                                                                               <  2:96414/251‑1 (MQ=255)
                                  ggCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGACGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTACCGCCCCGTCGATCTg                                                                                                                                                                                                             >  2:96415/1‑251 (MQ=255)
                                    ccAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGAt                                                                                                                                                                                                           <  1:127711/251‑1 (MQ=255)
                                                 ccacccccCTGCAAGATAAACGGCTTGAAGCCATACTTAATGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGCAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATAGGCGCTGGGCACTTTGACGGGGATATCGCGCGGGGTTTCGACGTACGATCCTAACGCCCAGTCGATCTGCTGTTCGTTACCg                                                                                                                                                                                                >  2:372131/1‑249 (MQ=255)
                                                  cacccaccTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCGCCCACAGCACGCGGATCTCTTGCGGATCGGTCTGGAAGTTGGACCAGATATCGTGGTGTAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGc                                                                                                                                                                                             <  1:62641/251‑1 (MQ=255)
                                                   acccaccTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCt                                                                                                                                                                                            <  1:354242/251‑1 (MQ=255)
                                                   acccaccTGACAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTTGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCTCCCATACGCAGCATATCGGCGCTGGTCATTTTGTCTGTGATATCGCGCGGGTTTTCGCCTTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCt                                                                                                                                                                                            <  1:332723/251‑1 (MQ=255)
                                                      caccTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCAAAGCACGCGGCTCTCTTGCGGATCGGCCTGGAAGTTTGACCCGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTAATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTc                                                                                                                                                                                          >  2:190462/1‑250 (MQ=255)
                                                         cTGCCAGATAAACGGCTTGACGCCATAATTCAGGCGCTCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGCTCGGACTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGTTGTTCGTTACCGTGCTTCGCat                                                                                                                                                                                      >  2:46555/1‑251 (MQ=255)
                                                           gCCAGATAATCGGCTTGAAGCCATTCTTCAGGCGAGCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATTCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCataa                                                                                                                                                                                    <  2:290951/251‑1 (MQ=255)
                                                              aGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGACGATCTGATGTTCGTTACCGTGCTTCGCataat                                                                                                                                                                                   >  1:446404/1‑249 (MQ=255)
                                                                          tGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGCTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTgg                                                                                                                                                                     >  1:165081/1‑251 (MQ=255)
                                                                            aaGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTggg                                                                                                                                                                    >  1:195759/1‑250 (MQ=255)
                                                                               ccATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTc                                                                                                                                                                <  1:39534/251‑1 (MQ=255)
                                                                                  tACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGACCGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATAACGGCGCTGGTCATTTTGTCGGTTATATCACGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGGTAGAGTAGTGGGAGTcgcc                                                                                                                                                             >  2:26357/1‑251 (MQ=255)
                                                                                   acttcctgctatctttgtttagtacccACAGCTGGCGGATATCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGGGTTTTAACGGGAGCACTACTTTCGCATTCAGCGCTTCACCCATACGCATCATAGCGGCGCTTGGCATGTTGTCGGTGATATTGTGCTGGTTTTCGCTGTACGATCCTAACTCCACGTCGTTCTGATGTACGTGACCTTGCTTCGTATAATAGGTATATTAGAGGCAGTcgccg                                                                                                                                                            <  2:263501/227‑1 (MQ=255)
                                                                                       cAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGTCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTGTCGCATTCAGCGCTTCACCCATACGCATCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGGCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCtg                                                                                                                                                         <  1:453798/250‑1 (MQ=255)
                                                                                         ggCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGAGATCGTGGTGGAACGGGATCACTCCTTTCTCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACTTCGATCTGATGGTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGa                                                                                                                                                      <  2:152729/251‑1 (MQ=255)
                                                                                            gagcTTTTATCAGCTCCCACAGCACGCGGATCGCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGTTGGAACGGGATCACTACTGTCGCATTCAGCGCTTCGCCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATAGCGCGCGGGTTGTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATAc                                                                                                                                                   <  1:432203/248‑1 (MQ=255)
                                                                                                    ttCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGCTTGACCAGATATCGTGGCGGAACGGGCTCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAgcc                                                                                                                                           >  2:156211/1‑251 (MQ=255)
                                                                                                     tCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAgccg                                                                                                                                          >  1:15787/1‑251 (MQ=255)
                                                                                                            ccACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTGGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCg                                                                                                                                   <  2:165081/251‑1 (MQ=255)
                                                                                                             cacaGCACGCGGCTCTCTTGCGGATCGGCCTGGCAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGt                                                                                                                                  >  2:322486/1‑251 (MQ=255)
                                                                                                                aGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGtttt                                                                                                                               >  1:455835/1‑251 (MQ=255)
                                                                                                                aGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGATAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGtttt                                                                                                                               >  1:321455/1‑251 (MQ=255)
                                                                                                                 gCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTg                                                                                                                              <  2:20630/251‑1 (MQ=255)
                                                                                                                       ggATCGCTTGCGGATCGGACTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGTAGCCGCCAGGCGTTTTGAACAgg                                                                                                                        <  2:300623/251‑1 (MQ=255)
                                                                                                                                  ggATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCTCCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACcgc                                                                                                             <  2:209369/251‑1 (MQ=255)
                                                                                                                                         ccTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCg                                                                                                      <  1:13116/251‑1 (MQ=255)
                                                                                                                                         ccTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCg                                                                                                      <  1:171627/251‑1 (MQ=255)
                                                                                                                                                    gACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATc                                                                                           >  2:93266/1‑251 (MQ=255)
                                                                                                                                                       cAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTgg                                                                                        >  1:44363/1‑251 (MQ=255)
                                                                                                                                                          atatCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAg                                                                                     <  2:44363/251‑1 (MQ=255)
                                                                                                                                                            accGTGGTGGAACTGGATCTCTACTTGCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTGGTCGGTTATCTCGCTCGTGTTTTCTGCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGTGGCGCATAATAGTTATAGTATGGGGAGTCGCCGCTTTGATACAGGGAGCCGCCAGGCGTTTTGAACATGTAGTTCACCGCGCGATCGCTTATGCCATCTGTCAGTa                                                                                   <  2:208395/249‑1 (MQ=255)
                                                                                                                                                              cGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACg                                                                                 >  1:107236/1‑251 (MQ=255)
                                                                                                                                                                                    tttCGCATTCAACGCTTCACCCATACGCAGCCTATCCGCGCTGGACATTTTGTCGGTGATATCCCGCGGGTTTTCGCCGTACGGTCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGCAGTGGGAGTCGCCGCTGTGATACAGGGAGCAGCCAGGCGATTTTAACAGGTAGTTCACCGCGCGATCGTCCCGGCCATCTGGCAGTACGCCAGCAGCTTTTTGTTCggccg                                                           >  2:386810/1‑249 (MQ=255)
                                                                                                                                                                                       cGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCggcaggc                                                         >  1:164059/1‑250 (MQ=255)
                                                                                                                                                                                          aTTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAggg                                                     >  2:346666/1‑251 (MQ=255)
                                                                                                                                                                                            tCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGt                                                    <  1:337064/250‑1 (MQ=255)
                                                                                                                                                                                             cAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTg                                                   >  1:154160/1‑250 (MQ=255)
                                                                                                                                                                                               gcgcTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTACTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCAGCCAGGCGTTTTGAACAGGGAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTGTTTGATCGGCAGGCAGGGTgctc                                                >  1:456397/1‑248 (MQ=255)
                                                                                                                                                                                                     cccccATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGc                                           <  2:164059/248‑1 (MQ=255)
                                                                                                                                                                                                       cccATACGCAGCATATCGGCGCTGGTCTTTTTGTCGGTGATAGCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTGACCGTGCGTCGCATAATAGTTAGAGTAGGGGGAGTCGCCGCTGTGATACAGGGAGGCGCCAGGCGTTTTGAGCGGGTAGTTCACCGCGCGATCGTCCATGCCAGCTGGCAGTACGCGAGCCGCGTTTTGATCGTCAGGCTGGGTGTTTAGTGCAGt                                        <  2:57001/251‑1 (MQ=255)
                                                                                                                                                                                                       cccATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTCACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGtttt                                                                                                                               >  2:79491/1‑164 (MQ=255)
                                                                                                                                                                                                       cccATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGtttt                                                                                                                               <  1:79491/164‑1 (MQ=255)
                                                                                                                                                                                                       cccATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGt                                        <  1:93266/251‑1 (MQ=255)
                                                                                                                                                                                                         cATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGACTAATAGTTAGGGTAGTGGGAGTCGCCGCTGTGCTACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGCTCggcaggccgagtgagcagggccgaac                                      >  2:408212/1‑232 (MQ=255)
                                                                                                                                                                                                             cGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTc                                  >  1:221647/1‑251 (MQ=255)
                                                                                                                                                                                                                  cagatCGGCGCTGGTCATTTTGTCGTGGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCAGGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAg                             <  2:319750/248‑1 (MQ=255)
                                                                                                                                                                                                                     atCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGACCAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCAt                          >  1:399477/1‑251 (MQ=255)
                                                                                                                                                                                                                     atCGGCGCTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCa                           <  2:195759/250‑1 (MQ=255)
                                                                                                                                                                                                                         gcgcTGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAg                      <  1:75504/251‑1 (MQ=255)
                                                                                                                                                                                                                             tGGTCATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCTTTACCGTGCTTCGCCTAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAgcgc                   >  1:204027/1‑250 (MQ=255)
                                                                                                                                                                                                                                 cATTTTGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAGCGCATGa               >  1:308215/1‑250 (MQ=255)
                                                                                                                                                                                                                                    tttGTCGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAGCGCATGAAtt            <  1:344011/250‑1 (MQ=255)
                                                                                                                                                                                                                                        gcGGTGATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAGCGCATGAATTTCAAt       <  2:107236/250‑1 (MQ=255)
                                                                                                                                                                                                                                         cGGTGATAGCGGGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCGGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAGCGCATGAATTTCAATg      <  2:324935/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gggATATCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGAGCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGAGACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCTTCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTc                                  <  1:64493/219‑1 (MQ=255)
                                                                                                                                                                                                                                            ggATATCGCGCGTGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCAGTGATACCGGGAGCCGCCATGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTGTTTGAGCTGCAGGCAGTGTGATCAGTGCAGTACGGTCGAAAGCAGCAAGCGCATGACTTTCAATGTCt   <  2:178936/250‑1 (MQ=255)
                                                                                                                                                                                                                                             gATATCGCGCGGGTTTTCACCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGCGTAGTGGGAGTCGCCGCTGTAATACAGGGAGCCGCCAGGCGTTTTGAACAGGAAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTCATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGCAAGCATCAAGCGCATGAATTTCAATGTCtt  >  2:221634/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      
CGGATAGTGGTACTCGAAGTTGTCTTTATCCAGCGGCCAGGTAAATTTGCCACCCACCTGCCAGATAAACGGCTTGAAGCCATACTTCAGGCGATCTTTTTTCATCTCCCACAGCACGCGGATCTCTTGCGGATCGGCCTGGAAGTTTGACCAGATATCGTGGTGGAACGGGATCACTACTTTCGCATTCAGCGCTTCACCCATACGCAGCATATCGGCGCTGGTCATTTTGTCGGTGATACCGCGCGGGTTTTCGCCGTACGATCCTAACGCCACGTCGATCTGATGTTCGTTACCGTGCTTCGCATAATAGTTAGAGTAGTGGGAGTCGCCGCTGTGATACAGGGAGCCGCCAGGCGTTTTGAACAGGTAGTTCACCGCGCGATCGTCCATGCCATCTGGCAGTACGCCAGCCGCTTTTTGATCGGCAGGCAGGGTGATCAGTGCAGTACGGTCGAAAGCATCAAGCGCATGAATTTCAATGTCTT  >  NC_000913/4418621‑4419108

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