Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 462,592 A→T 8.3% intergenic (‑128/+13) pqqF ← / ← davA coenzyme PQQ synthesis protein F/5‑aminopentanamidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947462,5920AT 8.3% 231.4 / 8.3 80intergenic (‑128/+13)pqqF/davAcoenzyme PQQ synthesis protein F/5‑aminopentanamidase
Reads supporting (aligned to +/- strand):  ref base A (33/40);  new base T (4/3);  total (37/43)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TTGGCGGATGGCGTCAGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCAGCAAG  >  NC_002947/462444‑462737
                                                                                                                                                    |                                                                                                                                                    
ttGGCGGATGGCGTCAGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCa                                                                                                                                                      >  3:367600/1‑149 (MQ=255)
  ggCGGATGGCGTCAGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCa                                                                                                                                                      >  7:8650/1‑147 (MQ=255)
   gCGGATGGCGTCAGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAg                                                                                                                                                   >  7:350625/1‑149 (MQ=255)
   gCGGATGGCGTCAGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAg                                                                                                                                                   >  6:149090/1‑149 (MQ=255)
          gCGTCAGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTAGGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGccc                                                                                                                                             <  8:205634/148‑1 (MQ=255)
             tCAGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCTATGCGGcccattc                                                                                                                                         <  7:97556/149‑5 (MQ=255)
               aGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCCGCTATGCGGcccattcgc                                                                                                                                       <  5:377816/149‑7 (MQ=255)
                      ggggCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTa                                                                                                                                <  6:348625/149‑1 (MQ=255)
                             cGGTGTCTGCTCTCCATGCTTACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGGCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGGATAGCGGCCCCTGTCAGCCTTTACGCAg                                                                                                                           <  5:288631/147‑1 (MQ=255)
                             cGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAg                                                                                                                           <  8:8650/147‑1 (MQ=255)
                             cGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGt                                                                                                                         >  3:198236/1‑149 (MQ=255)
                             cGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGt                                                                                                                         >  2:112985/1‑149 (MQ=255)
                               gTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTg                                                                                                                        >  4:184975/1‑148 (MQ=255)
                                         tcCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTg                                                                                                               <  3:139506/147‑1 (MQ=255)
                                         tcCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTg                                                                                                               <  5:548/147‑1 (MQ=255)
                                         tcCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGgc                                                                                                             >  6:268531/1‑149 (MQ=255)
                                          ccATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGgcg                                                                                                            >  2:124039/1‑149 (MQ=255)
                                           cATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTgg                                                                                                              <  4:367600/146‑1 (MQ=255)
                                           cATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGgcg                                                                                                            >  6:390499/1‑148 (MQ=255)
                                              gCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTc                                                                                                                  >  5:188519/1‑139 (MQ=255)
                                              gCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTc                                                                                                                  <  6:188519/139‑1 (MQ=255)
                                                tAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATc                                                                                                      >  1:118724/1‑149 (MQ=255)
                                                 aaCCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCg                                                                                                    >  6:309540/1‑149 (MQ=255)
                                                 aaCCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCg                                                                                                    >  2:277891/1‑149 (MQ=255)
                                                 aaCCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCg                                                                                                    >  8:303559/1‑149 (MQ=255)
                                                  aCCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGc                                                                                                                       >  3:169235/1‑130 (MQ=255)
                                                  aCCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGc                                                                                                                       <  4:169235/130‑1 (MQ=255)
                                                  aCCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCg                                                                                                    <  7:154427/148‑1 (MQ=255)
                                                  aCCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCg                                                                                                    <  3:6947/148‑1 (MQ=255)
                                                   ccGATTTGTTTGCAGGCCCGGCCTTTTCGCGTGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGt                                                                                                  <  8:350625/149‑1 (MQ=255)
                                                    cGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGt                                                                                                  >  4:85966/1‑148 (MQ=255)
                                                     gATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTg                                                                                                                        <  3:234249/126‑1 (MQ=255)
                                                     gATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTg                                                                                                                        >  4:234249/1‑126 (MQ=255)
                                                          gtttgCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTa                                                                                                                                <  4:128358/113‑1 (MQ=255)
                                                          gtttgCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTa                                                                                                                                >  3:128358/1‑113 (MQ=255)
                                                             tgCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑ggg                                                                                       >  3:310473/1‑148 (MQ=255)
                                                              gCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGa                                                                                      >  7:35917/1‑148 (MQ=255)
                                                              gCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGa                                                                                      >  3:279270/1‑148 (MQ=255)
                                                                   cccGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGcc                                                                                >  2:280700/1‑149 (MQ=255)
                                                                   cccGCCCTTTTCGCCTGTATGCCGCCGCTATCAAAGCCACCGCAAATCCCGTGGGAGCGCGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGcc                                                                                <  4:310473/149‑1 (MQ=255)
                                                                    ccGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAATTCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGccc                                                                               <  1:176136/149‑1 (MQ=255)
                                                                    ccGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGccc                                                                               >  7:363242/1‑149 (MQ=255)
                                                                     cGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCt                                                                              <  1:124039/149‑1 (MQ=255)
                                                                      ggCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTg                                                                             <  4:134704/149‑1 (MQ=255)
                                                                         cTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTgcgc                                                                          <  1:343129/149‑1 (MQ=255)
                                                                          ttttCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGcc                                                                         >  3:3151/1‑149 (MQ=255)
                                                                             tCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGcccc                                                                       <  3:85966/148‑1 (MQ=255)
                                                                              cgcgGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCt                                                                      <  5:149090/148‑1 (MQ=255)
                                                                               gcgGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTg                                                                     <  8:35917/148‑1 (MQ=255)
                                                                                 gggTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGcac                                                                  <  4:279270/149‑1 (MQ=255)
                                                                                 gggTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGcac                                                                  >  5:389407/1‑149 (MQ=255)
                                                                                    tATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGcaccac                                                               >  6:117525/1‑149 (MQ=255)
                                                                                      tGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAACCGCCCTGCGCCCCTGCACCAccc                                                             >  6:266405/1‑149 (MQ=255)
                                                                                        gcgcCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCg                                                            >  8:233696/1‑148 (MQ=255)
                                                                                          gcCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCt                                                          <  6:303753/148‑1 (MQ=255)
                                                                                          gcCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCt                                                          >  5:303753/1‑148 (MQ=255)
                                                                                          gcCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTc                                                         <  5:268531/149‑1 (MQ=255)
                                                                                               tATGAAAGGCAGCGCAAATCCCGTGTGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCt                                                    <  4:198236/149‑1 (MQ=255)
                                                                                               tATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTAc                                                                                                                               <  3:42101/77‑1 (MQ=255)
                                                                                               tATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTAc                                                                                                                               >  4:42101/1‑77 (MQ=255)
                                                                                                 tGAAAGGCAGCGCAAATCCCGTGGGAGCGGGTGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGaa                                                                                     <  3:304442/114‑1 (MQ=255)
                                                                                                 tGAAAGGCAGCGCAAATCCCGTGGGAGCGGGTGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGaa                                                                                     >  4:304442/1‑114 (MQ=255)
                                                                                                 tGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTc                                                   <  2:81602/148‑1 (MQ=255)
                                                                                                   taaGGCATCGCAAATCCCGTGGGATCGGGCGTGCCCGCGTATTGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCATGTGCATCTCCTTGCGGAAATCTGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCg                                                  <  1:280700/147‑1 (MQ=255)
                                                                                                   aaaGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCCCCCGCTCATGCTCGAg                                                >  6:206397/1‑149 (MQ=255)
                                                                                                        cAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTg                                           >  1:114653/1‑149 (MQ=255)
                                                                                                         aGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGc                                          <  1:277891/149‑1 (MQ=255)
                                                                                                            gcAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCcaa                                       <  8:363242/149‑1 (MQ=255)
                                                                                                             cAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCcaac                                      <  5:117525/149‑1 (MQ=255)
                                                                                                                 tCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACt                                  <  3:184975/149‑1 (MQ=255)
                                                                                                                   ccGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTgg                                <  2:118724/149‑1 (MQ=255)
                                                                                                                        ggAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTc                           <  6:389407/149‑1 (MQ=255)
                                                                                                                           gCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCATGTATTCGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCa                          <  1:112985/149‑1 (MQ=255)
                                                                                                                           gCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATc                        <  5:390499/149‑1 (MQ=255)
                                                                                                                             gggCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATccc                      <  4:3151/149‑1 (MQ=255)
                                                                                                                                 gTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGctct                                            <  2:371806/123‑1 (MQ=255)
                                                                                                                                       gcgAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCAt            >  1:261100/1‑149 (MQ=255)
                                                                                                                                          caggggCCGCTATGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGc         >  7:296911/4‑149 (MQ=255)
                                                                                                                                             gggCCGCTATGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGTCCGGCCATg           >  7:183566/1‑144 (MQ=255)
                                                                                                                                             gggCCGCTATGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGTCCGGCCATg           <  8:183566/144‑1 (MQ=255)
                                                                                                                                             gggCCGCTATGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCAg      >  5:307564/1‑149 (MQ=255)
                                                                                                                                             gggCCGCTATGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCAg      >  5:302364/1‑149 (MQ=255)
                                                                                                                                                    aTAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCAGCAAg  >  2:368488/1‑146 (MQ=255)
                                                                                                                                                    aTAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCAGCAAg  <  1:368488/146‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTGGCGGATGGCGTCAGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCAGCAAG  >  NC_002947/462444‑462737

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: