Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,187,550 C→T A734V (GCG→GTG)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,187,5500CT99.0% 365.4 / NA 96A734V (GCG→GTG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (48/46);  minor base G (0/1);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGT  >  NC_000913/4187402‑4187695
                                                                                                                                                    |                                                                                                                                                 
atCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTgg                                                                                                                                                 <  2:406393/151‑1 (MQ=255)
atCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGCGCGCGTGGTTCTGTgg                                                                                                                                                 >  2:930026/1‑151 (MQ=255)
   tGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGcac                                                                                                                                              <  1:352929/151‑1 (MQ=255)
   tGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGcac                                                                                                                                              <  1:838619/151‑1 (MQ=255)
   tGGGCTGCGGCGAACGATCGGGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGcac                                                                                                                                              <  2:929908/151‑1 (MQ=255)
      gCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAg                                                                                                                                            <  2:596283/150‑1 (MQ=255)
        tGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGAtt                                                                                                                                         <  1:259785/151‑1 (MQ=255)
            gcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATTGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGt                                                                                                                                      >  1:808743/1‑150 (MQ=255)
             cgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTc                                                                                                                                     <  2:160773/150‑1 (MQ=255)
             cgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTc                                                                                                                                     <  2:954298/150‑1 (MQ=255)
             cgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTc                                                                                                                                     <  2:526953/150‑1 (MQ=255)
                 cGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCtt                                                                                                                                >  2:924360/1‑151 (MQ=255)
                   aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTg                                                                                                                               >  2:300610/1‑150 (MQ=255)
                   aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTg                                                                                                                               >  2:455415/1‑150 (MQ=255)
                    tCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCT‑‑GGCACAGATTCGTCAGCTTGCggg                                                                                                                           <  1:989301/151‑1 (MQ=255)
                     cGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTg                                                                                                                               >  2:247214/1‑148 (MQ=255)
                     cGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTg                                                                                                                               >  1:378160/1‑148 (MQ=255)
                        gtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTa                                                                                                                         >  2:406398/1‑151 (MQ=255)
                          aTCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATg                                                                                                                       >  2:792797/1‑151 (MQ=255)
                          aTCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATg                                                                                                                       <  1:237881/151‑1 (MQ=255)
                            ccAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCg                                                                                                                     <  2:37038/151‑1 (MQ=255)
                            ccAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCg                                                                                                                     <  1:978807/151‑1 (MQ=255)
                             cAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGt                                                                                                                    <  2:574871/151‑1 (MQ=255)
                             cAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGt                                                                                                                    <  2:126461/151‑1 (MQ=255)
                             cAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGt                                                                                                                    <  2:368435/151‑1 (MQ=255)
                             cAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGt                                                                                                                    <  2:51350/151‑1 (MQ=255)
                                 gCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGt                                                                                                                 >  1:455356/1‑150 (MQ=255)
                                 gCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGt                                                                                                                 >  2:553720/1‑150 (MQ=255)
                                  cGATGTTGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGt                                                                                                                 <  2:965556/149‑1 (MQ=255)
                                      gatgGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTAAACAGCATCTACATGATGGCCGACTCCGGGGCGCGTGGGTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATg                                                                                                           >  2:282423/1‑151 (MQ=255)
                                       atgGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATgg                                                                                                          >  1:143293/1‑151 (MQ=255)
                                        tgGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATgg                                                                                                          >  2:271014/1‑150 (MQ=255)
                                               ccTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCg                                                                                                  >  1:621369/1‑151 (MQ=255)
                                                   cAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATg                                                                                              <  1:785893/151‑1 (MQ=255)
                                                     aaCTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATgg                                                                                             >  1:669244/1‑150 (MQ=255)
                                                     aaCTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGAAGGGTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATgg                                                                                             >  2:308/1‑150 (MQ=255)
                                                       cTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCt                                                                                           >  2:179789/1‑150 (MQ=255)
                                                           aaCCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCcatc                                                                                      <  1:491699/151‑1 (MQ=255)
                                                             ccGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCcatcat                                                                                    >  1:520929/1‑151 (MQ=255)
                                                             ccGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCcatcat                                                                                    <  1:163239/151‑1 (MQ=255)
                                                               gTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCg                                                                                  >  2:805593/1‑151 (MQ=255)
                                                               gTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCg                                                                                  >  1:11397/1‑151 (MQ=255)
                                                                    tAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACg                                                                             >  1:353029/1‑151 (MQ=255)
                                                                        cGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCaa                                                                         <  1:62481/151‑1 (MQ=255)
                                                                             cGGCCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCAc                                                                     >  1:271581/1‑150 (MQ=255)
                                                                               gTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACcg                                                                   >  2:150890/1‑150 (MQ=255)
                                                                                 cAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGa                                                                >  2:223351/1‑151 (MQ=255)
                                                                                   ggAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGaa                                                               >  1:57736/1‑150 (MQ=255)
                                                                                       gagaAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTc                                                           >  2:438020/1‑150 (MQ=255)
                                                                                       gagaAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTc                                                           >  2:804220/1‑150 (MQ=255)
                                                                                        agaAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTc                                                           <  1:547110/149‑1 (MQ=255)
                                                                                         gaAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCgg                                                        >  1:48897/1‑150 (MQ=255)
                                                                                            gCAGGTTTCCTTAAACAGCCTCTACATGATGGCCGCCTACGGTGCGAGTGGTTCTGTGGCCCAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGCTGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGa                                                      <  2:21067/150‑1 (MQ=255)
                                                                                             cAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGa                                                      <  1:255687/149‑1 (MQ=255)
                                                                                               ggTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGt                                                  >  1:169763/1‑151 (MQ=255)
                                                                                                gTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAg                                                    <  1:485379/148‑1 (MQ=255)
                                                                                                   tCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGa                                              >  2:429362/1‑151 (MQ=255)
                                                                                                   tCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGa                                              >  2:552244/1‑151 (MQ=255)
                                                                                                   tCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGa                                              >  1:190638/1‑151 (MQ=255)
                                                                                                   tCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGa                                              >  2:864367/1‑151 (MQ=255)
                                                                                                     cTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAAc                                            <  2:57744/151‑1 (MQ=255)
                                                                                                     cTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAAc                                            <  2:377171/151‑1 (MQ=255)
                                                                                                     cTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAAc                                            <  1:285701/151‑1 (MQ=255)
                                                                                                      ttCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACg                                           <  1:78913/151‑1 (MQ=255)
                                                                                                             gCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCa                                    >  1:554615/1‑151 (MQ=255)
                                                                                                              cATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTc                                      <  1:442743/148‑1 (MQ=255)
                                                                                                              cATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTc                                      <  2:125337/148‑1 (MQ=255)
                                                                                                              cATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAg                                   <  2:464969/151‑1 (MQ=255)
                                                                                                              cATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAg                                   >  1:78597/1‑151 (MQ=255)
                                                                                                                 cTATATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTa                                 >  1:806306/1‑150 (MQ=255)
                                                                                                                 cTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTa                                 >  2:323194/1‑150 (MQ=255)
                                                                                                                 cTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTa                                 >  2:151710/1‑150 (MQ=255)
                                                                                                                 cTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTa                                 >  1:50898/1‑150 (MQ=255)
                                                                                                                 cTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTa                                 >  1:1030037/1‑150 (MQ=255)
                                                                                                                 cTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTAc                                <  1:689325/151‑1 (MQ=255)
                                                                                                                  tACATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACt                               <  2:611051/151‑1 (MQ=255)
                                                                                                                   aCATGATGGCCGGCTCCGGTGCGCGTGGTTCTGGGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACtt                              <  2:320716/151‑1 (MQ=255)
                                                                                                                   aCATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACtt                              >  2:992685/1‑151 (MQ=255)
                                                                                                                   aCATGATGGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACtt                              <  2:577796/151‑1 (MQ=255)
                                                                                                                      tgatgGCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAt                           >  1:1054538/1‑151 (MQ=255)
                                                                                                                          ggCCGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAtctc                        >  2:830585/1‑150 (MQ=255)
                                                                                                                             cGACTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccacc                    <  2:455535/151‑1 (MQ=255)
                                                                                                                               aCTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccaccca                  >  1:782777/1‑151 (MQ=255)
                                                                                                                                cTCCGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAACCGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccacccac                 <  1:496914/151‑1 (MQ=255)
                                                                                                                                  ccGGTGCGCGTGGTTCTGTGGCGCAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACgg               <  1:1023321/151‑1 (MQ=255)
                                                                                                                                  ccGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACgg               <  2:795346/151‑1 (MQ=255)
                                                                                                                                  ccGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACgg               <  1:49198/151‑1 (MQ=255)
                                                                                                                                   cGGTGCGCGTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGt              >  1:770519/1‑151 (MQ=255)
                                                                                                                                         gcgTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGGGCTCGt        >  2:800757/1‑151 (MQ=255)
                                                                                                                                          cgTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTa       <  2:777963/151‑1 (MQ=255)
                                                                                                                                          cgTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTa       <  1:18707/151‑1 (MQ=255)
                                                                                                                                          cgTGGTTCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTa       <  1:1051101/151‑1 (MQ=255)
                                                                                                                                               ttCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCCCCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGt  <  2:830948/151‑1 (MQ=255)
                                                                                                                                               ttCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAgg   <  1:485862/150‑1 (MQ=255)
                                                                                                                                               ttCTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAgg   <  1:429204/150‑1 (MQ=255)
                                                                                                                                                 cTGTGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGt  >  2:441622/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGT  >  NC_000913/4187402‑4187695

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