Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 460,937 G→C 6.0% P510A (CCC→GCC) ‡ pqqF ← coenzyme PQQ synthesis protein F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947460,9370GC 6.0% 264.0 / 6.0 83P510A (CCC→GCC) ‡pqqFcoenzyme 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 = 4.93e-01

GTATCCTCGAAGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTG  >  NC_002947/460789‑461082
                                                                                                                                                    |                                                                                                                                                 
gTATCCTCGAAGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATggg                                                                                                                                                   <  3:410523/149‑1 (MQ=255)
     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:227163/1‑149 (MQ=255)
                   ccatgccatgCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAgc                                                                                                                                >  7:119602/1‑149 (MQ=255)
                         catgccaCTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTg                                                                                                                          >  3:233829/1‑149 (MQ=255)
                         catgccaCTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTg                                                                                                                          >  5:193153/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                                                                                               >  7:145887/1‑149 (MQ=255)
                                                    tGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGc                                                                                               >  2:153675/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                                                                                >  8:138459/1‑147 (MQ=255)
                                                                     gTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATcaaggcc                                                                                >  1:112737/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                                                                <  8:391604/149‑1 (MQ=255)
                                                                                   aggcaCCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTg                                                                >  6:50579/1‑149 (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                                          >  7:175924/1‑149 (MQ=255)
                                                                                                         ccccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAg                                          <  1:153675/149‑1 (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                  >  2:57007/1‑149 (MQ=255)
                                                                                                                                 ttGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTa                  >  8:15572/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:87599/1‑149 (MQ=255)
                                                                                                                                              tGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCActg     >  7:160080/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)
                                                                                                                                                    |                                                                                                                                                 
GTATCCTCGAAGCCTGCGGCCATGCCATGCCACTCGGTGTGCAGCCACAGGCTGTCGAGTTGCTGCTGGGTCAGCGTGCAGGCAGGCACCGCTGGGGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTG  >  NC_002947/460789‑461082

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: