Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 3,815,825 C→T intergenic (‑57/+9) pyrE ← / ← rph orotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi3,815,8250CT95.8% 456.1 / 5.8 143intergenic (‑57/+9)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base T (61/75);  total (65/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

Only 100 of 143 total aligned reads displayed.

TTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCT  >  NC_000913_3_bme_pgi/3815678‑3815971
                                                                                                                                                    |                                                                                                                                                  
ttttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCt                                                                                                                                                  <  2:456579/150‑1 (MQ=255)
 tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTATCAACAAAAAAAAAAAAGGCGACTCATCAGTCGTCtt                                                                                                                                                 <  2:1252750/150‑1 (MQ=255)
  ttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTa                                                                                                                                                <  3:759518/150‑1 (MQ=255)
     cgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTaaaa                                                                                                                                             <  1:898826/150‑1 (MQ=255)
       cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCtttttttt                                                                                                                                           <  3:654847/150‑7 (MQ=255)
       cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAAt                                                                                                                                           <  1:1581025/150‑1 (MQ=255)
         cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCa                                                                                                                                         <  3:391164/150‑1 (MQ=255)
             tttCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGttt                                                                                                                                     <  4:644895/150‑1 (MQ=255)
                  gCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGTCGACTCATCAGTCGTCTTAAAATTCAGTTTGcccg                                                                                                                                >  3:202434/1‑148 (MQ=255)
                   cGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAgc                                                                                                                               >  4:938079/1‑150 (MQ=255)
                      aaaCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCCCATCTTACTTTTCTACAGACAAAAAAAAGGCGACCCATCAGTCGTCTTAAAAATCAGTTTGCCAGcgcc                                                                                                                            >  1:1020786/1‑150 (MQ=255)
                       aaCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGcgccg                                                                                                                           >  4:834002/1‑150 (MQ=255)
                        acgcGCCAAACTTTAACACCTGCTTGCGAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTATAAAAAAAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGcgccgc                                                                                                                          <  2:715523/147‑1 (MQ=255)
                        aCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCCTCTTAAAAATCAGTTTGCCAGcgccgc                                                                                                                          >  2:437907/1‑150 (MQ=255)
                         cTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCgccgcc                                                                                                                         >  3:824458/1‑150 (MQ=255)
                             ccAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTCCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACACATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCAGCCTTCt                                                                                                                     >  1:1026631/1‑150 (MQ=255)
                             ccAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCt                                                                                                                     >  2:502604/1‑150 (MQ=255)
                             ccAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCt                                                                                                                     <  4:583749/150‑1 (MQ=255)
                             ccAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCt                                                                                                                     <  3:834002/150‑1 (MQ=255)
                              ccaaCTGTAGCACCTGCTTGCTAAGCGCAAATTCAATAACCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAAACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTg                                                                                                                    <  2:154497/148‑1 (MQ=255)
                               aaaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGc                                                                                                                   <  2:1118170/150‑1 (MQ=255)
                                   tttAACACCTGCTTGCTAAGCGCAAGTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCg                                                                                                               <  2:354207/150‑1 (MQ=255)
                                    ttAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGc                                                                                                              <  2:693922/150‑1 (MQ=255)
                                     tAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCt                                                                                                             >  2:1437026/1‑150 (MQ=255)
                                      aaCACCTGCTTGCTAAGAGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTATACAGAAAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTa                                                                                                            <  2:214991/150‑1 (MQ=255)
                                       acacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAc                                                                                                           >  3:145226/1‑150 (MQ=255)
                                       acacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAc                                                                                                           >  3:1036548/1‑150 (MQ=255)
                                         cccTGCTTGCTACGCACAAAGCCAAGAAACGTACGCTAACATGTTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACaa                                                                                                         <  1:517412/149‑1 (MQ=255)
                                         acCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACaa                                                                                                         >  2:1528749/1‑150 (MQ=255)
                                         acCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACaa                                                                                                         <  1:191522/150‑1 (MQ=255)
                                                 gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAgc                                                                                                                               <  2:1033467/120‑1 (MQ=255)
                                                   tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgatt                                                                                                   <  1:683844/146‑1 (MQ=255)
                                                   tAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgatt                                                                                                   >  2:683844/1‑146 (MQ=255)
                                                    aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgatt                                                                                              <  3:792318/150‑1 (MQ=255)
                                                     aGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgattc                                                                                             <  1:725201/150‑1 (MQ=255)
                                                          aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTc                                                                                        >  2:1651425/1‑150 (MQ=255)
                                                            aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCgg                                                                                      <  4:385472/150‑1 (MQ=255)
                                                              tCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGc                                                                                    <  1:298864/150‑1 (MQ=255)
                                                                aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc                                                                                                  <  1:1756270/134‑1 (MQ=255)
                                                                aaTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCa                                                                                  >  3:417658/1‑150 (MQ=255)
                                                                 aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCag                                                                                 >  2:376114/1‑150 (MQ=255)
                                                                          cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaacaa                                                                        <  1:1586032/150‑1 (MQ=255)
                                                                            cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGa                                                                      <  1:1702457/150‑1 (MQ=255)
                                                                             tGATATGGTTTCATGACTTGTCGCGTCCTTGTGCTTTTCTACATACAAAAAAAAGGCGGCTCAGCAGTCGTCTTAAAAGTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGAt                                                                     <  4:1031477/150‑1 (MQ=255)
                                                                                 atGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTc                                                                                        >  3:596463/1‑127 (MQ=255)
                                                                                  tGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTa                                                                >  3:1068564/1‑150 (MQ=255)
                                                                                   ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                                               <  4:804766/150‑1 (MQ=255)
                                                                                         aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgatt                                                                                              >  2:1314390/1‑113 (MQ=255)
                                                                                         aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTc                                                         >  3:578946/1‑150 (MQ=255)
                                                                                         aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTc                                                         <  4:145226/150‑1 (MQ=255)
                                                                                           gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc                                                                                                  <  1:1004386/108‑1 (MQ=255)
                                                                                           gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCa                                                        <  2:382367/150‑1 (MQ=255)
                                                                                            ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATg                                                      >  2:31038/1‑150 (MQ=255)
                                                                                            ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATg                                                      >  2:338358/1‑150 (MQ=255)
                                                                                                       cATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                                           <  1:674693/150‑1 (MQ=255)
                                                                                                         tCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGa                                                                   <  2:610343/124‑1 (MQ=255)
                                                                                                         tCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGa                                                                   >  1:610343/1‑124 (MQ=255)
                                                                                                         tCTTACGTTTCTACGGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGTCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCg                                         <  4:436558/150‑1 (MQ=255)
                                                                                                          cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCaga                                                                                >  1:1043538/1‑110 (MQ=255)
                                                                                                          cTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGc                                        >  3:78942/1‑150 (MQ=255)
                                                                                                           ttactttactACAGACAAAAAAAAGGCAACACATCAGGCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGTGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGcc                                       <  1:499602/150‑1 (MQ=255)
                                                                                                           ttACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgatt                                                                                                   >  1:780020/1‑90 (MQ=255)
                                                                                                           ttACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgatt                                                                                                   <  2:780020/90‑1 (MQ=255)
                                                                                                               ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTccc                                                                                           >  4:346521/1‑94 (MQ=255)
                                                                                                               ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTccc                                                                                           <  3:346521/94‑1 (MQ=255)
                                                                                                               ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCt                                   >  1:494322/1‑150 (MQ=255)
                                                                                                               ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCt                                   <  4:1068564/150‑1 (MQ=255)
                                                                                                                  tCTACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc                                <  4:578946/150‑1 (MQ=255)
                                                                                                                    tACAGACAAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCtt                                     <  2:1710394/144‑1 (MQ=255)
                                                                                                                    tACAGACAAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCtt                                     >  1:1710394/1‑144 (MQ=255)
                                                                                                                    tACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTc                                                         <  4:244108/123‑1 (MQ=255)
                                                                                                                    tACAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTc                                                         >  3:244108/1‑123 (MQ=255)
                                                                                                                     aCAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGc                                                                                    <  2:146376/95‑1 (MQ=255)
                                                                                                                     aCAGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGc                                                                                    >  1:146376/1‑95 (MQ=255)
                                                                                                                       aGACAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc                           <  1:1651425/150‑1 (MQ=255)
                                                                                                                         aCAAAAACAAGGTGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCt                         <  4:1098623/150‑1 (MQ=255)
                                                                                                                         aCAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcga                                                                                                >  1:1246504/1‑79 (MQ=255)
                                                                                                                         aCAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc                                                                                                  <  3:973715/77‑1 (MQ=255)
                                                                                                                         aCAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAAc                                               <  2:1692359/128‑1 (MQ=255)
                                                                                                                         aCAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCt                         >  2:480927/1‑150 (MQ=255)
                                                                                                                          cAAAAAAAATGCGACTCATCAGTCGCCTTAAAACTCAGCTTGCCAGTGCCGCCTTCTGCGTCGCTACAATGGATTCGATTGCCCTCGGACCAGAGCCAACGGGATGTGTAGCTCTTCTTGGGTGAACGGCTCGCCTTGTGCCGTCCCCTg                        <  2:1026631/150‑1 (MQ=255)
                                                                                                                          cAAAAAAAAGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                        >  2:562513/1‑150 (MQ=255)
                                                                                                                            aaaaaaaaGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGc                       >  1:1230909/1‑150 (MQ=255)
                                                                                                                             aaaaaaaGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCa                      <  4:78942/150‑1 (MQ=255)
                                                                                                                             aaaaaaaGGCGACTCATCAGTCGTCTTAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCAc                     >  2:824545/1‑150 (MQ=255)
                                                                                                                                aaaaGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACtt                   <  2:6175/150‑1 (MQ=255)
                                                                                                                                   aGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCaa                <  1:31038/150‑1 (MQ=255)
                                                                                                                                   aGGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCaa                <  1:1736245/150‑1 (MQ=255)
                                                                                                                                    ggCGACTCATCAGTCGTCTTAAAAATGAGTTTGCCAGCGCCGCCTTCTGCGTCGGTACAATGGATTCGATTCCCCTCGGGCCAGAGCCGACAAGATGAGTTGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAat               <  2:945193/150‑1 (MQ=255)
                                                                                                                                    ggCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACACTGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAat               <  1:1351334/150‑1 (MQ=255)
                                                                                                                                     gCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGc                                 >  2:117424/1‑131 (MQ=255)
                                                                                                                                     gCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGc                                 <  1:117424/131‑1 (MQ=255)
                                                                                                                                     gCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatg              <  1:480927/150‑1 (MQ=255)
                                                                                                                                     gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatg              >  2:1347916/1‑150 (MQ=255)
                                                                                                                                       gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgat            <  1:1347916/150‑1 (MQ=255)
                                                                                                                                       gACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAatgat            >  4:216382/1‑150 (MQ=255)
                                                                                                                                          tcatcaGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATgcg         >  3:309182/1‑150 (MQ=255)
                                                                                                                                            atcaGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCa                                                        <  3:694126/101‑1 (MQ=255)
                                                                                                                                            atcaGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCa                                                        >  4:694126/1‑101 (MQ=255)
                                                                                                                                                 tCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCt  >  1:1664004/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCT  >  NC_000913_3_bme_pgi/3815678‑3815971

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