Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,536,027 C→A 7.0% A332A (GCC→GCA pncB → nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,536,0270CA 7.0% 268.7 / 6.8 86A332A (GCC→GCApncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (38/42);  new base A (3/3);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-01

GCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAGCCGCCAAGACCCA  >  NC_002947/5535882‑5536176
                                                                                                                                                 |                                                                                                                                                     
gCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTgc                                                                                                                                                   <  2:1013447/150‑1 (MQ=255)
     cgcCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGa                                                                                                                                              <  2:775732/150‑1 (MQ=255)
          tGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCAt                                                                                                                                         <  2:394770/150‑1 (MQ=255)
            aTTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCa                                                                                                                                       >  1:750625/1‑150 (MQ=255)
             ttCGGGCGAGCCGGTGGCCTGGGCGGAGAAGGCCATTGCCCATTGCCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTATCCGACGACCTCAACCTAACCCGCTCGCTGGAAAGTTTCTCCGCCCTTCGCGGACGCACCaa                                                                                                                                      >  2:1709642/1‑150 (MQ=255)
             ttCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCaa                                                                                                                                      >  1:587329/1‑150 (MQ=255)
             ttCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCaa                                                                                                                                      >  2:800941/1‑150 (MQ=255)
              tCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAAc                                                                                                                                     <  2:1251942/150‑1 (MQ=255)
                 ggCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTc                                                                                                                                  <  1:1046834/150‑1 (MQ=255)
                  gCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCa                                                                                                                                 <  2:1758238/150‑1 (MQ=255)
                   cGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAg                                                                                                                                <  1:325202/150‑1 (MQ=255)
                       ccGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCttt                                                                                                                            <  2:518143/150‑1 (MQ=255)
                        cGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTg                                                                                                                           >  2:1084571/1‑150 (MQ=255)
                            aGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGAAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCAt                                                                                                                       <  2:923163/150‑1 (MQ=255)
                             gCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATc                                                                                                                      <  2:673994/150‑1 (MQ=255)
                              ccTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCg                                                                                                                     >  2:1157593/1‑150 (MQ=255)
                                tGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAAc                                                                                                                                     <  2:1116399/132‑1 (MQ=255)
                                tGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAAc                                                                                                                                     >  1:1116399/1‑132 (MQ=255)
                                tGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGc                                                                                                                   >  2:871286/1‑150 (MQ=255)
                                 gggCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCa                                                                                                                  <  1:800941/150‑1 (MQ=255)
                                            gCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTTGCATCGGCACCAacctgacc                                                                                                       >  1:6643/1‑150 (MQ=255)
                                              cATTGCCCATTACCAGAAACTGGGTATCGACCCGATTACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAacctgacctg                                                                                                     <  1:1077159/150‑1 (MQ=255)
                                              cATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGa                                                                                                                                              >  1:517483/1‑109 (MQ=255)
                                              cATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGa                                                                                                                                              <  2:517483/109‑1 (MQ=255)
                                              cATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAacctgacctg                                                                                                     >  1:277880/1‑150 (MQ=255)
                                                   cccATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCaa                                                                                                                                      <  2:1217468/112‑1 (MQ=255)
                                                   cccATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCaa                                                                                                                                      >  1:1217468/1‑112 (MQ=255)
                                                   cccATTACCAGAAACTAGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACa                                                                                                >  1:475847/1‑150 (MQ=255)
                                                    ccATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACAt                                                                                               >  2:920767/1‑150 (MQ=255)
                                                          ccAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAgg                                                                                         >  2:996809/1‑150 (MQ=255)
                                                               aaCTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTgg                                                                                    <  2:587329/150‑1 (MQ=255)
                                                                 cTGGGTATCGACCTGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGcc                                                                                  <  1:1266596/150‑1 (MQ=255)
                                                                    ggTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGCCCGACGACATCCCAGGTGTGGccccc                                                                               >  2:451323/1‑150 (MQ=255)
                                                                     gTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCa                                                                              <  1:533544/150‑1 (MQ=255)
                                                                          gACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACg                                                                                                                                    >  1:1207816/1‑91 (MQ=255)
                                                                          gACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACg                                                                                                                                    <  2:1207816/91‑1 (MQ=255)
                                                                            cccGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACAt                                                                       >  1:1856719/1‑150 (MQ=255)
                                                                              cGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCg                                                                     >  2:1901793/1‑150 (MQ=255)
                                                                               gATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAacctgacct                                                                                                      <  1:264064/116‑1 (MQ=255)
                                                                               gATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAacctgacct                                                                                                      >  2:264064/1‑116 (MQ=255)
                                                                                aTGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTg                                                                   <  2:1388757/150‑1 (MQ=255)
                                                                                   aCCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTg                                                                   >  2:575514/1‑147 (MQ=255)
                                                                                   aCCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTg                                                                   <  1:575514/147‑1 (MQ=255)
                                                                                    ccAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGTATCAACGTCAGCTTTGGCATCGGCAc                                                                                                                 >  1:635345/1‑100 (MQ=255)
                                                                                    ccAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGTATCAACGTCAGCTTTGGCATCGGCAc                                                                                                                 <  2:635345/100‑1 (MQ=255)
                                                                                     cAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTaa                                                              >  1:1688726/1‑150 (MQ=255)
                                                                                     cAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTaa                                                              >  2:447708/1‑150 (MQ=255)
                                                                                     cAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTaa                                                              >  2:1491635/1‑150 (MQ=255)
                                                                                       aGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCGGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCCTGAACATCGTGCTCaaaa                                                            <  2:182416/150‑1 (MQ=255)
                                                                                          cgatGGTGTTCTTCAACAGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGTCCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGa                                                         <  2:1504876/147‑1 (MQ=255)
                                                                                           gCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGAc                                                        <  2:6643/150‑1 (MQ=255)
                                                                                               gtgtTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGAc                                                    >  1:1719872/1‑150 (MQ=255)
                                                                                                tgtTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGACCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACt                                                   <  1:920767/150‑1 (MQ=255)
                                                                                                 gtTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTg                                                  <  2:1688726/150‑1 (MQ=255)
                                                                                                    ttCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCaa                                               >  1:144582/1‑150 (MQ=255)
                                                                                                       cGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGc                                                                                   >  2:148598/1‑111 (MQ=255)
                                                                                                       cGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGc                                                                                   <  1:148598/111‑1 (MQ=255)
                                                                                                       cGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTGTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACgg                                            >  1:1590101/1‑150 (MQ=255)
                                                                                                         aCGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCATGAACATc                                                                      >  2:770075/1‑121 (MQ=255)
                                                                                                         aCGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCATGAACATc                                                                      <  1:770075/121‑1 (MQ=255)
                                                                                                              tgccgatgccaaagctgacgttgatgcgTCCGCGCAGGGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTg                                                  <  1:1465284‑M1/109‑1 (MQ=255)
                                                                                                              tgccgatgccaaagctgacgttgatgcgTCCGCGCAGGGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTg                                                  >  2:1465284‑M1/29‑137 (MQ=255)
                                                                                                              cTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCg                                     >  1:458679/1‑150 (MQ=255)
                                                                                                               tCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTc                                                                                                                                  <  1:225022/56‑1 (MQ=255)
                                                                                                               tCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTc                                                                                                                                  >  2:225022/1‑56 (MQ=255)
                                                                                                                cAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGt                                   <  1:1157593/150‑1 (MQ=255)
                                                                                                                 aaCCTCCCCCGCTCGCTGGGAATTGTCCGCGCCCTGCGCGGACGCATCAACGTACGCTTTGGTCTAGCAAAGCACGGGACCTGCGACATCCCCGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTg                                  <  1:1304608/150‑1 (MQ=255)
                                                                                                                    cTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCa                                                                                           <  1:109484/90‑1 (MQ=255)
                                                                                                                    cTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCa                                                                                           >  2:109484/1‑90 (MQ=255)
                                                                                                                         ccGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGAt                          >  2:1931033/1‑150 (MQ=255)
                                                                                                                           gctcgctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGAtct                        >  2:721141/1‑150 (MQ=255)
                                                                                                                             tcgctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATctct                      <  1:1569415/150‑1 (MQ=255)
                                                                                                                              cgctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTg                     <  2:144582/150‑1 (MQ=255)
                                                                                                                               gctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGAtctc                       <  1:689323/147‑1 (MQ=255)
                                                                                                                               gctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGAtctc                       >  2:689323/1‑147 (MQ=255)
                                                                                                                               acgttgatgcgTCCGCGCAGGGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAacctgac                                                                                                        >  1:886493‑M1/12‑66 (MQ=255)
                                                                                                                               acgttgatgcgTCCGCGCAGGGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAacctgac                                                                                                        <  2:886493‑M1/55‑1 (MQ=255)
                                                                                                                                ctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGAc                   <  2:277880/150‑1 (MQ=255)
                                                                                                                                ctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGAc                   <  2:475847/150‑1 (MQ=255)
                                                                                                                                   gAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGaa                <  2:458679/150‑1 (MQ=255)
                                                                                                                                      tgcgTCCGCGCAGGGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGa                    <  2:1820471‑M1/139‑1 (MQ=255)
                                                                                                                                      tgcgTCCGCGCAGGGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGa                    >  1:1820471‑M1/5‑143 (MQ=255)
                                                                                                                                       ttttCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAgccg            >  2:1326720/1‑150 (MQ=255)
                                                                                                                                         ttCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCAACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAgccgcc          <  1:1901793/150‑1 (MQ=255)
                                                                                                                                         ttCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACAGCCTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAgccgcc          <  1:451323/150‑1 (MQ=255)
                                                                                                                                               GCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAGCCGCCAAGAcc    >  1:1931033‑M1/1‑150 (MQ=255)
                                                                                                                                                 ccTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACCTCGTGCTTAAAATGACCGACTGCACCGGCGCGCCGGTGGCCAAGATCTCTGACGAAGCCGCCAAGACCCa  <  2:773355/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAGCCGCCAAGACCCA  >  NC_002947/5535882‑5536176

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