Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,186,605 A→C H419P (CAC→CCC)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,186,6050AC97.5% 407.6 / 0.4 118H419P (CAC→CCC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base C (54/61);  total (57/61)
Fisher's exact test for biased strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 121 total aligned reads displayed.

AGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGG  >  NC_000913_3_bme_pgi/4186458‑4186748
                                                                                                                                                   |                                                                                                                                               
agaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAAcccc                                                                                                                                               <  4:505530/150‑1 (MQ=255)
   aaaTGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCgg                                                                                                                                            <  2:889817/150‑1 (MQ=255)
    aaTGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGt                                                                                                                                           >  1:185118/1‑150 (MQ=255)
         cACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTActgc                                                                                                                                      <  2:185118/150‑1 (MQ=255)
             ggggCTGTGCAGACCGTTCATCTGCGGCACGCTGTAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAATATGGTTTATCGCGAAGACGCTGTCGTTTGGGTTATCCTGGAGTAATTTATCCGCGAACCCCCGGTACTGCTGac                                                                                                                                  <  4:136292/147‑2 (MQ=255)
             ggAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACcccc                                                                                                                                              >  3:483312/1‑138 (MQ=255)
             ggAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACcccc                                                                                                                                              <  4:483312/138‑1 (MQ=255)
              gAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAAc                                                                                                                                 >  1:1016043/1‑150 (MQ=255)
                 cTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCTTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGt                                                                                                                              >  2:1033449/1‑150 (MQ=255)
                 cTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGt                                                                                                                              >  4:602433/1‑150 (MQ=255)
                    ttCCAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCa                                                                                                                           >  2:69342/1‑150 (MQ=255)
                    ttCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCa                                                                                                                           <  1:879560/150‑1 (MQ=255)
                     tCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCAc                                                                                                                          <  1:1033449/150‑1 (MQ=255)
                      cAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCGCCATCAAGGATGAGACGAAAAACGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCAcc                                                                                                                         <  1:597914/150‑1 (MQ=255)
                      cAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCCTGCTACCACAACTACGGCTGTGAAGAAAATGGTTGAACGCGATGAATCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCAcc                                                                                                                         <  1:597821/150‑1 (MQ=255)
                       aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGccccc                                                                                                                        >  4:429428/1‑149 (MQ=255)
                            gTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGCAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGActct                                                                                                                   >  1:268788/1‑150 (MQ=255)
                            gTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGActct                                                                                                                   <  3:490902/150‑1 (MQ=255)
                             ttCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTg                                                                                                                  >  4:704991/1‑150 (MQ=255)
                             gtCATCTACGGCAAGCTGGAACGGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTg                                                                                                                  <  2:44510/149‑1 (MQ=255)
                                 tCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCAcc                                                                                                              <  2:331692/150‑1 (MQ=255)
                                 tCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCAcc                                                                                                              >  3:511541/1‑150 (MQ=255)
                                 tCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCAcc                                                                                                              >  2:93121/1‑150 (MQ=255)
                                  cTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCg                                                                                                             >  4:603402/1‑150 (MQ=255)
                                        cAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTggg                                                                                                       <  1:513909/150‑1 (MQ=255)
                                        cAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTggg                                                                                                       >  2:1309433/1‑150 (MQ=255)
                                             tGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATcc                                                                                                  >  3:139190/1‑150 (MQ=255)
                                              ggAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGtcggtg                                                                                                 <  4:197949/150‑6 (MQ=255)
                                                   tGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAt                                                                                            >  1:222416/1‑150 (MQ=255)
                                                     cGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAttt                                                                                          <  3:234625/150‑1 (MQ=255)
                                                      gTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTg                                                                                         <  1:995405/150‑1 (MQ=255)
                                                         gTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAAc                                                                                      <  3:704991/150‑1 (MQ=255)
                                                           cTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCg                                                                                    >  2:446158/1‑150 (MQ=255)
                                                             tGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGt                                                                                  <  1:446158/150‑1 (MQ=255)
                                                              gCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTa                                                                                 <  2:780144/150‑1 (MQ=255)
                                                                 accaccATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGAc                                                                                                                      >  2:524511/1‑110 (MQ=255)
                                                                 accaccATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGAc                                                                                                                      <  1:524511/110‑1 (MQ=255)
                                                                 accaccATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTg                                                                              >  1:441935/1‑150 (MQ=255)
                                                                  ccaccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGa                                                                             >  2:1124610/1‑150 (MQ=255)
                                                                   caccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGAt                                                                            >  2:1034675/1‑150 (MQ=255)
                                                                     ccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTActgc                                                                                                                                      <  4:244844/90‑1 (MQ=255)
                                                                     ccaTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTActgc                                                                                                                                      >  3:244844/1‑90 (MQ=255)
                                                                       aTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATcc                                                                                                  <  1:912262/124‑1 (MQ=255)
                                                                       aTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGaa                                                                        <  4:645862/150‑1 (MQ=255)
                                                                       aTTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGAATTTGAACCGGTACTGATCGaa                                                                        >  3:695439/1‑150 (MQ=255)
                                                                        ttAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAg                                                                       <  4:139190/150‑1 (MQ=255)
                                                                        ttAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAg                                                                       >  4:307823/1‑150 (MQ=255)
                                                                          aaaGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGt                                                                     <  3:602433/150‑1 (MQ=255)
                                                                              cTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAg                                                                 <  3:307823/150‑1 (MQ=255)
                                                                              cTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAg                                                                 <  3:15765/150‑1 (MQ=255)
                                                                                     gaaAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTa                                                                                                                                          <  1:390830/70‑1 (MQ=255)
                                                                                      aaaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAg                                                                 >  4:750144/1‑142 (MQ=255)
                                                                                      aaaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAg                                                                 <  3:750144/142‑1 (MQ=255)
                                                                                       aaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGc                                                        >  2:302702/1‑150 (MQ=255)
                                                                                       aaaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGc                                                        >  2:302580/1‑150 (MQ=255)
                                                                                        aaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCt                                                       >  4:264925/1‑150 (MQ=255)
                                                                                        aaTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCt                                                       <  2:759538/150‑1 (MQ=255)
                                                                                         aTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGAct                                                                                                                     <  3:146397/87‑1 (MQ=255)
                                                                                         aTGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGAct                                                                                                                     >  4:146397/1‑87 (MQ=255)
                                                                                           ggTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCAc                                                                                                                          >  4:774627/1‑80 (MQ=255)
                                                                                             ttGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCAcc                                                  <  1:69342/150‑1 (MQ=255)
                                                                                               gAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGt                                                                                  <  4:70685/116‑1 (MQ=255)
                                                                                               gAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGt                                                                                  >  3:70685/1‑116 (MQ=255)
                                                                                                aGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGc                                               <  2:441935/150‑1 (MQ=255)
                                                                                                 gcgcGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCt                                              <  1:93121/150‑1 (MQ=255)
                                                                                                  cgcgAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAtt                                                                                           <  2:531432/104‑1 (MQ=255)
                                                                                                         aagCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTActgct                                                                                                                                     >  1:192861/1‑55 (MQ=255)
                                                                                                         aagCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTActgct                                                                                                                                     <  2:192861/55‑1 (MQ=255)
                                                                                                         aagCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTtgtg                                      <  1:1034675/150‑1 (MQ=255)
                                                                                                          agCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGc                                     <  4:637339/150‑1 (MQ=255)
                                                                                                             tGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                                                                                        >  3:286043/1‑80 (MQ=255)
                                                                                                             tGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                                                                                        <  4:286043/80‑1 (MQ=255)
                                                                                                              gTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGAATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCa                                 <  4:695439/150‑1 (MQ=255)
                                                                                                               tCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTg                                                                                                         >  1:308613/1‑77 (MQ=255)
                                                                                                               tCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTg                                                                                                         <  2:308613/77‑1 (MQ=255)
                                                                                                                 gTTTGGGATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCatat                              >  1:964056/1‑150 (MQ=255)
                                                                                                                     gggATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACg                          <  4:511541/150‑1 (MQ=255)
                                                                                                                      ggATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGa                                                                                                                       <  2:1266938/56‑1 (MQ=255)
                                                                                                                      ggATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGa                                                                                                                       >  1:1266938/1‑56 (MQ=255)
                                                                                                                       gATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCg                                                                                    <  2:797929/90‑1 (MQ=255)
                                                                                                                       gATATCCTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCg                                                                                    >  1:797929/1‑90 (MQ=255)
                                                                                                                            ccTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACtt                   >  3:522575/1‑150 (MQ=255)
                                                                                                                             cTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTc                  <  3:264925/150‑1 (MQ=255)
                                                                                                                             cTGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTc                  >  1:1313095/1‑150 (MQ=255)
                                                                                                                              tGGACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCg                 >  4:529492/1‑150 (MQ=255)
                                                                                                                                gACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGc                                  >  2:1192753/1‑131 (MQ=255)
                                                                                                                                gACGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGAt               >  1:103232/1‑150 (MQ=255)
                                                                                                                                 aCGAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATg              >  2:1243560/1‑150 (MQ=255)
                                                                                                                                   gAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGt            >  3:398132/1‑150 (MQ=255)
                                                                                                                                   gAAGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGt            <  2:1016043/150‑1 (MQ=255)
                                                                                                                                     aGTTATCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGa          >  1:31095/1‑150 (MQ=255)
                                                                                                                                         aTCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAAc                                                                                      >  1:637882/1‑70 (MQ=255)
                                                                                                                                         aTCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAg      >  2:25800/1‑150 (MQ=255)
                                                                                                                                          tCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGt            <  1:28230/143‑1 (MQ=255)
                                                                                                                                          tCCGCGAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGt            >  2:28230/1‑143 (MQ=255)
                                                                                                                                            cgcgAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATcc                                                                                                  >  1:125677/1‑55 (MQ=255)
                                                                                                                                            cgcgAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATcc                                                                                                  <  2:125677/55‑1 (MQ=255)
                                                                                                                                             gcgAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                                                                                        <  3:794934/48‑1 (MQ=255)
                                                                                                                                             gcgAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATgg  >  2:1203885/1‑150 (MQ=255)
                                                                                                                                               gAACCCCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGAt               >  1:1176283/1‑135 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGAAGGTAAAGCTATCCAGCTGCACCCGCTGGTTTGTGCGGCATATAACGCCGACTTCGATGGTGACCAGATGG  >  NC_000913_3_bme_pgi/4186458‑4186748

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: