Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,670,018 T→G 16.1% intergenic (+71/‑96) prpF → / → prpD aconitate isomerase/2‑methylcitrate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,670,0180TG16.1% 147.2 / 21.3 67intergenic (+71/‑96)prpF/prpDaconitate isomerase/2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (26/30);  new base G (5/6);  total (31/36)
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.26e-01

AACCAAGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCT  >  NC_002947/2669872‑2670163
                                                                                                                                                  |                                                                                                                                                  
aaCCAAGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAAt                                                                                                                                                <  1:1395218/150‑1 (MQ=255)
  ccaccgCAATCACGAGCCGTAGGGCTCGCGTGCTGATGGAGGGCGGGGTTCGCCTGCCCGGCGATAGCTTCTAAAGCACCACGGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATca                                                                                                                                              <  2:1465823/145‑1 (MQ=38)
   cAAGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATcac                                                                                                                                             >  2:2402507/1‑150 (MQ=255)
    aaGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATcaca                                                                                                                                            >  1:3829372/1‑150 (MQ=255)
    aaGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATcaca                                                                                                                                            >  1:3829453/1‑150 (MQ=255)
          aTCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGttt                                                                                                                                      <  2:55666/150‑1 (MQ=255)
             aTGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGaa                                                                                                                                   <  1:1181334/150‑1 (MQ=255)
                  cggttgTGCTGGGGTGCTGTTGTTGGGGTGGGTTGGCGTGCCTGGCTTTAGCTTCTAAAGCACCACTGGCGCTGGGCTTAAATAGGGGCCGATTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtg                                                                                                                              <  2:3892444/145‑1 (MQ=255)
                   cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                             <  1:567818/150‑1 (MQ=255)
                   cGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgt                                                                                                                             <  1:923701/150‑1 (MQ=255)
                    gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAGAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGGTAATCACAGGGTTTGAATGCggtg                                                                                                                            >  2:1218957/1‑150 (MQ=255)
                    gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                            <  2:2174848/150‑1 (MQ=255)
                    gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                            <  2:3786005/150‑1 (MQ=255)
                    gTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCtgtg                                                                                                                            <  2:3786038/150‑1 (MQ=255)
                           tCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATccc                                                                                                                     <  2:4184231/150‑1 (MQ=255)
                            cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                    <  2:3967605/150‑1 (MQ=255)
                            cgcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCACGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCt                                                                                                                    <  1:2969328/150‑1 (MQ=255)
                             gcgTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtg                                                                                                                   <  2:4945747/150‑1 (MQ=255)
                               gTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCtgtg                                                                                                                 <  1:1036396/150‑1 (MQ=255)
                                tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAAACACAGGGTTTGAATGCTGTGATATCCCTGTgg                                                                                                                >  2:672740/1‑150 (MQ=255)
                                tGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTGTGCAGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTgg                                                                                                                <  1:1454793/150‑1 (MQ=17)
                                tGCTGATGGAGGGCTGGATTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTGTGCAGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTgg                                                                                                                <  1:1454883/150‑1 (MQ=17)
                                  cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCTGCTGTGCAGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                              <  1:4143176/150‑1 (MQ=17)
                                  cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGa                                                                                                              <  2:3236995/150‑1 (MQ=255)
                                  cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGa                                                                                                              <  2:1843656/150‑1 (MQ=255)
                                  cTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGa                                                                                                              <  2:3785186/150‑1 (MQ=18)
                                   tGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGGGTGGCCCCGCCCCCACGTAAAAACACGGGGTGTATTCGGTGATATAACGGTggggg                                                                                                             >  1:3004068/1‑148 (MQ=255)
                                    gATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGc                                                                                                            <  2:579754/150‑1 (MQ=18)
                                      tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                          >  2:1411470/1‑150 (MQ=255)
                                      tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCgg                                                                                                          <  1:1923985/150‑1 (MQ=255)
                                      tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCgg                                                                                                          <  2:643933/150‑1 (MQ=18)
                                      tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCGCCCACACGGAAAAACAAGGTTTGAATGCGGGGGTATCACAGGGGGAggcg                                                                                                          >  1:4764847/1‑147 (MQ=255)
                                      tGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCGTGCCCGCCCCCCCCGGGAAAACCAGGGGTTGTATGCTGTGTTATCCCGGTGGGGGCgg                                                                                                          >  1:665987/1‑150 (MQ=255)
                                          gggCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCaa                                                                                                      <  2:254573/150‑1 (MQ=255)
                                          gggCTGGGTTCGCGTGCCTGGCGAGAGCTTCTACCGCACCACTGGCGCTCGGCGTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCaa                                                                                                      <  2:1537744/150‑1 (MQ=39)
                                              tGGGTTCGCGTGCCTGGCGATAGCTTCTACAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGccc                                                                                                  <  2:4714332/150‑1 (MQ=255)
                                              tGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCGCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATACCGGTGGGAGCGGGCAAGccc                                                                                                  >  2:1462491/1‑150 (MQ=255)
                                                     gcgTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCGCCCACAGGTAATCCCCGGGGTTGTATGCTGTGATATACCTGTGGGAGCGGGCAAGCCCGCGGAgg                                                                                           >  1:4544343/1‑150 (MQ=38)
                                                                 aTAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCCCTCCCACAGGTAATCACCGGGTTTTTATGCTGTGCTATCCCTGGGGGTGCGGGCAAGCCCGCGAAGGGCTGCAcagccg                                                                               >  2:889539/1‑148 (MQ=38)
                                                                      ttCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCGCCCACAGGTAATCACAGGGGTTGGATTCTGTGATATACCGGTGGGAGCGGGCCAGCCCGCGAAGGGCGGCACAGCAGCCCCg                                                                          >  1:1155064/1‑150 (MQ=37)
                                                                       tCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACCGGGTTGGAAAGCTGTGATAACCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCgg                                                                         >  1:5117480/1‑150 (MQ=255)
                                                                            aGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCCCTCCCACACGTAAAAAAAGGGTTTGGATTCTGTGATATACCGGTGGGAGCGGGCAAGCCCGCCAATGGCTGCACAGCAGCCCCCGccaat                                                                    >  1:1197104/1‑146 (MQ=37)
                                                                                 caggggCGCTCGGCTTCATGAGGGGCCGGTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGGTAATCACAGGTTTTTATTGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGaa                                                               <  2:3004068/146‑1 (MQ=37)
                                                                                   cTGGCGCTCGTATTAAAGAGGGGCCGATTTCCGCCAATTCGCGGGCTTGCCCGATCCCACTGTTAAAAACTGTTTTTTAATGATTTTATATCCCTTTGGGAGCGGCAAACCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAgg                                                             <  1:370960/150‑1 (MQ=255)
                                                                                     ggCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACTGGTAATCACTGGTTTTTATTGCTGTTATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGa                                                           <  2:2184877/150‑1 (MQ=38)
                                                                                      gcgcggggCTTCCAGCGGGGCCGCTTGGCGGCCCTCCGCGGGCTTGCCCGCTCCCACCGGTAAGCCCAGGTTTTGATTGCTGTGATATCCGTGGGGGAGCGGGCACGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAg                                                          <  2:3996693/144‑1 (MQ=37)
                                                                                        gcTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGt                                                                                                                                        <  1:33689/70‑1 (MQ=255)
                                                                                        gcTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGt                                                                                                                                        >  2:33689/1‑70 (MQ=255)
                                                                                               ttCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGGTTGAATGCTGTGATATCCCGGTGGGAGCGGGCAAGCCCGCCAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATg                                                 >  2:4811106/1‑150 (MQ=255)
                                                                                                      ggggCCGCTTTGCGGCCCTTCGCGGGCTTGCCCCCTCCCCCAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGCGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCCCATCTGAAGGGGGCAGGATAAGAGCGCCa                                          >  2:4386015/1‑150 (MQ=25)
                                                                                                              tttGCGGCCCTTCGCGGGCTTGCCCGCTCCCACCGGTAAACACAGGGTTTGAATGCTGCGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGAc                                  <  1:795108/150‑1 (MQ=255)
                                                                                                                            gggggtggCCCGCTCCCACAGGTAATCACAGGTTTTGAAGGCTGTGATAGCCCGGTGGGTGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGccc                    <  1:4131071/143‑1 (MQ=255)
                                                                                                                                 ttGCCCGCTCCCACCGGGATATCACAGCATTCCAACCCTGTGAGTACCGGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGACGGGAGCTGGATATGAGCGCCCACGTAGACCTCCACGATCGCCCGGagg               >  2:4630581/1‑148 (MQ=17)
                                                                                                                                 ttGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt               >  1:3767301/1‑150 (MQ=255)
                                                                                                                                 ttGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt               >  1:1921356/1‑150 (MQ=255)
                                                                                                                                 ttGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt               >  2:1219243/1‑150 (MQ=255)
                                                                                                                                 ttGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGAtt               >  1:2773215/1‑150 (MQ=255)
                                                                                                                                    cccGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACg            >  1:935063/1‑150 (MQ=255)
                                                                                                                                     ccGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa           >  1:765637/1‑150 (MQ=18)
                                                                                                                                     ccGCTCCCACAGGGATATCACAGCATTCAAACCCTGTGATTACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa           >  1:1547585/1‑150 (MQ=18)
                                                                                                                                     ccGCGCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGa           >  1:4528213/1‑150 (MQ=255)
                                                                                                                                      cGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGAc          >  1:1936115/1‑150 (MQ=255)
                                                                                                                                      cGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGAc          >  2:2443164/1‑150 (MQ=255)
                                                                                                                                       gCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGAcc         >  1:2977806/1‑150 (MQ=255)
                                                                                                                                           ccACAGGTAATCACAGGGTTTGAATGCTGAGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGt     >  1:3772365/1‑150 (MQ=255)
                                                                                                                                              ccggTAACCCCCGGGTTTTAATGCTTTGATATCCCTGTGGGAGCGGGAAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCt  <  1:1644361/148‑1 (MQ=255)
                                                                                                                                              ccggTAAACACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCTAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCt  <  1:1644459/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AACCAAGGCAATCATGAGCCGTAGTGCTCGCGTGCTGATGGAGGGCTGGGTTCGCGTGCCTGGCGATAGCTTCTAAAGCACCACTGGCGCTCGGCTTCAAGAGGGGCCGCTTTGCGGCCCTTCGCGGGCTTGCCCGCTCCCACAGGTAATCACAGGGTTTGAATGCTGTGATATCCCTGTGGGAGCGGGCAAGCCCGCGAAGGGCTGCACAGCAGCCCCGGCAATTCTGAAGGGAGCTGGATATGAGCGCCAACGTAGACCTCAACGATCGCCCGGATTACGACCGGGTGCT  >  NC_002947/2669872‑2670163

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