Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,667,309 G→T 9.7% intergenic (+36/+35) nhaP → / ← PP_4975 NhaP‑type Na+(K+)/H+ antiporter/acyl‑CoA thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,667,3090GT 9.7% 332.4 / 13.8 109intergenic (+36/+35)nhaP/PP_4975NhaP‑type Na+(K+)/H+ antiporter/acyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base G (46/52);  new base T (5/6);  total (51/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

GGCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGTCGACAGTC  >  NC_002947/5667160‑5667458
                                                                                                                                                      |                                                                                                                                                     
gacGAGGAACGGGACCTGGGGCTGTCGATCACCTACATCGTGGTGGTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTg                                                                                                                                                       >  2:681013/3‑150 (MQ=255)
   gAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAAc                                                                                                                                                    >  2:2303650/1‑150 (MQ=255)
    aGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGAATCTGGTGCAGGGGTTGAGAATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAAcc                                                                                                                                                   >  1:1011486/1‑150 (MQ=255)
       aaCGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCcgc                                                                                                                                                <  2:2140773/150‑1 (MQ=255)
          gggACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGaa                                                                                                                                             <  2:1825721/150‑1 (MQ=255)
          gggACCTGCTGCTGGCGAGCACCTACATCGTGGTGTTGTCGGCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGaa                                                                                                                                             <  2:450462/150‑1 (MQ=255)
           ggACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAag                                                                                                                                            >  2:1941919/1‑150 (MQ=255)
            gACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGCTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGaatg                                                                                                                                           <  2:849215/150‑3 (MQ=255)
              ccTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGCTTCTGGTGCAGGGGTTGAGCACTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTGGAGAATGCCGGGGCCGCTTTGCGGCCGATTCGCGGGTGAACCCGCGAAGAgg                                                                                                                                         <  2:141211/150‑1 (MQ=255)
              ccTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGCGTTGCGAGTCCCAGGGTGGTGGCGCACGGTCAGCGCCCACCCTTCAGCATCGCGGACGAGCTTGGCGTCCCATTCGCGGGTTCACCCGCGaatggg                                                                                                                                         <  1:2280261/150‑5 (MQ=255)
               ctgctgctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGc                                                                                                                                        >  2:1028313/1‑150 (MQ=255)
                 gctgctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCg                                                                                                                                      <  2:1920335/150‑1 (MQ=255)
                   tgctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCAGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGggcggc                                                                                                                                    >  2:1570728/1‑150 (MQ=255)
                                cTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAgg                                                                                                                       <  1:265542/150‑1 (MQ=255)
                                cTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAgg                                                                                                                       >  1:549701/1‑150 (MQ=255)
                                  ccATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGGGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGt                                                                                                                     <  2:1879891/149‑1 (MQ=255)
                                     tCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTa                                                                                                                  <  1:1709427/150‑1 (MQ=255)
                                     tCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTa                                                                                                                  >  2:1857127/1‑150 (MQ=255)
                                       gtggtgTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTAtt                                                                                                                <  1:1605965/150‑1 (MQ=255)
                                              tgtCGTCGATCCTCGTGCAGGCTTTGAGCCTTGGGCGCGTGTTGAGCAAGGTCAGCGCCCCGCCCTGAGCATGCCGGGGGCGCTCTCCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGc                                                                                                         <  1:1064060/150‑1 (MQ=255)
                                               gtcgtcGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGcc                                                                                                        <  1:1488172/150‑1 (MQ=255)
                                                  gtcgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGaa                                                                                                     <  2:1571713/150‑1 (MQ=255)
                                                   tcgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAAt                                                                                                    >  1:618869/1‑150 (MQ=255)
                                                   tcgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAAt                                                                                                    >  1:1854401/1‑150 (MQ=255)
                                                    cgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATc                                                                                                   <  2:1854401/150‑1 (MQ=255)
                                                      aTTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCg                                                                                                 >  1:1126841/1‑150 (MQ=255)
                                                      aTTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCg                                                                                                 <  1:2303650/150‑1 (MQ=255)
                                                             tGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAgg                                                                                                                                         <  1:867079/103‑1 (MQ=255)
                                                             tGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAgg                                                                                                                                         >  2:867079/1‑103 (MQ=255)
                                                                  gggTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTAtt                                                                                                                <  2:1164994/123‑1 (MQ=255)
                                                                  gggTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTAtt                                                                                                                >  1:1164994/1‑123 (MQ=255)
                                                                          cATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTAc                                                                             >  1:1080631/1‑150 (MQ=255)
                                                                           aTTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACt                                                                            >  2:1170257/1‑150 (MQ=255)
                                                                            ttGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATCGGGAAATTAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACtt                                                                           <  1:1570728/150‑1 (MQ=255)
                                                                             tGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTg                                                                          <  1:1941919/150‑1 (MQ=255)
                                                                                gCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGAtc                                                                       >  2:2142374/1‑150 (MQ=255)
                                                                                   ggtggtGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTTCACGTACTTGATCTCg                                                                    >  2:479211/1‑150 (MQ=255)
                                                                                   ggtggtGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCg                                                                    <  1:2267692/150‑1 (MQ=255)
                                                                                   ggtggtGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCg                                                                    <  1:1997508/150‑1 (MQ=255)
                                                                                    gtggtgCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTgg                                                                                         >  1:288649/1‑128 (MQ=255)
                                                                                    gtggtgCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTgg                                                                                         <  2:288649/128‑1 (MQ=255)
                                                                                    gtggtgCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCgg                                                                   >  1:1476691/1‑150 (MQ=255)
                                                                                     tggtgCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGt                                                                  >  1:1431184/1‑150 (MQ=255)
                                                                                     tggtgCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCCGGTGAACCCGAGAAGAGTCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGt                                                                  <  1:1450155/150‑1 (MQ=255)
                                                                                       gtgCGCAAGGTCAGCGCCCAGCCTTTAGATTGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTgc                                                                >  1:315754/1‑150 (MQ=255)
                                                                                        tgCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGt                                                                  >  2:1472686/1‑147 (MQ=255)
                                                                                        tgCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGt                                                                  <  1:1472686/147‑1 (MQ=255)
                                                                                              aaGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGt                                                                                                                     >  2:347654/1‑91 (MQ=255)
                                                                                              aaGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTTCCGAACCCGGGAAGTTGTCCTGCACGTACTTGATCTCGGTGCGCCCAt                                                          <  1:1857127/150‑1 (MQ=255)
                                                                                                ggTCGGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATgg                                                        <  2:1445620/150‑1 (MQ=255)
                                                                                                 gTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCgg                                                                                                                                     <  2:1330522/72‑1 (MQ=255)
                                                                                                 gTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCgg                                                                                                                                     >  1:1330522/1‑72 (MQ=255)
                                                                                                 gTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATggg                                                       <  2:1126841/150‑1 (MQ=255)
                                                                                                 gTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATggg                                                       >  2:463799/1‑150 (MQ=255)
                                                                                                      cgcCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGcc                                                  <  2:618869/150‑1 (MQ=255)
                                                                                                       gcCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCg                                                 >  1:2099543/1‑150 (MQ=255)
                                                                                                        cccAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAAc                                                                                                                           <  2:1029715/75‑1 (MQ=255)
                                                                                                        cccAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAAc                                                                                                                           >  1:1029715/1‑75 (MQ=255)
                                                                                                         ccAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCAt                                                          <  1:1382458/139‑1 (MQ=255)
                                                                                                         ccAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCAt                                                          >  2:1382458/1‑139 (MQ=255)
                                                                                                           aGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCaa                                             <  2:2099543/150‑1 (MQ=255)
                                                                                                               ttGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCa                                         >  1:1199420/1‑150 (MQ=255)
                                                                                                                    aaTGCCGGGGCCGCTTTGCGGCCCATTCGCGTGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGt                                                                                        >  2:1996223/1‑98 (MQ=255)
                                                                                                                    aaTGCCGGGGCCGCTTTGCGGCCCATTCGCGTGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGt                                                                                        <  1:1996223/98‑1 (MQ=255)
                                                                                                                       gCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCgg                                                                                                >  1:318088/1‑87 (MQ=255)
                                                                                                                       gCCGGGGCCGCTTTGCGGCCCATTCGCGGGGGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCgg                                                                                                <  2:318088/87‑1 (MQ=255)
                                                                                                                        ccGGCGCCGCTTTGCGGCCCATTCGCGGGTGAACCGGCGAAGAGGCCGGCACACGCAACCAGGGTCTATTCGACTGCCGAATCCGAGAACTGGTCCAGCACGTACTTGATCCCGGTGCCCCCATGGGGTGCAGGCAAACCATCCTCGccc                                >  2:1292205/1‑150 (MQ=255)
                                                                                                                         cGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTAc                                                                             >  1:642307/1‑104 (MQ=255)
                                                                                                                         cGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTAc                                                                             <  2:642307/104‑1 (MQ=255)
                                                                                                                         cGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGt                                                                  <  1:1995558/115‑1 (MQ=255)
                                                                                                                         cGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGt                                                                  >  2:1995558/1‑115 (MQ=255)
                                                                                                                             gCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTc                                                                                       <  2:1896921/90‑1 (MQ=255)
                                                                                                                             gCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTc                                                                                       >  1:1896921/1‑90 (MQ=255)
                                                                                                                               cGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGa                         >  1:1255507/1‑150 (MQ=255)
                                                                                                                                gcaaagCGGCCCCTTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGAc                        <  1:1226586/145‑1 (MQ=255)
                                                                                                                                gCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGAc                        >  1:50230/1‑150 (MQ=255)
                                                                                                                                  tttGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGa                      >  1:790944/1‑150 (MQ=255)
                                                                                                                                       ggCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCt                                                                                     >  1:1326467/1‑82 (MQ=255)
                                                                                                                                       ggCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCt                                                                                     <  2:1326467/82‑1 (MQ=255)
                                                                                                                                         cccATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCAt               >  2:1869295/1‑150 (MQ=255)
                                                                                                                                         cccATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCAt               <  2:1476691/150‑1 (MQ=255)
                                                                                                                                          ctcttCGCGGGTTCACCCGCGAATAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCg                                                                    >  1:1600752/4‑96 (MQ=38)
                                                                                                                                          ctcttCGCGGGTTCACCCGCGAATAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCg                                                                    <  2:1600752/93‑1 (MQ=37)
                                                                                                                                          ccATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATc              >  1:29172/1‑150 (MQ=255)
                                                                                                                                           tcttCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCt                                                                                     <  2:1083443/76‑1 (MQ=37)
                                                                                                                                           tcttCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCt                                                                                     >  1:1083443/3‑78 (MQ=25)
                                                                                                                                           tcttCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCt             >  1:562948/3‑150 (MQ=255)
                                                                                                                                             ttCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCAt                                                          <  2:1900665/103‑1 (MQ=255)
                                                                                                                                             ttCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCAt                                                          >  1:1900665/1‑103 (MQ=255)
                                                                                                                                             ttCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTg           >  2:2319723/1‑150 (MQ=255)
                                                                                                                                             ttCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTg           <  1:463799/150‑1 (MQ=255)
                                                                                                                                               cgcgGGTGAACCCGCGAAGAGGCCGGCACAGGCAAAAAGGGTCTCTTTCACCGCCGAAACCGGCAACTGTTCCTGCACGTACTGTATCTCGGTGCTCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGTc         <  2:1431184/150‑1 (MQ=255)
                                                                                                                                                gcgGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCaa                                             >  2:1429530/1‑113 (MQ=255)
                                                                                                                                                gcgGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCaa                                             <  1:1429530/113‑1 (MQ=255)
                                                                                                                                                gcgGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATc                                       >  1:1162965/1‑119 (MQ=255)
                                                                                                                                                gcgGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATc                                       <  2:1162965/119‑1 (MQ=255)
                                                                                                                                                 ggggTTAACCCGCGAAGAGGCCGGCACAGTCAACCAGGTTCTATTCGACTGCCGAATCCGGCAATTTGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGGTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGTCGa       <  1:624181/149‑1 (MQ=255)
                                                                                                                                                   ggTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTg                                                                 <  2:1842121/90‑1 (MQ=255)
                                                                                                                                                   ggTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGTCGACa     <  1:1028313/150‑1 (MQ=255)
                                                                                                                                                      gAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGTCGACAGTc  <  1:412539/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GGCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGTCGACAGTC  >  NC_002947/5667160‑5667458

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