Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,851,577 G→T 9.1% W207C (TGG→TGT relA → ATP:GTP 3'‑pyrophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,851,5770GT 9.1% 399.3 / 19.6 132W207C (TGG→TGTrelAATP:GTP 3'‑pyrophosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (51/69);  new base T (6/6);  total (57/75)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

Only 100 of 137 total aligned reads displayed.

CGCGTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAA  >  NC_002947/1851428‑1851726
                                                                                                                                                     |                                                                                                                                                     
cgcgTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTgg                                                                                                                                                       >  1:2658725/1‑150 (MQ=255)
cgcgTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCATGGTTATTTTGCACAGCAAGTgg                                                                                                                                                       <  2:2283328/150‑1 (MQ=255)
       cgcTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTgg                                                                                                                                                <  2:184884/150‑1 (MQ=255)
       cgcTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTgg                                                                                                                                                >  1:766496/1‑150 (MQ=255)
        gcTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGa                                                                                                                                               >  1:2136340/1‑150 (MQ=255)
          tGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTTAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAg                                                                                                                                             <  1:1073844/150‑1 (MQ=255)
          tGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCACGTGGGAGCTGGAAg                                                                                                                                             >  1:2636343/1‑150 (MQ=255)
            aTCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGAc                                                                                                                                           <  2:1604511/150‑1 (MQ=255)
            aTCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGAc                                                                                                                                           <  1:1072790/150‑1 (MQ=255)
            aTCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGAc                                                                                                                                           >  1:2578182/1‑150 (MQ=255)
            aTCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCAGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGGATCGATCACATCAAGTGGGAGCTGGAAGAc                                                                                                                                           <  2:894657/150‑1 (MQ=255)
             tCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGCGGGAGCGGGAAGAcc                                                                                                                                          >  2:2859171/1‑150 (MQ=255)
              cAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCt                                                                                                                                         <  2:1046873/150‑1 (MQ=255)
              cAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCt                                                                                                                                         >  2:921774/1‑150 (MQ=255)
              cAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCAGCGGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCt                                                                                                                                         >  2:178915/1‑150 (MQ=255)
                aGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGt                                                                                                                                       <  2:2929603/150‑1 (MQ=255)
                 gCTGGCCGAACGCACGTGCTCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTc                                                                                                                                      <  1:921774/150‑1 (MQ=255)
                     gCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTc                                                                                                                                  <  2:1515238/150‑1 (MQ=255)
                      ccGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTcc                                                                                                                                 <  2:1962181/150‑1 (MQ=255)
                      ccGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTcc                                                                                                                                 <  2:1698740/150‑1 (MQ=255)
                            gCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTAcc                                                                                                                           >  1:2394264/1‑150 (MQ=255)
                             cACGTGCGCAATCCGGGCGGTCAAGTCCGCAGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGTGAGCTGGAAGACCTGTCCTTCCGCTACCt                                                                                                                          <  2:1777988/150‑1 (MQ=255)
                              aCGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTg                                                                                                                         >  1:1090344/1‑150 (MQ=255)
                                 tGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGt                                                                                                                                       >  2:618811/1‑133 (MQ=255)
                                 tGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGt                                                                                                                                       <  1:618811/133‑1 (MQ=255)
                                    gcAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGTGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGccc                                                                                                                   >  1:2863156/1‑150 (MQ=255)
                                    gcAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGccc                                                                                                                   <  2:1090344/150‑1 (MQ=255)
                                    gcAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGccc                                                                                                                   <  2:2658725/150‑1 (MQ=255)
                                      aaTCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGa                                                                                                                 <  2:924432/150‑1 (MQ=255)
                                       aTCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGAc                                                                                                                <  1:2117530/150‑1 (MQ=255)
                                        tCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTgg                                                                                                                                                >  2:1489924/1‑117 (MQ=255)
                                         ccGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCa                                                                                                              <  1:2737853/150‑1 (MQ=255)
                                                gTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAg                                                                                                       <  2:2136340/150‑1 (MQ=255)
                                                 tCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGc                                                                                                      >  1:2435270/1‑150 (MQ=255)
                                                   aaGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTggg                                                                                                                                                      <  1:2728246/100‑1 (MQ=255)
                                                   aaGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTggg                                                                                                                                                      >  2:2728246/1‑100 (MQ=255)
                                                   aaGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAg                                                                                                    >  2:2707379/1‑150 (MQ=255)
                                                          ccgATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCa                                                                                             >  1:2536262/1‑150 (MQ=255)
                                                          ccgATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCa                                                                                             >  1:1152411/1‑150 (MQ=255)
                                                               gACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTg                                                                                        <  1:441256/150‑1 (MQ=255)
                                                                 cGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGtt                                                                                      <  1:2707379/150‑1 (MQ=255)
                                                                 cGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGtt                                                                                      <  2:2435270/150‑1 (MQ=255)
                                                                   agaAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGc                                                                                    <  2:2536262/150‑1 (MQ=255)
                                                                         gCCTGCGTGTCGCGCGTGAAGTGCTCGACATCTACGCGCCGCTCGCGCACCGTCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAgc                                                                              <  2:2093183/150‑1 (MQ=255)
                                                                           cTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAgcgg                                                                            >  2:1188561/1‑150 (MQ=255)
                                                                            tGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAgcggc                                                                           <  1:1376822/150‑1 (MQ=255)
                                                                                 gtCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGCGTATCGATCACATCAAGTGGTAGCTGGAAGACCTGTCCTTCCACTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTg                                                                      <  1:523545/150‑1 (MQ=255)
                                                                                   cgcgcgTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGa                                                                    <  1:1392425/150‑1 (MQ=255)
                                                                                    gcgcgTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGAc                                                                   >  1:87298/1‑150 (MQ=255)
                                                                                     cgcgTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGAcc                                                                  <  2:1260127/150‑1 (MQ=255)
                                                                                      gcgTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATc                                                                                                 <  1:571279/118‑1 (MQ=255)
                                                                                      gcgTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATc                                                                                                 >  2:571279/1‑118 (MQ=255)
                                                                                          gAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGc                                                                                    <  2:2572821/127‑1 (MQ=255)
                                                                                          gAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGc                                                                                    >  1:2572821/1‑127 (MQ=255)
                                                                                          gAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAg                                                             >  2:2832180/1‑150 (MQ=255)
                                                                                          gAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCGACCTGGAGCCCGACCAGTACAAGCAGATCGCCAACCGTGTGCATGAGCGGCGGCTGGACCGCGAg                                                             <  1:1331430/150‑1 (MQ=255)
                                                                                           aaGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAgc                                                            <  1:1687543/150‑1 (MQ=255)
                                                                                            aGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGTGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAgcg                                                           >  1:2187009/1‑150 (MQ=255)
                                                                                             gtgtTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGTGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGtt                                                                                      <  1:1218693/122‑1 (MQ=255)
                                                                                             gtgtTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGTGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGtt                                                                                      >  2:1218693/1‑122 (MQ=255)
                                                                                               gtTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCtt                                                        <  1:1188561/150‑1 (MQ=255)
                                                                                                ttCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGTGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTtc                                                       >  2:841959/1‑150 (MQ=255)
                                                                                                ttCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAgcgc                                                          >  2:2112907/1‑147 (MQ=255)
                                                                                                ttCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAgcgc                                                          <  1:2112907/147‑1 (MQ=255)
                                                                                                ttCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCt                                                         <  2:3008423/148‑1 (MQ=255)
                                                                                                ttCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCt                                                         >  1:3008423/1‑148 (MQ=255)
                                                                                                ttCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTtc                                                       >  2:290782/1‑150 (MQ=255)
                                                                                                ttCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTtc                                                       >  1:1784663/1‑150 (MQ=255)
                                                                                                ttCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGATGTTACATGAGCGGCGGCTGGACCCCGAGCGCTtc                                                       >  2:388975/1‑150 (MQ=255)
                                                                                                   gACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGTGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTtcatc                                                    >  1:3029622/1‑150 (MQ=255)
                                                                                                        cTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCCGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGa                                               >  2:495926/1‑150 (MQ=255)
                                                                                                              gccgcTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGtgat                                         <  1:976298/150‑1 (MQ=255)
                                                                                                              gccgcTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGGTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGtgat                                         <  1:2090379/150‑1 (MQ=255)
                                                                                                               ccgcTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGtgatg                                        <  2:87298/150‑1 (MQ=255)
                                                                                                                  cTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGACCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAAc                                     <  1:904525/150‑1 (MQ=255)
                                                                                                                     gcgcACCGCCTGGGTATCGATCACATCAAGTGTGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAg                                  <  2:2863156/150‑1 (MQ=255)
                                                                                                                       gcACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCt                                >  2:593678/1‑150 (MQ=255)
                                                                                                                        cACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTg                               >  1:2898931/1‑150 (MQ=255)
                                                                                                                             ccTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGtgat                                         >  2:1482024/1‑135 (MQ=255)
                                                                                                                             ccTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGtgat                                         <  1:1482024/135‑1 (MQ=255)
                                                                                                                                gggTATCGATCACATCAAGTGGGAGGTCGAACACCTGCCCTTCCGCTACCTGGAGCGCGACCAGTACAAGCAGTTCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATCAACCAGCTGCAGAACGa                       <  2:232359/150‑1 (MQ=255)
                                                                                                                                 ggTATCGATCACATCAAGTGTGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCGCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAg                      <  1:1428811/150‑1 (MQ=255)
                                                                                                                                   tATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGtt                    >  1:1503571/1‑150 (MQ=255)
                                                                                                                                       gATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCaa                                                                                            >  1:2209023/1‑74 (MQ=255)
                                                                                                                                       gATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCaa                                                                                            <  2:2209023/74‑1 (MQ=255)
                                                                                                                                       gATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCggcggc                                                                        <  1:1685520/94‑1 (MQ=255)
                                                                                                                                       gATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCggcggc                                                                        >  2:1685520/1‑94 (MQ=255)
                                                                                                                                       gATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAg                      >  1:2250450/1‑144 (MQ=255)
                                                                                                                                       gATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAg                      <  2:2250450/144‑1 (MQ=255)
                                                                                                                                        aTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCggcggc                                                                        >  2:2691531/1‑93 (MQ=255)
                                                                                                                                        aTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTTCAGAACGAGTTGCTCg               >  2:2331385/1‑150 (MQ=255)
                                                                                                                                        aTCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCg               >  1:263262/1‑150 (MQ=255)
                                                                                                                                         tCACATCAAGTGTGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGAc                                                                                                                <  2:1287765/52‑1 (MQ=255)
                                                                                                                                             aTCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACt          >  1:2908161/1‑150 (MQ=255)
                                                                                                                                                aaGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGc       <  2:263262/150‑1 (MQ=255)
                                                                                                                                                 aGTGGGATCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCg      <  1:178915/150‑1 (MQ=255)
                                                                                                                                                 aGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAAc                                     <  2:1777447/119‑1 (MQ=255)
                                                                                                                                                 aGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAAc                                     >  1:1777447/1‑119 (MQ=255)
                                                                                                                                                  gTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGt     >  1:2061344/1‑150 (MQ=255)
                                                                                                                                                     ggAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGaa  >  2:1137119/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGTGGCGCTGATCAAGCTGGCCGAACGCACGTGCGCAATCCGGGCGGTCAAGTCCGCCGATGACGAGAAACGCCTGCGTGTCGCGCGTGAAGTGTTCGACATCTACGCGCCGCTCGCGCACCGCCTGGGTATCGATCACATCAAGTGGGAGCTGGAAGACCTGTCCTTCCGCTACCTGGAGCCCGACCAGTACAAGCAGATCGCCAAGCTGTTGCATGAGCGGCGGCTGGACCGCGAGCGCTTCATCAGCGACGTGATGAACCAGCTGCAGAACGAGTTGCTCGCCACTGGCGTGAA  >  NC_002947/1851428‑1851726

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