Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 460,938 G→C 6.0% I509M (ATC→ATG pqqF ← coenzyme PQQ synthesis protein F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947460,9380GC 6.0% 263.1 / 6.0 83I509M (ATC→ATGpqqFcoenzyme PQQ synthesis protein F
Reads supporting (aligned to +/- strand):  ref base G (39/39);  new base C (3/2);  total (42/41)
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.97e-01

CTCGAAGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACT  >  NC_002947/460794‑461086
                                                                                                                                                |                                                                                                                                                    
cTCGAAGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATc                                                                                                                                                   <  1:408787/148‑1 (MQ=255)
cTCGAAGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCa                                                                                                                                                  <  1:79437/149‑1 (MQ=255)
  cGAAGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAgg                                                                                                                                                >  3:145806/1‑149 (MQ=255)
   gAAGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAgg                                                                                                                                                <  8:245234/148‑1 (MQ=255)
       ccTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCgg                                                                                                                                           <  7:331876/149‑1 (MQ=255)
        cTGCGGCCATGCCATGCCACGCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGAGCAGGGCCggg                                                                                                                                          <  7:311146/149‑1 (MQ=255)
         tGCGGCCATGCCGTGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCgggg                                                                                                                                         >  8:291235/1‑149 (MQ=255)
          gCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCggggg                                                                                                                                        >  1:80623/1‑149 (MQ=255)
              ccatgccatgCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAgc                                                                                                                                    >  7:119602/1‑149 (MQ=255)
              ccatgccatgCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAgc                                                                                                                                    >  7:227163/1‑149 (MQ=255)
                    catgccaCTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTg                                                                                                                              >  5:193153/1‑149 (MQ=255)
                    catgccaCTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTg                                                                                                                              >  3:233829/1‑149 (MQ=255)
                        ccaCTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAg                                                                                                                           <  6:193153/148‑1 (MQ=255)
                         caCTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAggg                                                                                                                         <  7:238609/149‑1 (MQ=255)
                             cGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGTGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGcagc                                                                                                                     <  8:275156/149‑1 (MQ=255)
                              ggTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCACGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGcag                                                                                                                      <  6:52421/147‑1 (MQ=255)
                               gtgtGCAGCCACAGGCTGTCGAGTTGATGCTGGGTCAGCGTGCAGGCAGGCACCGCTGTGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCATGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGcagcca                                                                                                                   <  2:116256/149‑1 (MQ=255)
                                  tgCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGcagccagc                                                                                                                 <  8:278460/148‑1 (MQ=255)
                                    cAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAAc                                                                                                              >  2:367650/1‑149 (MQ=255)
                                       ccACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTgg                                                                                                           <  8:119602/149‑1 (MQ=255)
                                        cacaGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCGCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTgg                                                                                                           >  1:355529/1‑148 (MQ=255)
                                        cacaGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGc                                                                                                          <  3:134525/149‑1 (MQ=255)
                                            ggCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGcagcca                                                                                                                   >  2:434235/1‑136 (MQ=255)
                                            ggCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGcagcca                                                                                                                   <  1:434235/136‑1 (MQ=255)
                                            ggCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGt                                                                                                      <  4:420901/149‑1 (MQ=255)
                                             gCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTg                                                                                                     >  2:33922/1‑149 (MQ=255)
                                             gCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTg                                                                                                     <  5:356460/149‑1 (MQ=255)
                                             gCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCCCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCTGTGCAGGGCAGCCAGCAACTGGCAGGTg                                                                                                     >  7:126076/1‑149 (MQ=255)
                                              cTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTgg                                                                                                    <  8:227163/149‑1 (MQ=255)
                                               tGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGc                                                                                                   >  2:153675/1‑149 (MQ=255)
                                               tGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGc                                                                                                   >  7:145887/1‑149 (MQ=255)
                                                 tCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAg                                                                                                                           >  3:185942/1‑123 (MQ=255)
                                                 tCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAg                                                                                                                           <  4:185942/123‑1 (MQ=255)
                                                 tCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCtt                                                                                                 <  2:80623/149‑1 (MQ=255)
                                                   gAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTg                                                                                                <  4:145806/148‑1 (MQ=255)
                                                   gAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGa                                                                                               <  7:291235/149‑1 (MQ=255)
                                                      ttGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCa                                                                                            >  5:57584/1‑149 (MQ=255)
                                                        gctgctGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCAt                                                                                           <  1:33922/148‑1 (MQ=255)
                                                          tgctgGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATc                                                                                          >  2:395864/1‑147 (MQ=255)
                                                            ctgGGGCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTAGCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCGGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaagg                                                                                      <  3:206182/149‑1 (MQ=255)
                                                                gTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaaggcc                                                                                    >  1:112737/1‑147 (MQ=255)
                                                                gTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaaggcc                                                                                    >  8:138459/1‑147 (MQ=255)
                                                                   aGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTg                                                                                                            <  2:421974/120‑1 (MQ=255)
                                                                   aGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTg                                                                                                            >  1:421974/1‑120 (MQ=255)
                                                                   aGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaaggccaaggc                                                                               <  2:355529/149‑1 (MQ=255)
                                                                     cGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCtt                                                                             >  6:252264/1‑149 (MQ=255)
                                                                         caggcaggCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATCCCCATCAGGGcc                                                                                                                                             <  5:202541/81‑1 (MQ=38)
                                                                         caggcaggCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATCCCCATCAGGGcc                                                                                                                                             >  6:202541/1‑81 (MQ=38)
                                                                              aggcaCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTg                                                                    >  6:50579/1‑149 (MQ=255)
                                                                              aggcaCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTg                                                                    <  8:391604/149‑1 (MQ=255)
                                                                                gcaCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGcc                                                                  >  6:118877/1‑149 (MQ=255)
                                                                                gcaCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGcc                                                                  >  7:96402/1‑149 (MQ=255)
                                                                                     gCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCTGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCgcag                                                             <  8:126076/149‑1 (MQ=255)
                                                                                                  tACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTgg                                                 >  3:222182/1‑148 (MQ=255)
                                                                                                    ccccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCa                                               <  5:50579/148‑1 (MQ=255)
                                                                                                    ccccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAg                                              <  1:153675/149‑1 (MQ=255)
                                                                                                    ccccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAg                                              >  7:175924/1‑149 (MQ=255)
                                                                                                           ccGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGgcg                                       <  8:96402/149‑1 (MQ=255)
                                                                                                           ccGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGgcg                                       <  2:55018/149‑1 (MQ=255)
                                                                                                             gcgTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGgcg                                       >  6:278334/1‑147 (MQ=255)
                                                                                                             gcgTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGgcg                                       <  5:278334/147‑1 (MQ=255)
                                                                                                               gTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGt                                   <  5:89040/149‑1 (MQ=255)
                                                                                                               gTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCCTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGt                                   <  7:138459/149‑1 (MQ=255)
                                                                                                                tCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTg                                  <  6:57584/149‑1 (MQ=255)
                                                                                                                      aGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCa                                     >  7:266057/1‑140 (MQ=255)
                                                                                                                      aGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCa                                     <  8:266057/140‑1 (MQ=255)
                                                                                                                      aGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTg                            >  3:365824/1‑149 (MQ=255)
                                                                                                                            ttGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGTAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTa                      >  8:229649/1‑149 (MQ=255)
                                                                                                                            ttGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTa                      >  8:15572/1‑149 (MQ=255)
                                                                                                                            ttGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTa                      >  2:57007/1‑149 (MQ=255)
                                                                                                                            ttGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTGGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTa                      >  1:154325/1‑149 (MQ=255)
                                                                                                                                gCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTcc                  >  3:336552/1‑149 (MQ=255)
                                                                                                                                ccccggCCCTGATCCCCATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTcc                  >  1:419301/5‑149 (MQ=255)
                                                                                                                                 cAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTccc                 >  7:284972/1‑149 (MQ=255)
                                                                                                                                  ccggCCCTGATCCCCATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTccc                 >  3:277405/3‑148 (MQ=255)
                                                                                                                                  ccggCCCTGATCCCCATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTccc                 <  4:277405/146‑1 (MQ=255)
                                                                                                                                   gggCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCg               >  5:420618/1‑149 (MQ=255)
                                                                                                                                    ggCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCgg              <  2:112737/149‑1 (MQ=255)
                                                                                                                                        cTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCAct          <  2:154325/149‑1 (MQ=255)
                                                                                                                                         tGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCActg         >  7:160080/1‑149 (MQ=255)
                                                                                                                                         tGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCActg         >  7:87599/1‑149 (MQ=255)
                                                                                                                                           aTGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCActgct       <  4:233829/149‑1 (MQ=255)
                                                                                                                                            tGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCActgctg      >  2:344404/1‑149 (MQ=255)
                                                                                                                                             ggggATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCActgctg      >  8:379877/1‑148 (MQ=255)
                                                                                                                                                gATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACt  >  5:242541/1‑149 (MQ=255)
                                                                                                                                                gATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACt  >  5:416420/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CTCGAAGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACT  >  NC_002947/460794‑461086

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: