Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 460,473 T→C 5.7% A665A (GCA→GCG pqqF ← coenzyme PQQ synthesis protein F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947460,4730TC 5.7% 207.0 / 2.4 70A665A (GCA→GCGpqqFcoenzyme PQQ synthesis protein F
Reads supporting (aligned to +/- strand):  ref base T (30/36);  new base C (2/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

TGTGGCCAAGTGTGTCTGCCACGCCCATTCGGCGACATCGGCATTGGCCATGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTACCCGCAGGCGCTGGTAAACCGGCCCTTGCAGCAATTGGGCTAACAGGCG  >  NC_002947/460326‑460622
                                                                                                                                                   |                                                                                                                                                     
tgtgGCCAAGTGTGTCTGGCACGCCCATTCGGCGACATCGGCATTGGCCATGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCTAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGc                                                                                                                                                     <  2:157660/150‑1 (MQ=255)
tgtgGCCAAGTGTGTCTGCCACGCCCATTCGGCGACATCGGCATTGGCCATGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGc                                                                                                                                                     <  2:863403/150‑1 (MQ=255)
     ccAAGTGTGTCTGCCACGCCCATTCGGCGACATCGGCATTGGCCATGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCGGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGa                                                                                                                                                 <  2:1245026/149‑1 (MQ=255)
     ccAAGTGTGTCTGCCACGCCCATTCGGCGACATCGGCATTGGCCATGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCGCCTGTGCCTGAt                                                                                                                                                >  1:1298915/1‑150 (MQ=255)
     ccAAGTGTGTCTGCCACGCCCATTCGGCGACATCGGCATTGGCCATGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGa                                                                                                                                                 >  1:1245026/1‑149 (MQ=255)
      cAAGTGTGTCTGCCACGCCCATTCGGCGACATCGGCATTGGCCATGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATc                                                                                                                                               >  2:1023135/1‑150 (MQ=255)
             gtgtGCCACGCCCATTCGGCGACATCGGCATTGGCCATGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGCGCCTGATCGGTGTgg                                                                                                                                        >  1:1935594/4‑150 (MQ=255)
                tGCCACGCCCATTCGGCGACATCGGCATTGGCCATGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTggggg                                                                                                                                     <  2:1285684/150‑1 (MQ=255)
                 gCCACGCCCATTCGGCGACATCGGCATTGGCCATGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTgggggg                                                                                                                                    >  2:1297808/1‑150 (MQ=255)
                     cGCTCATTCGGCGACATCGGCATTGGCCATGGCGGGCTGGTCGAATTGCCCAGCCAATGTCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGGAGGTTGAGCAGGTGGTCGAGCAACTGTGCCTGATCGGTGTGGGGGGATTg                                                                                                                                <  2:802044/150‑1 (MQ=255)
                      gCCCATTCGGCGACATCGGCATTGGCCATGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGc                                                                                                                                                     >  1:318027/1‑128 (MQ=255)
                      gCCCATTCGGCGACATCGGCATTGGCCATGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGc                                                                                                                                                     <  2:318027/128‑1 (MQ=255)
                                   cATCGGCATTGGCCATGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTGAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAAcagca                                                                                                                  <  1:1484234/150‑1 (MQ=255)
                                   cATCGGCAGTGGCCATGGTGGGCTCGTCGAATTGCCCAGCCAATGCCGGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGATCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAAcagca                                                                                                                  <  1:15831/150‑1 (MQ=255)
                                       ggCATTGGCCATGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGccc                                                                                                              >  1:1248233/1‑150 (MQ=255)
                                        gCATTGGCCATGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCa                                                                                                             >  1:3378/1‑150 (MQ=255)
                                         cATAGGCCATGGCGGGCTCGTCGAATGGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCAc                                                                                                            <  1:478213/150‑1 (MQ=255)
                                                 aTGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGCGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGcccccc                                                                                                    <  2:1935594/150‑1 (MQ=255)
                                                          tcgtcgAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGa                                                                                           <  1:1122327/150‑1 (MQ=255)
                                                             tcgAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGa                                                                                        >  2:1466857/1‑150 (MQ=255)
                                                                    gCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGgcgc                                                                                 <  2:1258615/150‑1 (MQ=255)
                                                                     cccAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCt                                                                                >  1:496319/1‑150 (MQ=255)
                                                                     cccAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCt                                                                                >  1:106647/1‑150 (MQ=255)
                                                                       cAGCCAATGCCTGGCGAGCGGCCGGGTCACGGGGGACGCCATGGGGCAGCAGGCTGAGCACGTGGTCGAGGACCTGTGCCTGAACGGTGTGGGGGGGTTGCACGCCGAACCGCAGCCCACCGACGCCCTCGACCTGGAGAAAGGCGCTgc                                                                              >  2:456531/1‑149 (MQ=255)
                                                                       cAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGc                                                                                                                           <  2:553912/105‑1 (MQ=255)
                                                                       cAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGc                                                                                                                           >  1:553912/1‑105 (MQ=255)
                                                                        aGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGaa                                                                             >  1:1117577/1‑150 (MQ=255)
                                                                           cAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAacac                                                                          >  2:625926/1‑150 (MQ=255)
                                                                                  tGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAg                                                                   <  1:1023135/150‑1 (MQ=255)
                                                                                  tGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAg                                                                   <  1:1008987/150‑1 (MQ=255)
                                                                                    gCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCGGACGAAACGCGATGAACACCGGCTAGca                                                                 >  2:1288525/1‑149 (MQ=255)
                                                                                      gAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAGCAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGa                                                               <  2:1255571/150‑1 (MQ=255)
                                                                                      gAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGCGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCt                                                                                             >  1:1362863/1‑120 (MQ=255)
                                                                                      gAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGCGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCt                                                                                             <  2:1362863/120‑1 (MQ=255)
                                                                                       agaggacGGGTCTAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGTGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAg                                                              <  2:515183/144‑1 (MQ=255)
                                                                                        gCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGaa                                                                                                                       >  2:915561/1‑92 (MQ=255)
                                                                                        gCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGaa                                                                                                                       >  2:916159/1‑92 (MQ=255)
                                                                                        gCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGaa                                                                                                                       <  1:916159/92‑1 (MQ=255)
                                                                                        gCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGaa                                                                                                                       <  1:915561/92‑1 (MQ=255)
                                                                                        gCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCACGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGt                                                             <  1:8821/150‑1 (MQ=255)
                                                                                         cGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGtt                                                            <  2:1117577/150‑1 (MQ=255)
                                                                                           gCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGc                                                          >  2:496573/1‑150 (MQ=255)
                                                                                             cGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCaa                                                        <  1:1420343/150‑1 (MQ=255)
                                                                                               ggTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGt                                                      <  1:625926/150‑1 (MQ=255)
                                                                                               ggTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGt                                                      <  1:1297808/150‑1 (MQ=255)
                                                                                                   aaGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGAtt                                                                                                                                 >  1:1266107/1‑71 (MQ=255)
                                                                                                   aaGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGAtt                                                                                                                                 <  2:1266107/71‑1 (MQ=255)
                                                                                                   aaGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGAc                                                                                               >  1:1229578/1‑105 (MQ=255)
                                                                                                   aaGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGAc                                                                                               <  2:1229578/105‑1 (MQ=255)
                                                                                                    aGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTAcc                                                 <  2:1248233/150‑1 (MQ=255)
                                                                                                         gACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTACCCGCAg                                            <  1:496573/150‑1 (MQ=255)
                                                                                                              cATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGa                                                               <  2:282618/126‑1 (MQ=255)
                                                                                                              cATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGa                                                               >  1:282618/1‑126 (MQ=255)
                                                                                                               aTGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTACCCGCAGGCGCTg                                      <  1:1497495/150‑1 (MQ=255)
                                                                                                               aTGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTACCCGCAGGCGCTg                                      >  2:783893/1‑150 (MQ=255)
                                                                                                                          ggTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGt                                                             <  1:1732995/116‑1 (MQ=255)
                                                                                                                          ggTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGt                                                             >  2:1732995/1‑116 (MQ=255)
                                                                                                                                    gtggtCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGCCCCCCTCGCCCTGACGAACGGCGCTGAACACCGCTGAGCCGATTCTCAATTCTACCCGCAGGCGCTGGTAAACCGGCCCTTGCAGCaa                 >  2:218197/1‑150 (MQ=255)
                                                                                                                                    gtggtCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACAcc                                                                         >  2:1923336/1‑94 (MQ=255)
                                                                                                                                    gtggtCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACAcc                                                                         <  1:1923336/94‑1 (MQ=255)
                                                                                                                                    gtggtCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTACCCGCAGGCGCTGGTAAACCGGCCCTTGCAGCaa                 >  1:317703/1‑150 (MQ=255)
                                                                                                                                       gtCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTACCCGCAgg                                           >  2:979203/1‑121 (MQ=255)
                                                                                                                                       gtCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTACCCGCAgg                                           <  1:979203/121‑1 (MQ=255)
                                                                                                                                       gtCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTACCCGCAGGCGCTGGTAAACCgg                             <  2:655641/135‑1 (MQ=255)
                                                                                                                                       gtCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTACCCGCAGGCGCTGGTAAACCgg                             >  1:655641/1‑135 (MQ=255)
                                                                                                                                       gtCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTACCCGCAGGCGCTGGTAAACCGGCCCTTGCAGCAATTg              <  2:3378/150‑1 (MQ=255)
                                                                                                                                        tCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTACCCGCAGGCGCTGGTAAACCGGCCCTTGCAGCAATTgg             >  2:1535804/1‑150 (MQ=255)
                                                                                                                                          gAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTACCCGCAGGCGCTGGTAAACCGGCCCTTGCAGCAATTGGGc           >  1:1783487/1‑150 (MQ=255)
                                                                                                                                             cACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTACCCGCAGGCGCTGGTAAACCGGCCCTTGCAGCAATTGGGCTaa        <  2:496319/150‑1 (MQ=255)
                                                                                                                                                cTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTACCCGCAGGCGCTGGTAAACCGGCCCTTGCAGCAATTGGGCTAACAg     <  1:1466857/150‑1 (MQ=255)
                                                                                                                                                   tgCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACGCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTACCCGCAGGCGCTGGTAAACCGGCCCTTGCAGCAATTGGGCTAACAggcg  <  2:1298915/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGTGGCCAAGTGTGTCTGCCACGCCCATTCGGCGACATCGGCATTGGCCATGGCGGGCTCGTCGAATTGCCCAGCCAATGCCTGGCGAGCGGCCGGGTCAAGGGTGACGCCATGGCGCAGCAGGTTGAGCAGGTGGTCGAGCACCTGTGCCTGATCGGTGTGGGGGGATTGCACGCCGAACAGCAGCCCACCGACCCCCTCGACCTGACGAAAGGCGCTGAACACCGCGTAGCCGAGTTGCAAGTCTACCCGCAGGCGCTGGTAAACCGGCCCTTGCAGCAATTGGGCTAACAGGCG  >  NC_002947/460326‑460622

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