Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,614,754 A→C 7.9% I273M (ATT→ATG hldE ← bifunctional heptose 7‑phosphate kinase/heptose 1‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,614,7540AC 7.9% 174.7 / 7.1 69I273M (ATT→ATGhldEbifunctional heptose 7‑phosphate kinase/heptose 1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base C (4/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

CCAGCACGCCGCGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTGACCAGCAAGGCGCCC  >  NC_002947/5614606‑5614901
                                                                                                                                                    |                                                                                                                                                   
ccAGCACGCCGCGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCaa                                                                                                                                                    >  2:546740/1‑150 (MQ=255)
 cAGCACGCCGCGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAAt                                                                                                                                                   <  1:84668/150‑1 (MQ=255)
  aGCACGCCGCGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATg                                                                                                                                                  >  2:955840/1‑150 (MQ=255)
  aGCACGCCGCGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATg                                                                                                                                                  <  1:546740/150‑1 (MQ=255)
     aCGCCGCGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGcc                                                                                                                                               >  2:830291/1‑150 (MQ=255)
       gccgcGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCcgc                                                                                                                                             >  1:870773/1‑150 (MQ=255)
       gccgcGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCcgc                                                                                                                                             >  1:70302/1‑150 (MQ=255)
            gTTCGGAGACTTCCTCACGCTGGATGGCACGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTGGCCGAACACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCCATGGCCGCGGCaa                                                                                                                                        >  2:376938/1‑150 (MQ=255)
             ctCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAg                                                                                                                                       <  1:242376/149‑1 (MQ=255)
                ggAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTCGCCCGCGCCACCGCATGGGGCAGTGCCTCACCCGCGGCAATGGCCGCGGCAAggag                                                                                                                                    >  2:242376/1‑148 (MQ=255)
                ggAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGt                                                                                                                                    <  2:408392/150‑1 (MQ=255)
                    ccTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGaa                                                                                                                                >  2:131937/1‑150 (MQ=255)
                      ttCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAAt                                                                                                                              <  1:131937/150‑1 (MQ=255)
                      ttCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAAt                                                                                                                              <  2:801767/150‑1 (MQ=255)
                        ccTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCGGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGa                                                                                                                            <  2:690236/150‑1 (MQ=255)
                        ccTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGa                                                                                                                            <  1:465521/150‑1 (MQ=255)
                        ccTCACGCTGGATGGCCCGGCGCAGACCCGGGGCGACGCCTTGTCTCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGa                                                                                                                            <  1:804870/150‑1 (MQ=255)
                          tCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCCATTGCCGCgg                                                                                                                                           >  1:893548/1‑133 (MQ=255)
                          tCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCCATTGCCGCgg                                                                                                                                           <  2:893548/133‑1 (MQ=255)
                          tCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACa                                                                                                                          <  2:558785/150‑1 (MQ=255)
                          tCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACa                                                                                                                          <  2:266329/150‑1 (MQ=255)
                          tCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACa                                                                                                                          <  2:1084960/150‑1 (MQ=255)
                            aCGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGTCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGt                                                                                                                        <  2:1118168/150‑1 (MQ=255)
                                  gATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGcc                                                                                                                  >  1:872226/1‑150 (MQ=255)
                                        ccGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCAcc                                                                                                            >  2:434923/1‑150 (MQ=255)
                                          ggCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCgg                                                                                                          <  2:70302/150‑1 (MQ=255)
                                                gCTTAGAGGCGCTGATGGCAGCCGTACCCGGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGCCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACAt                                                                                                    <  2:198940/150‑1 (MQ=255)
                                                  tCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCg                                                                                                  >  2:1021018/1‑150 (MQ=255)
                                                  tCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGa                                                                                                 >  2:321331/1‑150 (MQ=255)
                                                    aGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGaa                                                                                                <  1:830291/150‑1 (MQ=255)
                                                      gggCGCTGATGGCAGCCGGACCCCGCTTGCCGACCAAAATGCCCGCTGCCAGGTTAGCCACGGCCACCGCATGGGGGAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAacc                                                                                              <  2:453108/150‑2 (MQ=255)
                                                          gcTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTc                                                                                          >  1:610561/1‑150 (MQ=255)
                                                                ggCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCCCCCGCGGGATTGGCCGCGGCAAGGGGGGAAATGACAGTATCGCCGGCACCGGTCGCATCGAACACTTCGCGTGc                                                                                    >  2:34950/1‑150 (MQ=255)
                                                                 gCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCaa                                                                                                                                        >  1:275038/1‑97 (MQ=255)
                                                                 gCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCaa                                                                                                                                        <  2:275038/97‑1 (MQ=255)
                                                                   aGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCcg                                                                                 >  2:203944/1‑150 (MQ=255)
                                                                   aGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCcg                                                                                 >  2:28225/1‑150 (MQ=255)
                                                                    gCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGGCACATCGAACACTTCGCGTGCCcgc                                                                                >  1:231847/1‑150 (MQ=255)
                                                                           ccAGCTTGCAGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGTGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCgcaggca                                                                         >  1:348727/1‑150 (MQ=255)
                                                                            cAGCTTGCCGAGCACTTGGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAAGGGCCGCGGCAGGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCgcaggcag                                                                        <  2:340322/150‑1 (MQ=255)
                                                                            cAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTcac                                                                                                      >  2:935377/1‑120 (MQ=255)
                                                                            cAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTcac                                                                                                      <  1:935377/120‑1 (MQ=255)
                                                                             aGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCCATTGCCGCgg                                                                                                                                           >  2:141719/1‑82 (MQ=255)
                                                                             aGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCCATTGCCGCgg                                                                                                                                           <  1:141719/82‑1 (MQ=255)
                                                                               cTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTg                                                                     <  1:610122/150‑1 (MQ=255)
                                                                                ttGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGc                                                                    >  1:17430/1‑150 (MQ=255)
                                                                                      aCCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCaa                                                                  <  1:887633/146‑1 (MQ=255)
                                                                                      aCCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCaa                                                                  >  2:887633/1‑146 (MQ=255)
                                                                                            aTGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGa                                                        <  1:1021018/150‑1 (MQ=255)
                                                                                                ccGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTACCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCaa                                                    <  1:433603/150‑1 (MQ=255)
                                                                                                          gTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAg                                          >  2:272518/1‑150 (MQ=255)
                                                                                                            tAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAggg                                        <  2:17430/150‑1 (MQ=255)
                                                                                                              gCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTc                                      >  1:522788/1‑150 (MQ=255)
                                                                                                               ccAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGa                                                                                                                            <  2:966241/63‑1 (MQ=255)
                                                                                                               ccAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGa                                                                                                                            >  1:966241/1‑63 (MQ=255)
                                                                                                               ccAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTcgc                                                                                        >  2:464905/1‑99 (MQ=255)
                                                                                                               ccAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTcgc                                                                                        <  1:464905/99‑1 (MQ=255)
                                                                                                                     ccACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTcatgccat                               >  2:537919/1‑150 (MQ=255)
                                                                                                                     ccACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTcatgccat                               <  2:876464/150‑1 (MQ=255)
                                                                                                                     ccACCGCAGGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAAGGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTcatgccat                               <  1:34950/150‑1 (MQ=255)
                                                                                                                       aCCGCATGGGGCAGGTCCTCACACGCGGCAATGGCCGCGACAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCACGTGCCCGCGCAGGCAGGTGCAACGCCGGCTAACCGACACACAATAAGGTaatgccatgc                             >  2:508202/1‑150 (MQ=255)
                                                                                                                        ccGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCt                            >  1:18808/1‑150 (MQ=255)
                                                                                                                            aTGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGccc                                                                                  <  2:599149/92‑1 (MQ=255)
                                                                                                                            aTGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGccc                                                                                  >  1:599149/1‑92 (MQ=255)
                                                                                                                                ggCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCgg                    <  1:537919/150‑1 (MQ=255)
                                                                                                                                   aGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTg                 <  1:272518/150‑1 (MQ=255)
                                                                                                                                   aGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTg                 >  1:600335/1‑150 (MQ=255)
                                                                                                                                    ggTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACCGTATCCCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTGa                <  2:870773/150‑1 (MQ=255)
                                                                                                                                          ttgccGCGGCCATTGCCGCGGCAAGGGTGGAAATGACAGTATCGACGGCAACGATCTCACTGAGCACTTCAGGGCACCGCGGGTTCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTGACCAGCa          >  1:689592/4‑150 (MQ=255)
                                                                                                                                                  gCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTGACCAGCAAGGCGccc  >  2:706491/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCAGCACGCCGCGTTCGGAGCCTTCCTCACGCTGGATGGCCCGGCGCAGCTCAGGGGCGCTGATGGCAGCCGTACCCAGCTTGCCGACCACAATGCCCGCTGCCAGGTTAGCCAGGGCCACCGCATGGGGCAGGTCCTCACCCGCGGCAATGGCCGCGGCAAGGGTGGAAATGACAGTATCGCCGGCACCGGTCACATCGAACACTTCGCGTGCCCGCGCAGGCAGGTGCAACGCCGGCTGACCAACACGCAACAGGGTCATGCCATGCTCACCGCGGGTGACCAGCAAGGCGCCC  >  NC_002947/5614606‑5614901

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: