Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 462,593 T→A 8.5% intergenic (‑129/+12) 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,5930TA 8.5% 231.0 / 10.2 82intergenic (‑129/+12)pqqF/davAcoenzyme PQQ synthesis protein F/5‑aminopentanamidase
Reads supporting (aligned to +/- strand):  ref base T (34/41);  new base A (4/3);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

GCGGATGGCGTCAGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCAGCAAG  >  NC_002947/462447‑462737
                                                                                                                                                  |                                                                                                                                                   
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                                                                                                                         >  2:112985/1‑149 (MQ=255)
                          cGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGt                                                                                                                         >  3:198236/1‑149 (MQ=255)
                            gTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTg                                                                                                                        >  4:184975/1‑148 (MQ=255)
                                      tcCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTg                                                                                                               <  5:548/147‑1 (MQ=255)
                                      tcCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTg                                                                                                               <  3:139506/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                                                                                                                  <  6:188519/139‑1 (MQ=255)
                                           gCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTc                                                                                                                  >  5:188519/1‑139 (MQ=255)
                                             tAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATc                                                                                                      >  1:118724/1‑149 (MQ=255)
                                              aaCCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCg                                                                                                    >  8:303559/1‑149 (MQ=255)
                                              aaCCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCg                                                                                                    >  6:309540/1‑149 (MQ=255)
                                              aaCCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCg                                                                                                    >  2:277891/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                                                                                                                                >  3:128358/1‑113 (MQ=255)
                                                       gtttgCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTa                                                                                                                                <  4:128358/113‑1 (MQ=255)
                                                          tgCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑ggg                                                                                       >  3:310473/1‑148 (MQ=255)
                                                           gCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGa                                                                                      >  3:279270/1‑148 (MQ=255)
                                                           gCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGa                                                                                      >  7:35917/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                                                                  >  5:389407/1‑149 (MQ=255)
                                                                              gggTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGcac                                                                  <  4:279270/149‑1 (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           <  8:183566/144‑1 (MQ=255)
                                                                                                                                          gggCCGCTATGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGTCCGGCCATg           >  7:183566/1‑144 (MQ=255)
                                                                                                                                          gggCCGCTATGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCAg      >  5:307564/1‑149 (MQ=255)
                                                                                                                                          gggCCGCTATGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCAg      >  5:302364/1‑149 (MQ=255)
                                                                                                                                                 aTAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCAGCAAg  <  1:368488/146‑1 (MQ=255)
                                                                                                                                                 aTAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCAGCAAg  >  2:368488/1‑146 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCGGATGGCGTCAGGCATGGGGGCTCCGGTGGCTGCTCTCCATGCTAACCGATTTGTTTGCAGGCCCGGCCTTTTCGCGGGTATGGCGCCGGTATGAAAGGCAGCGCAAATCCCGTGGGAGCGGGCGTGCCCGCGAATGGGCCGCATAGCGGCCCCTGTCAGCCTTTACGCAGGTGCAGCTCCTGGCGCAAATCGGTCAGGTA‑‑GGGAAACGCCCTGCGCCCCTGCACCACCCGCTCATGCTCGAGCTCTGCCAACAACTGGCACTCATCCCGCCCGGCCATGGCCAGCAAG  >  NC_002947/462447‑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: