Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 3,815,808 (A)8→7 intergenic (‑40/+26) 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,8010A.96.7% 687.5 / 17.6 152intergenic (‑33/+33)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base . (64/82);  total (66/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 160 total aligned reads displayed.

CAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCC  >  NC_000913_3_bme_pgi/3815652‑3815946
                                                                                                                                                     |                                                                                                                                                 
cAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACa                                                                                                                                                   <  3:128385/150‑1 (MQ=255)
 aGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaa                                                                                                                                                  >  1:598434/1‑150 (MQ=255)
      ggCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaaaa                                                                                                                                             <  3:434716/150‑1 (MQ=255)
            gaagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaaaaggc                                                                                                                                          >  2:1158763/1‑144 (MQ=255)
            gaagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaaaaggc                                                                                                                                          <  1:1158763/147‑4 (MQ=255)
            gaagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGAc                                                                                                                                      >  4:344824/1‑150 (MQ=255)
            gaagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTGGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGAc                                                                                                                                      >  2:1589490/1‑150 (MQ=255)
             aagaaATAGGGGCTTTTGCGCCCGGNTTTCAGCGTAAANTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACt                                                                                                                                     >  1:1392642/1‑150 (MQ=255)
             aagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACt                                                                                                                                     <  3:128277/150‑1 (MQ=255)
              agaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGACTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACtc                                                                                                                                    >  1:436413/1‑150 (MQ=255)
                aaaGAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACtcat                                                                                                                                  <  3:553705/150‑1 (MQ=255)
                    aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGt                                                                                                                              <  1:1191521/150‑1 (MQ=255)
                     ggggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTc                                                                                                                             >  1:1619413/1‑150 (MQ=255)
                      gggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCg                                                                                                                            >  1:21118/1‑150 (MQ=255)
                      gggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCg                                                                                                                            >  2:836922/1‑150 (MQ=255)
                          ttttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCtt                                                                                                                        <  3:570316/150‑1 (MQ=255)
                           tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTa                                                                                                                       >  2:197772/1‑150 (MQ=255)
                            ttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTGATGAGTCGCCtttt                                                                                                                      <  1:864425/150‑3 (MQ=255)
                            ttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTaa                                                                                                                      >  3:20185/1‑150 (MQ=255)
                            ttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTaa                                                                                                                      <  3:511441/150‑1 (MQ=255)
                              gcgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCACCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTaaaa                                                                                                                    <  1:578803/150‑1 (MQ=255)
                                  ccGGATTTCAGCGTAAAGTCGCGAAACTTTAACCCCTGGTTGCGAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTA‑TTAAAAAAAAAAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCa                                                                                                                <  2:240960/150‑1 (MQ=255)
                                  ccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCa                                                                                                                >  3:448679/1‑150 (MQ=255)
                                      aTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGttt                                                                                                            <  1:878168/150‑1 (MQ=255)
                                        ttCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGc                                                                                                          >  2:819302/1‑150 (MQ=255)
                                        ttCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGc                                                                                                          <  3:716258/150‑1 (MQ=255)
                                         tCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTGATGAGTCGCCttttt                                                                                                                     >  2:1267888/1‑135 (MQ=255)
                                          cAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCa                                                                                                        <  2:1619413/150‑1 (MQ=255)
                                            gCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGCTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCCTCAGTCGAATTCAAACTCAGTTTCCCAgc                                                                                                      >  4:874786/1‑150 (MQ=255)
                                             cGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTGATGAGTCGCCttttttt                                                                                                                   >  1:836938/1‑131 (MQ=255)
                                             cGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgcg                                                                                                     >  2:1421887/1‑150 (MQ=255)
                                             cGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgcg                                                                                                     >  3:18932/1‑150 (MQ=255)
                                                  aCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCgccgcc                                                                                                <  4:243048/150‑1 (MQ=255)
                                                     cGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCATTTTGCCAGCGCCGCCTTc                                                                                             >  3:790078/1‑150 (MQ=255)
                                                     cGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTc                                                                                             <  2:549518/150‑1 (MQ=255)
                                                      gCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGtt                                                                                                             >  1:1058151/1‑133 (MQ=255)
                                                            cTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCg                                                                                      <  2:849259/150‑1 (MQ=255)
                                                                aaCACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAc                                                                                  >  1:145104/1‑150 (MQ=255)
                                                                     cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATgg                                                                             >  3:62843/1‑150 (MQ=255)
                                                                      tgcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGga                                                                            <  4:20185/150‑1 (MQ=255)
                                                                        cttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgatt                                                                          <  1:1596022/150‑1 (MQ=255)
                                                                        cttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgatt                                                                          <  1:1024692/150‑1 (MQ=255)
                                                                        cttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgatt                                                                          >  4:696248/1‑150 (MQ=255)
                                                                         ttgctAAGCGCAAATAAAATCAACTGGCGCAAATATGATTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc                                                                         <  4:42329/150‑1 (MQ=255)
                                                                          tgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcg                                                                        >  1:1284182/1‑150 (MQ=255)
                                                                           gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcga                                                                       <  1:210269/150‑1 (MQ=255)
                                                                           gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc                                                                         <  1:581657/148‑1 (MQ=255)
                                                                           gctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc                                                                         >  2:581657/1‑148 (MQ=255)
                                                                            ctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcga                                                                       >  4:118836/1‑150 (MQ=255)
                                                                            ctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgat                                                                      >  2:22945/1‑150 (MQ=255)
                                                                            ctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgat                                                                      >  4:654381/1‑150 (MQ=255)
                                                                              aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgattc                                                                    <  3:773305/150‑1 (MQ=255)
                                                                               aGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTcc                                                                   <  1:403852/150‑1 (MQ=255)
                                                                                gcgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTccc                                                                  >  4:431691/1‑150 (MQ=255)
                                                                                 cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTcccc                                                                 >  1:1523140/1‑150 (MQ=255)
                                                                                 cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTcccc                                                                 >  4:830216/1‑150 (MQ=255)
                                                                                  gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTc                                                                                             >  1:676006/1‑121 (MQ=255)
                                                                                  gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTc                                                                                             <  2:676006/121‑1 (MQ=255)
                                                                                       ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGc                                                           <  4:448679/150‑1 (MQ=255)
                                                                                         cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCa                                                         <  1:1263305/150‑1 (MQ=255)
                                                                                           aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAANAt                                                                                                                  >  3:74750/1‑91 (MQ=255)
                                                                                           aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCGCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCaga                                                       >  1:1319700/1‑150 (MQ=255)
                                                                                             aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGc                                                     <  2:1523140/150‑1 (MQ=255)
                                                                                             aaaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGc                                                     <  1:1454480/150‑1 (MQ=255)
                                                                                              aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTCCTTTTCTACAGAC‑AAAAAAAGGCGTCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGcc                                                    <  2:1272830/150‑1 (MQ=255)
                                                                                                   gcgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCATTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaacaa                                               <  3:278033/150‑1 (MQ=255)
                                                                                                     gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGa                                             <  2:13310/150‑1 (MQ=255)
                                                                                                     gcTGATATGGTTTCATGCCTTCGCTCATCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGa                                             <  1:22945/150‑1 (MQ=255)
                                                                                                      cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGAt                                            <  3:830216/150‑1 (MQ=255)
                                                                                                      cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGAt                                            <  2:733842/150‑1 (MQ=255)
                                                                                                        gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGa                                          <  2:598434/150‑1 (MQ=255)
                                                                                                         atatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCg                                                                                         <  4:621232/102‑1 (MQ=255)
                                                                                                          tatGGTTTCATGCCTTCGCTCCTCATATTACTGTACTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCCGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGt                                        <  4:799442/150‑1 (MQ=255)
                                                                                                            tGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCtt                                                                                              <  2:798939/94‑1 (MQ=255)
                                                                                                            tGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCtt                                                                                              >  1:798939/1‑94 (MQ=255)
                                                                                                             ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGa                                                                                                                                       <  3:686351/52‑1 (MQ=255)
                                                                                                             ggTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGa                                                                                                                                       >  4:686351/1‑52 (MQ=255)
                                                                                                              gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGct                                    >  3:481683/1‑150 (MQ=255)
                                                                                                              gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGct                                    >  3:39169/1‑150 (MQ=255)
                                                                                                                  cATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAg                                         <  4:255783/141‑1 (MQ=255)
                                                                                                                   aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCg                                                                                      >  2:656414/1‑95 (MQ=255)
                                                                                                                   aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCa                               >  2:490596/1‑150 (MQ=255)
                                                                                                                   aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCa                               <  4:18932/150‑1 (MQ=255)
                                                                                                                      ccTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATgg                            >  2:268599/1‑150 (MQ=255)
                                                                                                                      ccTCCGCTCCTCATCTTCCTTTTCTCCAGAC‑AAAAAAAGGCGACTCATCAGCCTCCTTAAAAATCCGTTTGCCAGCGCCGCCTTCTGCGTCGCTACGATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTGTTCATgg                            <  2:611313/150‑1 (MQ=255)
                                                                                                                        tgcGCTCCTCATCTTACTTTTCTACCG‑CCAAAAAAAGGCGCCTCATCAGTCGCCTTAAAAATCCGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGt                          <  2:939254/148‑1 (MQ=255)
                                                                                                                           gCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGa                        >  1:1297806/1‑150 (MQ=255)
                                                                                                                               ctcATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGc                                                                                          >  2:916661/1‑79 (MQ=255)
                                                                                                                               ctcATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGc                                                                                          <  1:916661/79‑1 (MQ=255)
                                                                                                                                tcATCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCGGCGTAGTTGCTATTGAATTGATTGCCCTCGTGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt                  >  1:156147/1‑150 (MQ=255)
                                                                                                                                  aTCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCg                <  1:710694/150‑1 (MQ=255)
                                                                                                                                   tCTTACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGc               <  2:782995/150‑1 (MQ=255)
                                                                                                                                     ttACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCACCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCt             >  1:169155/1‑150 (MQ=255)
                                                                                                                                      tACTTTTCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTc                                <  1:783781/130‑1 (MQ=255)
                                                                                                                                         ttttCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg         >  1:1326587/1‑150 (MQ=255)
                                                                                                                                         ttttCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg         >  2:1539047/1‑150 (MQ=255)
                                                                                                                                            tCTACAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAg                                                                                                       >  1:920331/1‑53 (MQ=255)
                                                                                                                                               aCAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTATAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcc   <  2:1079543/150‑1 (MQ=255)
                                                                                                                                                cAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcc   <  2:1297806/150‑1 (MQ=255)
                                                                                                                                                cAGAC‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc  >  2:119074/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCC  >  NC_000913_3_bme_pgi/3815652‑3815946

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: ‑