Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 3,815,761 G→T P3Q (CCA→CAA)  pyrE ← orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi3,815,7610GT95.1% 290.1 / 4.5 82P3Q (CCA→CAA) pyrEorotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (35/43);  total (37/45)
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.89e-01

AACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGAT  >  NC_000913_3_bme_pgi/3815615‑3815904
                                                                                                                                                  |                                                                                                                                                
aaCGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGtt                                                                                                                                               <  3:323132/150‑1 (MQ=255)
 aCGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGttt                                                                                                                                              >  1:427704/1‑150 (MQ=255)
  cGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTc                                                                                                                                             <  3:658388/150‑1 (MQ=255)
   ggCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCa                                                                                                                                            <  4:385182/150‑1 (MQ=255)
       tAACAGTGCCAGATCGCGCCCGGTATTCAACAGCCCGGCGTTGAAGAACTAGGGGCTTTTGCGCCCGGATGTCAGGGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAGCTGGCGCTGATATGGTTTCATGcc                                                                                                                                        <  3:270697/150‑1 (MQ=255)
         aCAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCtt                                                                                                                                      <  1:472967/150‑1 (MQ=255)
         aCAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCtt                                                                                                                                      >  4:624224/1‑150 (MQ=255)
              gCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGctc                                                                                                                                 <  2:529038/150‑1 (MQ=255)
               ccAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGctcc                                                                                                                                <  1:77465/150‑1 (MQ=255)
                cAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGctcct                                                                                                                               >  2:919269/1‑150 (MQ=255)
                       cgcCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTa                                                                                                                        >  3:266458/1‑150 (MQ=255)
                        gcCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTAc                                                                                                                       >  4:366571/1‑150 (MQ=255)
                          ccGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACtt                                                                                                                     <  1:465356/150‑1 (MQ=255)
                            ggTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACtttt                                                                                                                   <  1:893850/150‑1 (MQ=255)
                              tATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCt                                                                                                                 <  4:178162/150‑1 (MQ=255)
                              tATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCt                                                                                                                 <  4:178166/150‑1 (MQ=255)
                                 tAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTCCa                                                                                                              <  2:464093/150‑1 (MQ=255)
                                      aGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaa                                                                                                        <  4:533168/150‑1 (MQ=255)
                                       gCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaa                                                                                                       <  4:670086/150‑1 (MQ=255)
                                        cccGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaa                                                                                                      >  2:721359/1‑150 (MQ=255)
                                        cccGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaa                                                                                                      >  2:721379/1‑150 (MQ=255)
                                           ggCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaaaa                                                                                                   <  1:427518/150‑1 (MQ=255)
                                                    cggttcggggCTTTTTGGCCCGGATTTCAGCGTTAACTCGCCAAAATTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCCCTTATATTGTTTCCTGTCTTCGCTCCTCATTTTTCTTTTCTAAAAAAAAAAAAAAAGGCGACt                                                                                           <  2:541594/144‑1 (MQ=255)
                                                       aTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcatc                                                                                       >  1:1028511/1‑150 (MQ=255)
                                                         aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtgatgag                                                                                     <  1:721359/150‑7 (MQ=255)
                                                         aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtgatgag                                                                                     <  1:721379/150‑7 (MQ=255)
                                                           gggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGc                                                                                                                                         >  4:231984/1‑97 (MQ=255)
                                                           gggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGc                                                                                                                                         <  3:231984/97‑1 (MQ=255)
                                                           gggCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTc                                                                                   >  2:927248/1‑150 (MQ=255)
                                                             gCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGCCTCATCAGTCGc                                                                                 >  1:493788/1‑150 (MQ=255)
                                                               ttttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCt                                                                               >  4:657387/1‑150 (MQ=255)
                                                               ttttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCt                                                                               >  1:85507/1‑150 (MQ=255)
                                                                  tGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaa                                                                            >  4:595940/1‑150 (MQ=255)
                                                                    cgcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaa                                                                          <  1:248343/150‑1 (MQ=255)
                                                                       ccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATc                                                                       >  4:174020/1‑150 (MQ=255)
                                                                        cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATg                                                                                                                                                  <  2:846922/75‑1 (MQ=255)
                                                                        cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATg                                                                                                                                                  >  2:599928/1‑75 (MQ=255)
                                                                        cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATg                                                                                                                                                  <  1:599928/75‑1 (MQ=255)
                                                                        cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATg                                                                                                                                                  >  1:846922/1‑75 (MQ=255)
                                                                        cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCa                                                                      <  2:427704/150‑1 (MQ=255)
                                                                        cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACACACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCa                                                                      <  1:955467/150‑1 (MQ=255)
                                                                        cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCa                                                                      <  1:597297/150‑1 (MQ=255)
                                                                                 gCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAg                                                             >  4:338432/1‑150 (MQ=255)
                                                                                  cGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgc                                                            >  2:80920/1‑150 (MQ=255)
                                                                                     aaaCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTc                                                                                                                                             <  1:210284/67‑1 (MQ=255)
                                                                                     aaaCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTc                                                                                                                                             >  2:210284/1‑67 (MQ=255)
                                                                                        cTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGctcctc                                                                                                                              <  2:317693/79‑1 (MQ=255)
                                                                                        cTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGctcctc                                                                                                                              >  1:317693/1‑79 (MQ=255)
                                                                                        cTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaa                                                                          <  4:462687/130‑1 (MQ=255)
                                                                                        cTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTaaaa                                                                          >  3:462687/1‑130 (MQ=255)
                                                                                          cGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCtt                                                    <  1:80920/150‑1 (MQ=255)
                                                                                           gCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTc                                                   <  3:657387/150‑1 (MQ=255)
                                                                                           gCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTc                                                   >  1:109763/1‑150 (MQ=255)
                                                                                             cAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTg                                                 <  3:166154/150‑1 (MQ=255)
                                                                                                  tttAACACCTGGTTGCTGAGCGCAAATGCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCg                                            <  4:312188/150‑1 (MQ=255)
                                                                                                   ttAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGc                                           >  1:752074/1‑150 (MQ=255)
                                                                                                    tAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCt                                          >  1:342497/1‑150 (MQ=255)
                                                                                                    tAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCt                                          >  2:131446/1‑150 (MQ=255)
                                                                                                     aaCACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTa                                         <  3:174020/150‑1 (MQ=255)
                                                                                                          cTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATg                                    >  1:538487/1‑150 (MQ=255)
                                                                                                           tgcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATgg                                   >  1:894945/1‑150 (MQ=255)
                                                                                                            gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGga                                  >  1:165701/1‑150 (MQ=255)
                                                                                                              ttgctAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgatt                                <  2:469813/150‑1 (MQ=255)
                                                                                                               tgctAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc                               <  1:927248/150‑1 (MQ=255)
                                                                                                                gctAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGc                                                                                                                                   >  1:349783/1‑50 (MQ=255)
                                                                                                                gctAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGc                                                                                                                                   <  2:349783/50‑1 (MQ=255)
                                                                                                                  tAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCa                                                                      >  1:1005514/1‑108 (MQ=255)
                                                                                                                  tAAGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCa                                                                      <  2:1005514/108‑1 (MQ=255)
                                                                                                                   aaGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGttt                                                                  >  2:193435/1‑111 (MQ=255)
                                                                                                                   aaGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGttt                                                                  <  1:193435/111‑1 (MQ=255)
                                                                                                                   aaGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgatt                           <  3:379960/150‑1 (MQ=255)
                                                                                                                   aaGCGCAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgatt                           <  4:266458/150‑1 (MQ=255)
                                                                                                                     gcgcAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTcc                         >  4:246341/1‑150 (MQ=255)
                                                                                                                      cgcAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTTATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTccc                        <  3:663264/150‑1 (MQ=255)
                                                                                                                      cgcAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTccc                        >  1:86168/1‑150 (MQ=255)
                                                                                                                       gcAAATTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGACGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTcccc                       >  2:960174/1‑150 (MQ=255)
                                                                                                                         aaaTTCAATAAACTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTc                     >  2:532708/1‑150 (MQ=255)
                                                                                                                                  aaaCTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGc                                                <  1:237422/114‑1 (MQ=255)
                                                                                                                                  aaaCTGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGc                                                >  2:237422/1‑114 (MQ=255)
                                                                                                                                      tGGCGCTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaa        <  3:98685/150‑1 (MQ=255)
                                                                                                                                        gcgcTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcat                                                                                        <  1:429382/68‑1 (MQ=255)
                                                                                                                                        gcgcTGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACtcat                                                                                        >  2:429382/1‑68 (MQ=255)
                                                                                                                                            tGATATTGTTTCATGCCTTGGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGAt  <  2:204888/150‑1 (MQ=255)
                                                                                                                                            tGATATTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGAt  <  2:496599/150‑1 (MQ=255)
                                                                                                                                              atatTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaaa                                                                                                    >  2:472076/1‑50 (MQ=255)
                                                                                                                                              atatTGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACaaaaaa                                                                                                    <  1:472076/50‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AACGGCCTAACAGTGCCAGATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGAT  >  NC_000913_3_bme_pgi/3815615‑3815904

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