Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,904,998 G→T 5.2% intergenic (‑165/+70) PP_4312 ← / ← PP_4313 hypothetical protein/peptidylprolyl isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,904,9980GT 5.2% 279.4 / 3.6 87intergenic (‑165/+70)PP_4312/PP_4313hypothetical protein/peptidylprolyl isomerase
Reads supporting (aligned to +/- strand):  ref base G (45/37);  new base T (3/2);  total (48/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

GCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGCCCCTACGC  >  NC_002947/4904851‑4905144
                                                                                                                                                   |                                                                                                                                                  
gCAGAGAATATCGGGTTTACCTGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAgg                                                                                                                                                   <  2:179828/149‑1 (MQ=255)
     gaATATCGGGTTTACCCGCGACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAg                                                                                                                                              <  2:186638/149‑1 (MQ=255)
     gaATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAg                                                                                                                                              <  8:79739/149‑1 (MQ=255)
         atCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGt                                                                                                                                            >  5:56341/1‑147 (MQ=255)
         atCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCCTTTGGTTCCTCCCTGGAGGTGAAAGCCGTTCCCTAGGCAAGGCCTAGGt                                                                                                                                            >  7:293841/1‑147 (MQ=255)
          tCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGTGTGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTTAATGCCGTTCCCTAGGCAAGGCCTAGGTTAg                                                                                                                                         <  1:20611/149‑1 (MQ=255)
          tCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAg                                                                                                                                         <  7:330707/149‑1 (MQ=255)
          tCGGGTTTACCCGCCACATCGCCCGGCAAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCCTTGCCTAgggaac                                                                                                                                         <  7:365300/149‑5 (MQ=255)
             ggTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAgggg                                                                                                                                      >  6:93741/1‑149 (MQ=255)
               tttACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCGTAGGTTAGGGGCt                                                                                                                                    >  7:5159/1‑149 (MQ=255)
                   cccGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAgg                                                                                                                                        <  7:380473/141‑1 (MQ=255)
                   cccGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAgg                                                                                                                                        >  8:380473/1‑141 (MQ=255)
                      gCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGtt                                                                                                                             >  4:300399/1‑149 (MQ=255)
                      gCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGtt                                                                                                                             <  6:57494/149‑1 (MQ=255)
                       ccACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTg                                                                                                                            <  3:8948/149‑1 (MQ=255)
                          caTCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTgcg                                                                                                                          >  4:134669/1‑148 (MQ=255)
                           aTCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTgcgc                                                                                                                         >  1:223864/1‑148 (MQ=255)
                            tCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCaa                                                                                                                       >  5:351447/1‑149 (MQ=255)
                            tCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCaa                                                                                                                       >  3:45845/1‑149 (MQ=255)
                              gCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTGCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAAtt                                                                                                                     <  1:418876/149‑1 (MQ=255)
                                ccGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATtgt                                                                                                                   >  1:26891/1‑149 (MQ=255)
                                ccGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATtgt                                                                                                                   >  4:206991/1‑149 (MQ=255)
                                          gCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGc                                                                                                          <  6:351447/148‑1 (MQ=255)
                                           cAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGc                                                                                                          >  7:77878/1‑147 (MQ=255)
                                            aGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAg                                                                                                       <  4:11143/149‑1 (MQ=255)
                                                   ccGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCg                                                                                                <  3:206991/149‑1 (MQ=255)
                                                        cAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGtt                                                                                            <  4:45845/148‑1 (MQ=255)
                                                        cAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTg                                                                                           >  5:166734/1‑149 (MQ=255)
                                                              aTGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGt                                                                                     <  5:414719/149‑1 (MQ=255)
                                                              aTGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGt                                                                                     <  3:235351/149‑1 (MQ=255)
                                                                gCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTg                                                                                    >  7:150300/1‑148 (MQ=255)
                                                                  cATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGaa                                                                                  >  5:77958/1‑148 (MQ=255)
                                                                       ttGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAGCCTCGGGTTGCAACGTGAACACTc                                                                             >  1:121/1‑148 (MQ=255)
                                                                       ttGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTc                                                                             >  6:369029/1‑148 (MQ=255)
                                                                            aaCGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAatcat                                                                       >  6:136231/1‑149 (MQ=255)
                                                                             aCGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAatcatc                                                                      <  5:369029/149‑1 (MQ=255)
                                                                                 gTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCcgc                                                                  >  8:82969/1‑149 (MQ=255)
                                                                                 gTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCcgc                                                                  <  1:226722/149‑1 (MQ=255)
                                                                                     gggTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGa                                                              <  5:292051/149‑1 (MQ=255)
                                                                                      ggTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAATACTCAATCATCCCGCGTGAg                                                             >  6:18204/1‑149 (MQ=255)
                                                                                      ggTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGa                                                              <  3:300399/148‑1 (MQ=255)
                                                                                              tGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCca                                                     <  6:166734/149‑1 (MQ=255)
                                                                                                 cccGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTc                                                       >  1:230526/1‑144 (MQ=255)
                                                                                                 cccGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTc                                                       <  2:230526/144‑1 (MQ=255)
                                                                                                  ccGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAatcatc                                                                      <  8:146544/128‑1 (MQ=255)
                                                                                                  ccGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAatcatc                                                                      >  7:146544/1‑128 (MQ=255)
                                                                                                  ccGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGa                                                              <  2:317320/136‑1 (MQ=255)
                                                                                                  ccGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGa                                                              >  1:317320/1‑136 (MQ=255)
                                                                                                  ccGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCcagcag                                                 >  6:249413/1‑149 (MQ=255)
                                                                                                  ccGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCca                                                     <  1:176564/145‑1 (MQ=255)
                                                                                                  ccGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCca                                                     >  2:176564/1‑145 (MQ=255)
                                                                                                        ccGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGAt                                           >  3:429738/1‑149 (MQ=255)
                                                                                                         cGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAGCCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGAtc                                          <  2:121/149‑1 (MQ=255)
                                                                                                         agTTTGGTTCCTGCCTGGAGGTTAATGCGGTTCCCTTGGCAATGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTTAGGTGCAAGAACCTCTGGTTGCAACGTGAACACTCAATCATCAGGAGTGAGCTCTTCCAGCAGTTctttc                                          <  8:5159/148‑5 (MQ=255)
                                                                                                          gTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAATACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGAtct                                         <  5:18204/149‑1 (MQ=255)
                                                                                                          gTTTGGTTCCTCCCTGGAGGTGAATGCAGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGAAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGAtct                                         >  2:397518/1‑149 (MQ=255)
                                                                                                               gTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAag                                    <  2:26891/149‑1 (MQ=255)
                                                                                                               gTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAag                                    <  7:82969/149‑1 (MQ=255)
                                                                                                                  cctccCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAaggt                                  >  8:314253/1‑148 (MQ=255)
                                                                                                                  cctccCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAaggt                                  >  6:82210/1‑148 (MQ=255)
                                                                                                                   ctccCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCCCGAaggtc                                 <  8:150300/148‑1 (MQ=255)
                                                                                                                      ccTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAaggtcag                               <  2:223864/147‑1 (MQ=255)
                                                                                                                      ccTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTAGATCTCGAaggtcaggt                             >  6:211134/1‑149 (MQ=255)
                                                                                                                        tGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCcag                                                    >  3:99376/1‑124 (MQ=255)
                                                                                                                        tGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCcag                                                    <  4:99376/124‑1 (MQ=255)
                                                                                                                          gAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCg                          >  1:180807/1‑148 (MQ=255)
                                                                                                                          gAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGa                         <  4:429738/149‑1 (MQ=255)
                                                                                                                           aGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGaa                        <  5:136231/149‑1 (MQ=255)
                                                                                                                                 aTGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTcggc                  <  7:314253/149‑1 (MQ=255)
                                                                                                                                  tGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGAtctc                                        >  8:94548/1‑126 (MQ=255)
                                                                                                                                  tGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGAtctc                                        <  7:94548/126‑1 (MQ=255)
                                                                                                                                  tGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTcggcg                 >  3:165506/1‑149 (MQ=255)
                                                                                                                                  tGCCGTTCCCTAGGACAGGCCTAGGTTCGGGGCTGGGTGGTGCGCAATTGTGTAAGGTGCTAGAACCTCGGGTTGCTACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTcggcg                 <  6:275162/149‑1 (MQ=255)
                                                                                                                                     cGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGa              >  5:148528/1‑149 (MQ=255)
                                                                                                                                     cGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGa              >  3:75612/1‑149 (MQ=255)
                                                                                                                                      ctaaccTAGGCCTTGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGAt             >  2:14369/5‑149 (MQ=255)
                                                                                                                                      ctaaccTAGGCCTTGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGAt             >  2:433897/5‑149 (MQ=255)
                                                                                                                                        aaccTAGGCCTTGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCGAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCg           >  5:290936/3‑149 (MQ=255)
                                                                                                                                          ccTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGcc         >  7:150906/1‑149 (MQ=255)
                                                                                                                                          ccTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGc          <  6:56341/148‑1 (MQ=255)
                                                                                                                                          ccTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGc          <  2:180807/148‑1 (MQ=255)
                                                                                                                                           cTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGccc        >  4:242715/1‑149 (MQ=255)
                                                                                                                                           cTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGccc        >  4:315029/1‑149 (MQ=255)
                                                                                                                                           cTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGccc        >  7:129314/1‑149 (MQ=255)
                                                                                                                                               gCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCAGTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGTGATCGCGCCTAc    >  8:330552/1‑149 (MQ=255)
                                                                                                                                               gCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGCCCCTAc    >  8:336473/1‑149 (MQ=255)
                                                                                                                                               gCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGCCCCTAc    >  8:73480/1‑149 (MQ=255)
                                                                                                                                                 aaGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGATCACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGCCCCTACGc  >  5:425012/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGCCCCTACGC  >  NC_002947/4904851‑4905144

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