Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,319 C→T V150V (GTG→GTA pyrE ← orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,815,3190CT98.4% 212.5 / NA 63V150V (GTG→GTApyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (32/30);  minor base G (0/1);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAG  >  NC_000913/3815075‑3815560
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
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)
                          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)
                                                          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)
                                                                                                                          agataggCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTgccgcc                                                                                                                   <  1:111541/245‑1 (MQ=255)
                                                                                                                           gataggcaatcagagttttgaGGGGGATGATAGGGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTGCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCt                                                                                                                  <  2:197223/231‑1 (MQ=255)
                                                                                                                            ataggCAATCAGATCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCtt                                                                                                                 <  2:346154/247‑1 (MQ=255)
                                                                                                                                gCAATCAGATCTTGCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCAcc                                                                                                             <  1:384498/251‑1 (MQ=255)
                                                                                                                                     cagatCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGCCCTTGTAACGCGCCACCAACCAGATTTCCGCCTTCCCCGTggg                                                                                                        >  2:240550/5‑250 (MQ=255)
                                                                                                                                      agatCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCCTCTACCAGCCTTACGCGTCCTTGTAACGCGCTACCCACCAGATTGCCGCCTTCACCGTGGTc                                                                                                       >  2:243618/4‑251 (MQ=255)
                                                                                                                                        atCTTTCAGGGTGATGATAGAGATCCCTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCt                                                                                                                                                <  2:241270/207‑1 (MQ=255)
                                                                                                                                        atCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCt                                                                                                                                                >  1:241270/2‑208 (MQ=255)
                                                                                                                                             tCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATTCCATCTCCTAACAGCATTACGTGTCCTTTTCACCCGCTACCACNCAGACTGCCACCCTCACCGTGGTCttttccc                                                                                                >  2:223825/1‑248 (MQ=255)
                                                                                                                                              cAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCCTCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGccc                                                                                               >  2:291027/1‑249 (MQ=255)
                                                                                                                                               aGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGcctc                                                                                              <  1:38870/251‑4 (MQ=255)
                                                                                                                                                 ggTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGcctctt                                                                                            <  1:329988/251‑1 (MQ=255)
                                                                                                                                                 ggTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGcctctt                                                                                            >  1:146882/1‑251 (MQ=255)
                                                                                                                                                   tgatgaTAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAAtct                                                                                                                                                                                                              >  1:254776/1‑135 (MQ=255)
                                                                                                                                                   tgatgaTAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAAtct                                                                                                                                                                                                              <  2:254776/135‑1 (MQ=255)
                                                                                                                                                     atgatAGAGATCACTTTGCAGTTGTAATCACGCTCACCTTCCTGAATCGCCGACATCTCGCCGCGCCCCAGTTCCTGACGATCAAACGAAATCAATACGCCAGCAAGCGTCGCGCCATCAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGAGGGGATCAAATCATCTACCAGCATTACGCGTCCCTTTAACGACCTCCCAACCAGATCGACCCCGTCACCGCGGACTTTTGcctatttccg                                                                                        >  2:266748/1‑242 (MQ=255)
                                                                                                                                                     atgatAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAACTCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCCCATCATCTCCAAGCATTAAGCGTCCTTGTAACGCGCTACCCACCAGATTGCCGCCTTCACCCGGGTCTTTTGCCTCTTTGCg                                                                                        >  2:88163/1‑251 (MQ=255)
                                                                                                                                                         tAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGACATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGACTCGCCGTTCCGGCGGTGATTACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTCCCCACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGttc                                                                                    >  2:88290/1‑250 (MQ=255)
                                                                                                                                                              aTCACTTTGCATGTGTAATCACGCTCAACGTCCTTAATCGCCGAAATCTCGCCTCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCTCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCCTTCCGGCGGTGATCACATCATCTACCTGCAGTACGCGTCCTTGTACCGCGCTGCCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCGCTTTTCGGTTGAAGCa                                                                               <  2:242623/251‑1 (MQ=255)
                                                                                                                                                                 aCTTTGCAGTTGTAATCACGCTCATCTTCCTGAATCGCCGAAAGTCCGCCGCGCCCCGGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAATCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGGATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTGTTGCCTCGTTGCGGTTAAAGCAGTa                                                                            <  2:84347/251‑1 (MQ=255)
                                                                                                                                                                     tGCAGTTGTAATCACGCTCAACTCCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCTGTTAAAGCAGTACgg                                                                         <  2:265619/250‑1 (MQ=255)
                                                                                                                                                                             tAATCACGCTCAACTTCATGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTTCCGCCTTCACCGTGGTCTTTTTCCCCTTTGCGGTTAAAGCAGTACGGcagggcc                                                                  >  2:367969/1‑246 (MQ=255)
                                                                                                                                                                               atcccGCTTAACTTCCTTAATCGCCGAAATCTCGCCGCGCCCCCGGTCCTGACTATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGTCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCTGTTAAAGCAGTACGGcaggtcaaggt                                                              <  1:35227/247‑1 (MQ=255)
                                                                                                                                                                                    ccTCAACTTCCTGAACCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGa                                                         <  1:423237/250‑1 (MQ=255)
                                                                                                                                                                                            tCCTGAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGGCAAGGGCATGATGCTcccc                                                 >  1:422325/1‑249 (MQ=255)
                                                                                                                                                                                                gAATCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAgc                                             <  2:146882/251‑2 (MQ=255)
                                                                                                                                                                                                 aaTCGCCGAAATCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGAGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGcagggccagggtatggtggcccgcccccg                                            >  2:375105/1‑232 (MQ=255)
                                                                                                                                                                                                  aTCGCCGAAATCTCCCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAgcgc                                           <  2:191190/251‑4 (MQ=255)
                                                                                                                                                                                                           aTCTCGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTGCTGGTGCTCCggccgcgccaccggggg                                  >  1:382223/1‑235 (MQ=255)
                                                                                                                                                                                                              ctcGCCGCGCCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGGCAAGATGGTCCGCCAGCGccccggcggggg                                >  2:327145/1‑241 (MQ=255)
                                                                                                                                                                                                                     cgcCCCCGTTCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAgtggtcgc                         <  1:9585/251‑4 (MQ=255)
                                                                                                                                                                                                                        ccccgTTCCGGACGATCAAGCGAAACCAATCCGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCAAATCGCCGTTCCGGCGGTGAGCTCATCATCTCCCAGCATTTCGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCCCCGTGGTCTGTTGCCTCTTGGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAgtggtcgcgac                      <  2:37360/248‑7 (MQ=255)
                                                                                                                                                                                                                              tCCTGACGATCAAGCGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGa                 <  2:43009/250‑1 (MQ=255)
                                                                                                                                                                                                                                            cGAAATCAATACGCCAGCAAGCGTCGCGCCATTAGCCTGAATAATCTCCATCGATTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCCTCTTTGCGGTTAAAGCAGTACGGCAGGTCAAGGTCATGATGCTCCGCCAGCGCCACGGCGGTAGTGGTCGCAATCGGGATCCCTTTGTAAGCAg  <  1:291027/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
GGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTAATCACGCTCAACTTCCTGAATCGCCGAAATCTCGCCGCGCCCGCGTTCCTGACGATCGAGCGAAATCAACACGCCAGCAAGCGTCGCGCCATTGGCCTGAATAATCTCCATCGACTCGCGAATCGCCGTTCCGGCGGTGATCACATCATCTACCAGCATTACGCGTCCTTGTAACGCGCTACCAACCAGATTGCCGCCTTCACCGTGGTCTTTTGCTTCTTTGCGGTTAAAGCAGTACGGCAGGTCCAGGTCGTGATGCTCCGCCAGTGCCACAGCGGTTGTGGTGGCAATCGGGATCCCTTTGTAAGCAG  >  NC_000913/3815075‑3815560

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