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,8930CT99.0% 359.3 / NA 102E212K (GAA→AAA) nusAtranscription termination/antitermination L factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (52/49);  minor base A (1/0);  total (53/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 102 total aligned reads displayed.

TGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTCGCGC  >  NC_000913_3_bme_pgi/3316748‑3317038
                                                                                                                                                 |                                                                                                                                                 
tGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATAc                                                                                                                                              <  2:593061/151‑1 (MQ=255)
     acGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAAc                                                                                                                                         <  2:557596/151‑1 (MQ=255)
       gcgcgcCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACa                                                                                                                                        >  3:256476/1‑150 (MQ=255)
           gcCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCg                                                                                                                                   >  1:330040/1‑151 (MQ=255)
               cGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCa                                                                                                                               >  1:256734/1‑151 (MQ=255)
                 cATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAg                                                                                                                              <  4:167890/150‑1 (MQ=255)
                 cATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGc                                                                                                                             >  1:66231/1‑151 (MQ=255)
                 cATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGc                                                                                                                             >  2:761184/1‑151 (MQ=255)
                     ccTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCAttt                                                                                                                         >  1:370806/1‑151 (MQ=255)
                     ccTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCAttt                                                                                                                         >  1:711568/1‑151 (MQ=255)
                      cTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCACCGATATACAGCACTTAAATACGGAACAGTTCGATCAGCATTTc                                                                                                                        >  1:277484/1‑151 (MQ=255)
                       tACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTcc                                                                                                                       >  1:739300/1‑151 (MQ=255)
                          gCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGc                                                                                                                    >  2:321000/1‑151 (MQ=255)
                          gCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGc                                                                                                                    <  1:880501/151‑1 (MQ=255)
                             aGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTg                                                                                                                 <  2:739300/151‑1 (MQ=255)
                             aGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTg                                                                                                                 <  2:366331/151‑1 (MQ=255)
                             aGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTg                                                                                                                 >  1:705750/1‑151 (MQ=255)
                              gCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTgg                                                                                                                >  1:587908/1‑151 (MQ=255)
                              gCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTg                                                                                                                 <  4:11048/150‑1 (MQ=255)
                                  cTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACg                                                                                                            <  1:265439/151‑1 (MQ=255)
                                  cTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACg                                                                                                            >  2:793936/1‑151 (MQ=255)
                                    aCCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGa                                                                                                           <  4:256476/150‑1 (MQ=255)
                                    aCCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAg                                                                                                          >  1:728834/1‑151 (MQ=255)
                                     ccGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAg                                                                                                          >  3:579935/1‑150 (MQ=255)
                                     ccGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGt                                                                                                         <  1:21316/151‑1 (MQ=255)
                                      cGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTg                                                                                                        >  2:560742/1‑151 (MQ=255)
                                       ggATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGa                                                                                                       <  1:436367/151‑1 (MQ=255)
                                       ggATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGa                                                                                                       <  2:910970/151‑1 (MQ=255)
                                        gATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGaa                                                                                                                                          >  3:98546/1‑115 (MQ=255)
                                        gATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGa                                                                                                       >  4:337619/1‑150 (MQ=255)
                                         aTCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACg                                                                                                     <  1:761184/151‑1 (MQ=255)
                                            gATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAAc                                                                                                  >  1:358508/1‑151 (MQ=255)
                                             aTACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGAc                                                                                                      >  3:419528/1‑146 (MQ=255)
                                             aTACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGAc                                                                                                      <  4:419528/146‑1 (MQ=255)
                                             aTACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACa                                                                                                 >  1:796507/1‑151 (MQ=255)
                                                 gTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTg                                                                                             <  2:14929/151‑1 (MQ=255)
                                                  tttATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGttcc                                                                                            <  2:587908/151‑3 (MQ=255)
                                                      tCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTgcgcg                                                                                         <  4:28273/150‑1 (MQ=255)
                                                        gTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATa                                                                                                                                               <  4:167251/94‑1 (MQ=255)
                                                        gTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATa                                                                                                                                               >  3:167251/1‑94 (MQ=255)
                                                        gTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCa                                                                                      >  1:288007/1‑151 (MQ=255)
                                                            gTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCCTTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCAcgcg                                                                                  <  2:778631/151‑1 (MQ=255)
                                                             ttttCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACgcgc                                                                                 <  1:684142/151‑1 (MQ=255)
                                                                   ccGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCa                                                                                                                            <  2:411574/102‑1 (MQ=255)
                                                                   ccGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCa                                                                                                                            >  1:411574/1‑102 (MQ=255)
                                                                   ccGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                                           >  2:796387/1‑151 (MQ=255)
                                                                    cgcgATTTTCGCACGAGAACCCGTATCGCTAGCCTCTGCTTTAATTTTAATCACTTCTTCGCCGTTTGCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACTAACAGTTGCGCGCCACGCGCTTCCggg                                                                          >  2:810488/1‑151 (MQ=255)
                                                                    cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                                           >  3:85096/1‑150 (MQ=255)
                                                                    cgcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCgg                                                                           >  4:398855/1‑150 (MQ=255)
                                                                      cgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGc                                                                         >  4:557889/1‑150 (MQ=255)
                                                                          tttCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAAc                                                                     <  3:157577/150‑1 (MQ=255)
                                                                               cACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACgg                                                                                                                                            <  2:487422/74‑1 (MQ=255)
                                                                               cACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACgg                                                                                                                                            >  1:487422/1‑74 (MQ=255)
                                                                                aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATag                                                              <  2:330040/151‑1 (MQ=255)
                                                                                aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATag                                                              <  1:744123/151‑1 (MQ=255)
                                                                                 cGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATaga                                                             >  2:333534/1‑151 (MQ=255)
                                                                                       cccGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGc                                                       >  2:208915/1‑151 (MQ=255)
                                                                                        ccGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAAAAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGc                                                       >  3:155101/1‑150 (MQ=255)
                                                                                          ggATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGa                                                                                                           >  3:263197/1‑96 (MQ=255)
                                                                                          ggATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGa                                                                                                           <  4:263197/96‑1 (MQ=255)
                                                                                           gATCGCGAGCCGCTGCTTTAATTTTAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCAc                                                    >  4:122271/1‑150 (MQ=255)
                                                                                           gATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACg                                                   <  1:766220/151‑1 (MQ=255)
                                                                                            aTCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGa                                                  >  2:944915/1‑151 (MQ=255)
                                                                                             tCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGaa                                                 <  2:711568/151‑1 (MQ=255)
                                                                                                   gCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGt                                           >  2:146441/1‑151 (MQ=255)
                                                                                                       ctgctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGcc                                       <  1:560742/151‑1 (MQ=255)
                                                                                                         gctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAg                                     <  2:728834/151‑1 (MQ=255)
                                                                                                          ctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAgg                                    >  1:768512/1‑151 (MQ=255)
                                                                                                          ctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAg                                     >  4:421219/1‑150 (MQ=255)
                                                                                                           tttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTc                                                                                                                                    >  1:768039/1‑54 (MQ=255)
                                                                                                           tttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTc                                                                                                                                    <  2:768039/54‑1 (MQ=255)
                                                                                                            ttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCgg                                                                                                                     <  3:404936/68‑1 (MQ=255)
                                                                                                            ttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCgg                                                                                                                     >  4:404936/1‑68 (MQ=255)
                                                                                                            ttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGc                                  <  1:208915/151‑1 (MQ=255)
                                                                                                            ttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGc                                  <  2:796507/151‑1 (MQ=255)
                                                                                                            ttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGc                                  <  1:482951/151‑1 (MQ=255)
                                                                                                            ttAATTCCAATCACTTCTTCGCCGATATATCGCCCTTTAATACGGAACAGATCGATCAGCATTTCCGGCTTGCAACGACTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGACCGCGGTCGCCACGGc                                  >  2:432778/1‑151 (MQ=255)
                                                                                                             tAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGGGACGAACAGTTGCGCGCCACGCGCGTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGc                                  >  4:298881/1‑150 (MQ=255)
                                                                                                                 ttCAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGaa                              <  4:579802/150‑1 (MQ=255)
                                                                                                                 ttCAATCACTTCCTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGAGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAg                             <  1:26286/151‑1 (MQ=255)
                                                                                                                   cAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGtt                           >  1:385311/1‑151 (MQ=255)
                                                                                                                   cAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGtt                           >  2:286385/1‑151 (MQ=255)
                                                                                                                   cAATCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGt                            <  3:421219/150‑1 (MQ=255)
                                                                                                                      tCACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTc                        <  1:146441/151‑1 (MQ=255)
                                                                                                                       cACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCa                       <  2:256734/151‑1 (MQ=255)
                                                                                                                       cACTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCa                       >  2:166809/1‑151 (MQ=255)
                                                                                                                        aCTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCa                       <  4:579935/150‑1 (MQ=255)
                                                                                                                        aCTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCAc                      <  1:321000/151‑1 (MQ=255)
                                                                                                                        aCTTCTTCGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCAc                      <  2:705750/151‑1 (MQ=255)
                                                                                                                           tcttcGCCGATTTCTGGCACTTTAATACGGAACAGTTc                                                                                                                                    <  1:644102/38‑1 (MQ=255)
                                                                                                                           tcttcGCCGATTTCTGGCACTTTAATACGGAACAGTTc                                                                                                                                    >  2:644102/1‑38 (MQ=255)
                                                                                                                            cttcGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTg                                                                                                                 <  1:706732/56‑1 (MQ=255)
                                                                                                                            cttcGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTg                                                                                                                 >  2:706732/1‑56 (MQ=255)
                                                                                                                             ttcGCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCgg                  <  3:398855/150‑1 (MQ=255)
                                                                                                                                gCCGATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgcag               >  3:65868/1‑150 (MQ=255)
                                                                                                                                 ccGATTTCTGGCACTTTAATACGGAACAGTTCGGTCAGCATTTCCGGCTTGGAACGAGAGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGGAGTGTTCACGCGgcagc              >  4:269095/1‑150 (MQ=255)
                                                                                                                                   gATTTCTGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCata           <  1:796387/151‑1 (MQ=255)
                                                                                                                                         tGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTcg     >  2:228184/1‑151 (MQ=255)
                                                                                                                                         tGGCACTTTAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTcg     >  2:405282/1‑151 (MQ=255)
                                                                                                                                            cACTTTCATTCGGGACAGTTCGTTCAGCATTTCCGGCTTGGAACGAGTCACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGGGCACGCCACGAACGCGGTCGCCAGGGCGGCAGTTTTCACGCGGCAGCATAGCTTCgcgc  <  2:370806/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTCGCGC  >  NC_000913_3_bme_pgi/3316748‑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: