Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 3,316,893 C→T E212K (GAA→AAA)  nusA ← transcription termination/antitermination L factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi3,316,8930CT98.2% 379.8 / ‑3.7 108E212K (GAA→AAA) nusAtranscription termination/antitermination L factor
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (51/55);  total (52/56)
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.99e-01

Only 100 of 108 total aligned reads displayed.

CTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTCGCGC  >  NC_000913_3_bme_pgi/3316747‑3317038
                                                                                                                                                  |                                                                                                                                                 
cTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAAt                                                                                                                                                <  4:311368/150‑1 (MQ=255)
   aaCACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACg                                                                                                                                             <  3:391017/150‑1 (MQ=255)
   aaCACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACg                                                                                                                                             <  3:3778/150‑1 (MQ=255)
     cacGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGa                                                                                                                                           >  1:213403/1‑150 (MQ=255)
      acGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGaa                                                                                                                                          <  4:738709/150‑1 (MQ=255)
         cgcgcCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAg                                                                                                                                       >  1:1242316/1‑150 (MQ=255)
           cgcCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGtt                                                                                                                                     >  2:878037/1‑150 (MQ=255)
            gcCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTc                                                                                                                                    <  2:213403/150‑1 (MQ=255)
                 gCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCa                                                                                                                               <  3:739768/150‑1 (MQ=255)
                   aTACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGc                                                                                                                             >  4:405733/1‑150 (MQ=255)
                   aTACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGc                                                                                                                             >  2:1190830/1‑150 (MQ=255)
                   aTACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGc                                                                                                                             >  1:1153241/1‑150 (MQ=255)
                     aCCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGa                                                                                                                                  <  1:121091/143‑1 (MQ=255)
                     aCCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGa                                                                                                                                  >  2:121091/1‑143 (MQ=255)
                         aCGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTcc                                                                                                                       >  2:873176/1‑150 (MQ=255)
                          cGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGtt                                                                                                                                     <  4:468510/135‑1 (MQ=255)
                          cGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGtt                                                                                                                                     >  3:468510/1‑135 (MQ=255)
                             aaGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCt                                                                                                                   >  2:874302/1‑150 (MQ=255)
                                   cTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAAc                                                                                                             <  1:874302/150‑1 (MQ=255)
                                     aCCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGa                                                                                                           >  1:916607/1‑150 (MQ=255)
                                           tCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACg                                                                                                     <  1:106152/150‑1 (MQ=255)
                                            cGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGa                                                                                                    <  4:210255/150‑1 (MQ=255)
                                             gATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGaa                                                                                                   <  2:868481/150‑1 (MQ=255)
                                                  gTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGtt                                                                                              >  3:135804/1‑150 (MQ=255)
                                                  gTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGtt                                                                                              >  2:832328/1‑150 (MQ=255)
                                                    ttATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGaa                                                                                                   >  2:1131284/1‑143 (MQ=255)
                                                    ttATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGaa                                                                                                   <  1:1131284/143‑1 (MQ=255)
                                                      aTCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATCTCAATCACTTCTTCGCCGATTTCTGGCACCTTAATACGGAACAGTTCGATCAGCATCTCCGCCTTGGAACGACGGACGAACaacgacaa                                                                                          >  2:684293/1‑143 (MQ=255)
                                                         gTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGcc                                                                                       >  1:620470/1‑150 (MQ=255)
                                                          ttGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTgcgcgc                                                                                        >  1:1183821/1‑148 (MQ=255)
                                                          ttGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCa                                                                                      >  3:121347/1‑150 (MQ=255)
                                                          gtGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTgcgcgc                                                                                        <  2:1183821/147‑1 (MQ=255)
                                                            ggTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCAcg                                                                                    >  1:379697/1‑150 (MQ=255)
                                                             gTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCAcgc                                                                                   >  3:414490/1‑150 (MQ=255)
                                                             gTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCAcgc                                                                                   <  4:362148/150‑1 (MQ=255)
                                                                ttCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCt                                                                                >  1:1198178/1‑150 (MQ=255)
                                                                     cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCt                                                                                                                   >  1:1093068/1‑110 (MQ=255)
                                                                     cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCt                                                                                                                   <  2:1093068/110‑1 (MQ=255)
                                                                     cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTgc                                                                                            <  1:556384/133‑1 (MQ=255)
                                                                     cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTgc                                                                                            >  2:556384/1‑133 (MQ=255)
                                                                     cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                                           >  4:48913/1‑150 (MQ=255)
                                                                     cgcgACTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                                           >  4:625936/1‑150 (MQ=255)
                                                                      gcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCggg                                                                          <  1:97183/150‑1 (MQ=255)
                                                                       cgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAg                                                                                                <  3:37647/127‑1 (MQ=255)
                                                                       cgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAg                                                                                                >  4:37647/1‑127 (MQ=255)
                                                                          ttttCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGaa                                                                      <  3:103418/150‑1 (MQ=255)
                                                                           tttCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAAc                                                                     <  3:48913/150‑1 (MQ=255)
                                                                             tCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGAt                                                                                                                                 >  4:507817/1‑88 (MQ=255)
                                                                             tCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGAt                                                                                                                                 <  3:507817/88‑1 (MQ=255)
                                                                                 aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGTCTTGGAACGAGTGACGaa                                                                                                   <  2:776466/114‑1 (MQ=255)
                                                                                 aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGaa                                                                                                   >  1:776466/1‑114 (MQ=255)
                                                                                 aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATa                                                               >  1:552675/1‑150 (MQ=255)
                                                                                 aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATa                                                               <  2:272715/150‑1 (MQ=255)
                                                                                 aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATa                                                               <  2:798751/150‑1 (MQ=255)
                                                                                  cGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATag                                                              <  4:514082/150‑1 (MQ=255)
                                                                                   gagaACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATaga                                                             <  4:135804/150‑1 (MQ=255)
                                                                                    agaACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATagag                                                            >  2:261169/1‑150 (MQ=255)
                                                                                      aaCCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGATTAGAGCa                                                          <  2:1242316/150‑1 (MQ=255)
                                                                                        cccGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGa                                                                                                               >  3:489461/1‑95 (MQ=255)
                                                                                        cccGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGa                                                                                                               <  4:489461/95‑1 (MQ=255)
                                                                                        cccGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACg                                                                    >  4:286381/1‑138 (MQ=255)
                                                                                        cccGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACg                                                                    <  3:286381/138‑1 (MQ=255)
                                                                                         ccGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGc                                                       <  3:645666/150‑1 (MQ=255)
                                                                                         ccGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGc                                                       <  1:873176/150‑1 (MQ=255)
                                                                                                 cgAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAAcg                                               >  3:75895/1‑150 (MQ=255)
                                                                                                   aGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAAcgcg                                             <  2:552675/150‑1 (MQ=255)
                                                                                                       gctgctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGc                                                                                                                             >  4:385769/1‑66 (MQ=255)
                                                                                                       gctgctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGaa                                                                 >  3:32439/1‑126 (MQ=255)
                                                                                                       gctgctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGaa                                                                 <  4:32439/126‑1 (MQ=255)
                                                                                                           ctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGTCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAg                                     >  2:691829/1‑150 (MQ=255)
                                                                                                           ctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAg                                     >  1:689656/1‑150 (MQ=255)
                                                                                                            tttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGCGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAgg                                    >  2:717627/1‑150 (MQ=255)
                                                                                                            tttAATTTCAATCACTTCTTCGCCGATTTCTGACACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGGGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAgg                                    <  1:522737/150‑1 (MQ=255)
                                                                                                               aaTTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGtt                                                                                              <  1:421920/89‑1 (MQ=255)
                                                                                                               aaTTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGtt                                                                                              >  2:421920/1‑89 (MQ=255)
                                                                                                                aTTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCgg                                <  3:405733/150‑1 (MQ=255)
                                                                                                                 tttCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                                           <  1:863905/106‑1 (MQ=255)
                                                                                                                 tttCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                                           >  2:863905/1‑106 (MQ=255)
                                                                                                                  ttCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGaa                              <  2:1153241/150‑1 (MQ=255)
                                                                                                                  ttCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGaa                              <  2:799169/150‑1 (MQ=255)
                                                                                                                  ttCAATCACTGCTGCGCCGAGTGCTGGCACTTTAATACAGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCTCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGaa                              <  2:630791/150‑1 (MQ=255)
                                                                                                                     aaTCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGtt                           >  4:554416/1‑150 (MQ=255)
                                                                                                                      aTCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGttt                          >  1:605024/1‑150 (MQ=255)
                                                                                                                      aTCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGttt                          <  1:832328/150‑1 (MQ=255)
                                                                                                                       tCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAg                                                                                                          >  4:475292/1‑69 (MQ=255)
                                                                                                                       tCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAg                                                                                                          <  3:475292/69‑1 (MQ=255)
                                                                                                                        cACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTc                        <  2:689656/150‑1 (MQ=255)
                                                                                                                              ttcGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCgg                  <  1:6441/150‑1 (MQ=255)
                                                                                                                               tcGCCGATTTCTTGCACTTTAATACGAAACGCGCCGATCATCACCGCCCACTGGTAACTAGTGACCAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgc                 <  3:644951/150‑1 (MQ=255)
                                                                                                                               tcGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgc                 <  4:192953/150‑1 (MQ=255)
                                                                                                                                 gCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgcag               <  4:121347/150‑1 (MQ=255)
                                                                                                                                 gCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgcag               >  3:399613/1‑150 (MQ=255)
                                                                                                                                 gCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgcag               >  2:224150/1‑150 (MQ=255)
                                                                                                                                  ccGATTTCTGGCACTTTAAAACGGAACAGGTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCCCGAACGCGGTCGCCAGGGCGGCAGTTTTCACGCGgcagc              >  4:456684/1‑150 (MQ=255)
                                                                                                                                    gATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCat            <  1:1190830/150‑1 (MQ=255)
                                                                                                                                       ttCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATc         >  1:370782/1‑150 (MQ=255)
                                                                                                                                          tGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATa                                                               <  3:533589/93‑1 (MQ=255)
                                                                                                                                          tGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATa                                                               >  4:533589/1‑93 (MQ=255)
                                                                                                                                          tGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGTCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTc      <  2:916607/150‑1 (MQ=255)
                                                                                                                                              aCTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTCgcgc  >  4:461224/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTCGCGC  >  NC_000913_3_bme_pgi/3316747‑3317038

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