Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,904,998 G→T 8.4% 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 8.4% 172.4 / 7.1 60intergenic (‑165/+70)PP_4312/PP_4313hypothetical protein/peptidylprolyl isomerase
Reads supporting (aligned to +/- strand):  ref base G (31/24);  new base T (3/2);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.01e-01

CAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGCCCCTACGC  >  NC_002947/4904852‑4905144
                                                                                                                                                  |                                                                                                                                                  
cAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTGGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGc                                                                                                                                                  <  5:265458/149‑1 (MQ=255)
 agagAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGAACGGCGCCCGTGCAGTGAATGACATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTACTCCCTGGAGCTGAATGCCGTTCCCTAGGCAAGGcc                                                                                                                                                 >  2:366038/1‑149 (MQ=255)
   agaATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGACATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGTCCTa                                                                                                                                               <  1:366038/149‑1 (MQ=255)
      atatCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGt                                                                                                                                            <  5:199653/149‑1 (MQ=255)
         tCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAg                                                                                                                                         >  8:268362/1‑149 (MQ=255)
         tCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAg                                                                                                                                         <  6:221483/149‑1 (MQ=255)
           gggTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCCTTGCCTAgggaacg                                                                                                                                        <  2:230102/148‑6 (MQ=255)
                 aCCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGCCAAGGCCTAGGTTAGGGGCTgg                                                                                                                                  <  4:172541/148‑1 (MQ=255)
                   ccGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGgtg                                                                                                                               >  5:266089/1‑149 (MQ=255)
                   ccGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGgtg                                                                                                                               >  4:11601/1‑149 (MQ=255)
                    cGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCCTTGCCTAgggaa                                                                                                                                          >  3:12297/1‑134 (MQ=255)
                    cGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCCTTGCCTAgggaa                                                                                                                                          <  4:12297/137‑4 (MQ=255)
                    cGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCAGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGgtgt                                                                                                                              >  8:177059/1‑149 (MQ=255)
                          aTCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGGGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTgcgc                                                                                                                         >  1:296493/1‑148 (MQ=255)
                          aTCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGGTAGGGGCTGGGTGTTgcgc                                                                                                                         >  3:51086/1‑148 (MQ=255)
                                   cTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTaa                                                                                                               <  7:268362/149‑1 (MQ=255)
                                      ctagccgggCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGt                                                                                                            <  8:166349/143‑1 (MQ=255)
                                          cAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCaa                                                                                                        <  2:225044/149‑1 (MQ=255)
                                           aGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCTAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAg                                                                                                       <  1:145559/149‑1 (MQ=255)
                                                       cAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTg                                                                                           >  5:185869/1‑149 (MQ=255)
                                                          tGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCa                                                                                         >  6:290661/1‑148 (MQ=255)
                                                           gAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAAc                                                                                       >  7:311674/1‑149 (MQ=255)
                                                             aTGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGt                                                                                     <  4:330585/149‑1 (MQ=255)
                                                                 cATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGACCCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAac                                                                                 >  7:341435/1‑149 (MQ=255)
                                                                 cATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAAtt                                                                                                                     <  2:122662/113‑1 (MQ=255)
                                                                 cATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTGAGGGGCTGGGTGTTGCGCAAtt                                                                                                                     >  1:122662/1‑113 (MQ=255)
                                                                      ttGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTc                                                                             >  1:330592/1‑148 (MQ=255)
                                                                           aaCGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGtt                                                                                                                             <  3:146851/95‑1 (MQ=255)
                                                                           aaCGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGtt                                                                                                                             >  4:146851/1‑95 (MQ=255)
                                                                           aaCGCGTGAGCGTCACAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGAAAGGCCTAGCTTAGGGGCTGGGTCATCCGCAATTGTGAAAGGTGCAAGAACCTCGGGTTCCAACGAGAACACTCAatcat                                                                       >  2:47598/1‑149 (MQ=255)
                                                                                   aGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTg                                                               <  5:290661/149‑1 (MQ=255)
                                                                                    gggTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGaa                                                                                  <  7:101574/129‑1 (MQ=255)
                                                                                    gggTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGaa                                                                                  >  8:101574/1‑129 (MQ=255)
                                                                                            gTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACCTGAACACTCAATCATCCCGCGTGAGCACTTcc                                                      <  8:311674/149‑1 (MQ=255)
                                                                                               ccccGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCcagc                                                   <  8:341435/149‑1 (MQ=255)
                                                                                                        cGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGAtc                                          <  2:296493/149‑1 (MQ=255)
                                                                                                        cGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGAtc                                          <  2:330592/149‑1 (MQ=255)
                                                                                                        cGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGAtc                                          >  3:306715/1‑149 (MQ=255)
                                                                                                               ttCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAagg                                   >  2:141426/1‑149 (MQ=255)
                                                                                                                 cctccCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAaggt                                  >  5:25562/1‑148 (MQ=255)
                                                                                                                  ctccCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAaggtc                                 <  4:306715/148‑1 (MQ=255)
                                                                                                                  ctccCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAaggtc                                 <  6:185869/148‑1 (MQ=255)
                                                                                                                  ctccCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAaggtc                                 <  3:11601/148‑1 (MQ=255)
                                                                                                                    cccTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAaggtcagg                              >  4:3816/1‑149 (MQ=255)
                                                                                                                               aaTGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTGAGGTCCGAATTcgg                   >  5:119598/1‑149 (MQ=255)
                                                                                                                               aaTGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTcgg                   >  7:186143/1‑149 (MQ=255)
                                                                                                                                  gCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCATTTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTcggcgg                >  4:7313/1‑149 (MQ=255)
                                                                                                                                    cGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGa              <  1:141426/149‑1 (MQ=255)
                                                                                                                                      taaccTAGGCCTTGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATc            >  5:357492/4‑149 (MQ=255)
                                                                                                                                      gtCCCTATGCAAGGCCTAGGTTCGTGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATc            <  6:266089/148‑1 (MQ=255)
                                                                                                                                       aaccTAGGCCTTGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAAACTCGGGTTGCAACGTGAACACTCAATCATCACGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCg           >  7:111358/3‑149 (MQ=255)
                                                                                                                                        cccTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGc          >  3:135492/1‑149 (MQ=255)
                                                                                                                                         ccTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGcc         >  5:372692/1‑149 (MQ=255)
                                                                                                                                         ccTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGc          <  7:177059/148‑1 (MQ=255)
                                                                                                                                          cTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGccc        >  3:162734/1‑149 (MQ=255)
                                                                                                                                           tAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGcccc       >  2:188916/1‑149 (MQ=255)
                                                                                                                                           tAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGcccc       >  4:378100/1‑149 (MQ=255)
                                                                                                                                             ggCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGCCCCTa     >  3:65585/1‑149 (MQ=255)
                                                                                                                                              gCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGCCCCTAc    >  1:93571/1‑149 (MQ=255)
                                                                                                                                                aaGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGCCCCTACGc  >  1:172334/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTGAGGGTCAGAGTGCCCCGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAAGGTGCAAGAACCTCGGGTTGCAACGTGAACACTCAATCATCCCGCGTGAGCACTTCCAGCAGTTCGATCTCGAAGGTCAGGTCCGAATTCGGCGGGATCGCCCCTACGC  >  NC_002947/4904852‑4905144

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