Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,851,755 A→G I267V (ATC→GTC)  relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,851,7550AG96.7% 293.3 / ‑1.4 92I267V (ATC→GTC) relAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (42/46);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

AGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTC  >  NC_002947/1851609‑1851899
                                                                                                                                                   |                                                                                                                                                
aGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCt                                                                                                                                               <  2:1663731/150‑1 (MQ=255)
   ccGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCAGACATCACCGGGCGGGGGAAACATGTCTAtt                                                                                                                                            >  2:2235040/1‑150 (MQ=255)
          gTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTg                                                                                                                                     >  1:1995972/1‑150 (MQ=255)
            aCAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGCGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGgc                                                                                                                                   <  2:1397294/150‑1 (MQ=255)
              aaGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGgcgc                                                                                                                                 <  2:1995972/150‑1 (MQ=255)
               cgCAGATCGCCAAGATCTTCCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCa                                                                                                                                >  2:2084731/2‑150 (MQ=255)
                gCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGAACGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCaa                                                                                                                               >  1:969968/1‑150 (MQ=255)
                   gATCGCCAAGCTGTTGCATGAGCTGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGAt                                                                                                                            >  1:1300783/1‑150 (MQ=255)
                       gCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAg                                                                                                                        >  2:1718798/1‑150 (MQ=255)
                       gCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTCGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCCAGATGCAg                                                                                                                        >  2:416602/1‑150 (MQ=255)
                         cAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAgcg                                                                                                                      >  1:48261/1‑150 (MQ=255)
                             cTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCaaa                                                                                                                  <  2:535336/150‑1 (MQ=255)
                              tGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAg                                                                                                                 <  2:969968/150‑1 (MQ=255)
                                ttGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTTATGAACCAGCTGCAGAACGAGTTGCTCGCCACGGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGc                                                                                                               <  2:413068/150‑1 (MQ=255)
                                  gCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCTATCTGGCGCAAGATGCAGCGCAAAGGCCt                                                                                                             >  1:452039/1‑150 (MQ=255)
                                   cATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTg                                                                                                            <  2:48261/150‑1 (MQ=255)
                                    aTGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAg                                                                                                                        <  1:549097/137‑1 (MQ=255)
                                    aTGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAg                                                                                                                        >  2:549097/1‑137 (MQ=255)
                                     tGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCTATCTGGCGCAAGATGCAGCGCAAAGGCCTGGa                                                                                                          <  2:452039/150‑1 (MQ=255)
                                      gAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAACCCTGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAg                                                                                                         <  1:2235040/150‑1 (MQ=255)
                                       aGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGt                                                                                                        <  1:1718798/150‑1 (MQ=255)
                                        gcggcgGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGtt                                                                                                       >  1:1887203/1‑150 (MQ=255)
                                        gcggcgGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGtt                                                                                                       <  2:2143642/150‑1 (MQ=255)
                                        gcggcgGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGtt                                                                                                       <  1:1505709/150‑1 (MQ=255)
                                              gcTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGACGGCATGCTAGTCGCTCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCa                                                                                                 >  2:1772673/1‑150 (MQ=255)
                                                tGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGa                                                                                               >  2:1299172/1‑150 (MQ=255)
                                                    ccGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTa                                                                                           >  2:262199/1‑150 (MQ=255)
                                                        gAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTacgac                                                                                       >  2:381351/1‑150 (MQ=255)
                                                           cgcTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGAcgtg                                                                                    <  1:712279/150‑1 (MQ=255)
                                                              ttCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGAcgtgcgt                                                                                 >  2:803079/1‑150 (MQ=255)
                                                              ttCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGAcgtgcgt                                                                                 >  2:659123/1‑150 (MQ=255)
                                                                catcaGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACgtgcgtgc                                                                               >  2:76024/1‑150 (MQ=255)
                                                                  tcaGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGAcgtgcgt                                                                                 <  2:76566/146‑1 (MQ=255)
                                                                  tcaGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGAcgtgcgt                                                                                 >  1:76566/1‑146 (MQ=255)
                                                                      cGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCgtgcggtgcg                                                                         >  1:1468224/1‑150 (MQ=255)
                                                                        aCGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGAGGTGcgcg                                                                       >  1:1785631/1‑150 (MQ=255)
                                                                        aCGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGGGCGTTGcgcg                                                                       >  1:1834199/1‑150 (MQ=255)
                                                                           tgatgaACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGc                                                                    <  2:1887203/150‑1 (MQ=255)
                                                                           tgatgaACCAGCTGAAGCACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAGGGTGATGGATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGc                                                                    <  1:396025/150‑1 (MQ=255)
                                                                            gatgaACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCACAAAGGCCTGGAGTTCAGCCAGATCTacgac                                                                                       >  1:516770/1‑130 (MQ=255)
                                                                            gatgaACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCACAAAGGCCTGGAGTTCAGCCAGATCTacgac                                                                                       <  2:516770/130‑1 (MQ=255)
                                                                               gaACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTg                                                                  <  1:38069/148‑1 (MQ=255)
                                                                               gaACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTg                                                                  >  2:38069/1‑148 (MQ=255)
                                                                               gaACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGt                                                                <  2:1300783/150‑1 (MQ=255)
                                                                                aaCCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTg                                                               >  1:534486/1‑150 (MQ=255)
                                                                                   cAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTgcgc                                                                        <  1:251167/138‑1 (MQ=255)
                                                                                   cAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTgcgc                                                                        >  2:251167/1‑138 (MQ=255)
                                                                                      cTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGaa                                                         <  2:1304504/150‑1 (MQ=255)
                                                                                       tGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGTAGGGCAATGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGaaa                                                        <  2:441239/150‑1 (MQ=255)
                                                                                         cAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATc                                                      <  1:2084731/150‑1 (MQ=255)
                                                                                           gagcGAGGTGCGCCCGACCGGCGGGAAGGCCCACATCAGCGGCCGGGCGAAACATGTCGATTCGATCTGGCGCAAGATGCCGCGCAAAGGCCGGGAGTTCAGCCAGATCTACGACCTTCGTGCGGTGCGCGTGCTGGTGCCGGAAATCcg                                                    <  1:2144598/147‑1 (MQ=255)
                                                                                           gAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCa                                                                                                                         >  1:765202/1‑81 (MQ=255)
                                                                                           gAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCa                                                                                                                         <  2:765202/81‑1 (MQ=255)
                                                                                              cGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGa                                                 <  1:76024/150‑1 (MQ=255)
                                                                                                 gTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTg                                              <  1:1299172/150‑1 (MQ=255)
                                                                                                    gCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTa                                           <  2:534486/150‑1 (MQ=255)
                                                                                                     cTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTac                                          >  2:1675333/1‑150 (MQ=255)
                                                                                                        gCCACTGGCGTGAAGGCCGACATCAGCGGACGGGCGAAACATGTCTATTCGATCTGGCGCAAGATACAGCGCAAAGGCCTGGAGCTCACCCAGATCTACGACGTTCGTGCGGTGCGCGTGCGGGTGCCGTCAATCCGCGACTGCTACGcc                                       >  2:399486/1‑150 (MQ=255)
                                                                                                          cACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGc                                                                    >  2:1763461/1‑119 (MQ=255)
                                                                                                          cACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGc                                                                    <  1:1763461/119‑1 (MQ=255)
                                                                                                          cACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGAGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACcgc                                     >  1:408121/1‑150 (MQ=255)
                                                                                                                  tGAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCa                             <  1:381351/150‑1 (MQ=255)
                                                                                                                   gAAGGCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCAt                            <  1:1974917/150‑1 (MQ=255)
                                                                                                                    aaGGCCGTCATCAGCCGCCGGGCGAAACAGGCCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCAGGAGGTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATc                           <  1:617847/150‑1 (MQ=255)
                                                                                                                       gCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTgg                                                                                                           <  2:140545/67‑1 (MQ=255)
                                                                                                                       gCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTgg                                                                                                           >  1:140545/1‑67 (MQ=255)
                                                                                                                       gCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACcgc                                     >  1:106142/1‑137 (MQ=255)
                                                                                                                       gCCGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACcgc                                     <  2:106142/137‑1 (MQ=255)
                                                                                                                        ccGGCAACAGCGGGCGGGGGAAACATGTCTATTCGATCTGGCGGAAGATGCAGGGCCAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGc                       <  2:1834199/150‑1 (MQ=255)
                                                                                                                        ccGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGc                       >  1:1731670/1‑150 (MQ=255)
                                                                                                                        ccGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGc                       <  1:1232781/150‑1 (MQ=255)
                                                                                                                         cGACATCAGTGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGca                      <  2:408121/150‑1 (MQ=255)
                                                                                                                         cGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCgtgcggtgcg                                                                         >  2:1652424/1‑99 (MQ=255)
                                                                                                                         cGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCgtgcggtgcg                                                                         <  1:1652424/99‑1 (MQ=255)
                                                                                                                         cGACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTCCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGca                      >  2:1515451/1‑150 (MQ=255)
                                                                                                                          gACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCt                                                                                            <  1:250444/79‑1 (MQ=255)
                                                                                                                          gACATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCt                                                                                            >  2:250444/1‑79 (MQ=255)
                                                                                                                           aCATCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGcaca                    <  1:659123/150‑1 (MQ=255)
                                                                                                                             aTCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACAcc                  >  2:1767651/1‑150 (MQ=255)
                                                                                                                              tCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCt                 <  1:262199/150‑1 (MQ=255)
                                                                                                                              tCAGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCt                 >  1:1487151/1‑150 (MQ=255)
                                                                                                                                aGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTgg                                <  2:1431696/133‑1 (MQ=255)
                                                                                                                                aGCGGCCGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTgg                                >  1:1431696/1‑133 (MQ=255)
                                                                                                                                      cGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTggcgg         <  1:2229586/150‑1 (MQ=255)
                                                                                                                                      cGGGCGAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTggcgg         >  2:1858542/1‑150 (MQ=255)
                                                                                                                                           gAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCatat    >  2:421365/1‑150 (MQ=255)
                                                                                                                                           gAAACATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCatat    >  2:421020/1‑150 (MQ=255)
                                                                                                                                              aaCATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGa                                                 <  1:1380773/103‑1 (MQ=255)
                                                                                                                                              aaCATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGa                                                 >  2:1380773/1‑103 (MQ=255)
                                                                                                                                              aaCATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTc  <  2:1731670/150‑1 (MQ=255)
                                                                                                                                              aaCATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTc  <  1:421365/150‑1 (MQ=255)
                                                                                                                                              aaCATGTCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTc  <  1:421020/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
AGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAAGGCCGACATCAGCGGCCGGGCGAAACATATCTATTCGATCTGGCGCAAGATGCAGCGCAAAGGCCTGGAGTTCAGCCAGATCTACGACGTGCGTGCGGTGCGCGTGCTGGTGCCGGAAATCCGCGACTGCTACACCGCGCTGGGCATCGTGCACACCTTGTGGCGGCATATTC  >  NC_002947/1851609‑1851899

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