Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,646,913 C→T 5.1% intergenic (‑43/+23) PP_3211 ← / ← PP_3212 ABC transporter ATP‑binding protein/Rieske 2Fe‑2S family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,646,9130CT 5.1% 496.4 / 9.4 156intergenic (‑43/+23)PP_3211/PP_3212ABC transporter ATP‑binding protein/Rieske 2Fe‑2S family protein
Reads supporting (aligned to +/- strand):  ref base C (80/68);  new base T (5/3);  total (85/71)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 158 total aligned reads displayed.

TGTAGGTTTTTTCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGACACAG  >  NC_002947/3646765‑3647062
                                                                                                                                                    |                                                                                                                                                     
  tAGGTTTTTTCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCTCCGGCGCa                                                                                                                                                    >  2:5048890/1‑150 (MQ=255)
  tAGGTTTTTTCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCa                                                                                                                                                    >  2:1104627/1‑150 (MQ=255)
      ttttttCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGTGCGCCgg                                                                                                                                                <  2:5018835/150‑1 (MQ=255)
      ttttttCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCgg                                                                                                                                                <  2:2989439/150‑1 (MQ=255)
       tttttCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCggg                                                                                                                                               >  1:254078/1‑150 (MQ=255)
       tttttCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCggg                                                                                                                                               <  1:423240/150‑1 (MQ=255)
            ccACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCGGGACGTCGGCCTTGGATAGGGTCTGCATCAACAGAGCGCGCGCCCCGTGGGTCATGAGCGGTGGTGTACGCCCCGGTGCGCCGGGGcata                                                                                                                                          <  2:1523453/150‑4 (MQ=255)
             cACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTg                                                                                                                                         <  1:1696931/150‑1 (MQ=255)
             cACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTg                                                                                                                                         <  1:2303603/150‑1 (MQ=255)
              aCCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTgg                                                                                                                                        <  1:519099/150‑1 (MQ=255)
               ccTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTggg                                                                                                                                       >  1:169319/1‑150 (MQ=255)
               ccTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTggg                                                                                                                                       >  2:330717/1‑150 (MQ=255)
                  gATTGGCGAAGAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCaa                                                                                                                                    <  2:1247538/150‑1 (MQ=255)
                  gATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCaa                                                                                                                                    <  2:2650079/150‑1 (MQ=255)
                         gAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAg                                                                                                                             >  2:3367823/1‑150 (MQ=255)
                             aaaGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGc                                                                                                                         >  1:3803066/1‑150 (MQ=255)
                             aaaGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGc                                                                                                                         <  2:1164219/150‑1 (MQ=255)
                              aaGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCt                                                                                                                        <  1:1698946/150‑1 (MQ=255)
                                ggggTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGa                                                                                                                      <  2:164456/150‑1 (MQ=255)
                                   gtgtCCACACCGCTGGTGGTCGGCAGGGCCGGCGGCGGGACGTCGGCCGTGGGTAGGGTCTGCATCGACATCGCGCTCTCGCCGTGGGTCATGAGCGGGGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCACGTGTCCGTGGCTGCggc                                                                                                                   >  2:4449996/1‑149 (MQ=255)
                                    tgtCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGac                                                                                                                  >  2:105442/1‑149 (MQ=255)
                                      tCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAaca                                                                                                                >  1:2713263/1‑150 (MQ=255)
                                      tCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAaca                                                                                                                >  2:1494062/1‑150 (MQ=255)
                                         acacCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACtt                                                                                                             >  2:3110701/1‑150 (MQ=255)
                                         acacCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACtt                                                                                                             >  1:503496/1‑150 (MQ=255)
                                         acacCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACtt                                                                                                             <  2:3567582/150‑1 (MQ=255)
                                          cacCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTc                                                                                                            <  2:2150296/150‑1 (MQ=255)
                                           acCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCg                                                                                                           >  2:456482/1‑150 (MQ=255)
                                           acCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCg                                                                                                           <  2:4203226/150‑1 (MQ=255)
                                              gATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGctc                                                                                                        >  2:3856059/1‑150 (MQ=255)
                                               aTGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGctct                                                                                                       >  2:1056529/1‑150 (MQ=255)
                                                tGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCtctc                                                                                                      <  1:1117799/150‑1 (MQ=255)
                                                  gtgtTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGc                                                                                                    >  1:2979630/1‑150 (MQ=255)
                                                     ttCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCAt                                                                                                 <  2:1930424/150‑1 (MQ=255)
                                                       cGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATgg                                                                                               <  1:3375040/150‑1 (MQ=255)
                                                         gCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGc                                                                                             >  1:4547279/1‑150 (MQ=255)
                                                         gCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGc                                                                                             >  1:3139180/1‑150 (MQ=255)
                                                         gCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGc                                                                                             >  2:3570411/1‑150 (MQ=255)
                                                         gCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGc                                                                                             <  1:838421/150‑1 (MQ=255)
                                                          cAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCt                                                                                            >  1:2301642/1‑150 (MQ=255)
                                                               ccTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAg                                                                                       >  2:1653423/1‑150 (MQ=255)
                                                                 tGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCa                                                                                     <  1:1562409/150‑1 (MQ=255)
                                                                   cGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGa                                                                                   <  1:1557912/150‑1 (MQ=255)
                                                                     gCTGGACGTCGGCCTTGGATAGGGTCTGCATCTACATAGCGCTCTCCCCGTGGGTCATGTGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGAcc                                                                                 <  1:2863160/150‑1 (MQ=255)
                                                                     gCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGAcc                                                                                 >  2:4109465/1‑150 (MQ=255)
                                                                     gCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGAcc                                                                                 >  2:2776242/1‑150 (MQ=255)
                                                                      cTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGAccc                                                                                >  2:4568987/1‑150 (MQ=255)
                                                                        ggACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCag                                                                              <  2:199985/150‑1 (MQ=255)
                                                                         gACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCagt                                                                             >  1:5081861/1‑150 (MQ=255)
                                                                               ggCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGc                                                                       <  2:4603551/150‑1 (MQ=255)
                                                                                gCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGcc                                                                      <  1:3708682/150‑1 (MQ=255)
                                                                                   ttGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGAc                                                                   >  2:3498446/1‑150 (MQ=255)
                                                                                   ttGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGAc                                                                   >  1:4008232/1‑150 (MQ=255)
                                                                                     ggATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACgg                                                                 <  2:4389224/150‑1 (MQ=255)
                                                                                     ggATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACgg                                                                 <  2:1558439/150‑1 (MQ=255)
                                                                                      gATAGGGTCTGCATCAACATAGAGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGGAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGa                                                                >  1:2848139/1‑150 (MQ=255)
                                                                                       aTAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGAt                                                               >  2:3337864/1‑150 (MQ=255)
                                                                                            gTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATaacca                                                          >  2:1923562/1‑150 (MQ=255)
                                                                                                  aTCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCgg                                                    >  1:4716806/1‑150 (MQ=255)
                                                                                                  aTCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCgg                                                    >  2:2645099/1‑150 (MQ=255)
                                                                                                    cAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGc                                                  >  1:1541602/1‑150 (MQ=255)
                                                                                                       cATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTg                                               <  1:4759111/150‑1 (MQ=255)
                                                                                                        aTAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGt                                              >  1:616758/1‑150 (MQ=255)
                                                                                                         tAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTc                                             <  1:164532/150‑1 (MQ=255)
                                                                                                         tAGCGCTCTCCCCGTGGGACATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTc                                             <  2:4633223/150‑1 (MQ=255)
                                                                                                         tAGCGCGCCCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCTCACCGGGGCGTTTGGCAAGTGTCAGTGGCTGAGGAACACGTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTc                                             <  1:763261/150‑1 (MQ=255)
                                                                                                             gcgctcCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGAt                                         <  2:574672/147‑1 (MQ=255)
                                                                                                             gcgcgccccGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGAt                                         <  1:3467810/145‑1 (MQ=255)
                                                                                                             gcTCGCCCCGTGTGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGAt                                         <  1:3467813/150‑1 (MQ=255)
                                                                                                              ctctCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATg                                        >  1:4327166/1‑150 (MQ=255)
                                                                                                              ctctCCCCGTAGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATg                                        <  1:2786173/150‑1 (MQ=255)
                                                                                                               tctcCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATgg                                       >  1:4735408/1‑150 (MQ=255)
                                                                                                               tctcCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATgg                                       >  2:3882857/1‑150 (MQ=255)
                                                                                                                 tcCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGACTTCAAGc                                                                                      >  1:987842/1‑101 (MQ=255)
                                                                                                                   cccGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGtt                                                                         >  2:565814/1‑112 (MQ=255)
                                                                                                                   cccGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCAtg                                   >  1:4772123/1‑150 (MQ=255)
                                                                                                                   cccGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCAtg                                   >  2:749037/1‑150 (MQ=255)
                                                                                                                         ggTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGt                             >  1:2611241/1‑150 (MQ=255)
                                                                                                                          gTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGtt                            <  1:1408768/150‑1 (MQ=255)
                                                                                                                           tCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTc                           >  1:2751295/1‑150 (MQ=255)
                                                                                                                               gAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCAcc                       >  1:4752816/1‑150 (MQ=255)
                                                                                                                               gAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCAcc                       <  1:1186116/150‑1 (MQ=255)
                                                                                                                                        ccaacgCCCCGGTGCGCCGGGGCATTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGcggcgg              >  1:110074/6‑150 (MQ=255)
                                                                                                                                        ccaaTGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGcggcgg              >  2:960147/4‑150 (MQ=255)
                                                                                                                                         gtATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCggg             >  2:3893903/1‑150 (MQ=255)
                                                                                                                                          tATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGgt            <  2:4525731/150‑1 (MQ=255)
                                                                                                                                          tATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGgt            <  2:511346/150‑1 (MQ=255)
                                                                                                                                          tATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGgt            >  1:1015424/1‑150 (MQ=255)
                                                                                                                                            tGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGgtgt          >  1:4455118/1‑150 (MQ=255)
                                                                                                                                            cgCCCCGGTGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGgtgt          >  1:1081344/2‑150 (MQ=255)
                                                                                                                                             gCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTc         <  2:3803066/150‑1 (MQ=255)
                                                                                                                                             gCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTc         <  2:2713263/150‑1 (MQ=255)
                                                                                                                                             gCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTc         <  2:4772123/150‑1 (MQ=255)
                                                                                                                                              ccccGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCAt                                                                                                 <  2:1975016/61‑1 (MQ=255)
                                                                                                                                              ccccGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCAt                                                                                                 >  1:1975016/1‑61 (MQ=255)
                                                                                                                                              ccccGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCg        >  2:3970503/1‑150 (MQ=255)
                                                                                                                                               cccGGCGCGCCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGa       >  2:28850/1‑150 (MQ=255)
                                                                                                                                               cccGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCTGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGa       >  1:4645666/1‑150 (MQ=255)
                                                                                                                                               cccGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGa       >  2:2732891/1‑150 (MQ=255)
                                                                                                                                                ccGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGac      >  1:4077669/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGTAGGTTTTTTCCACCTGATTGGCGAACAAAGGGGTGTCCACACCGATGGTGTTCGGCAGGGCCTGCGGCTGGACGTCGGCCTTGGATAGGGTCTGCATCAACATAGCGCTCTCCCCGTGGGTCATGAGCGGTGGTGTATGCCCCGGCGCACCGGGGCGTTGGGCAAGTGTCAGTGGCTGAGGAACACTTCGCTCTCGCCATGGGCTTCAAGCAGACCCAGTAGTTGCCGACGGATAACCAACCCGGGCTTGTCCGATGGCATGTGCAGTTCCACCCTGCGGCGGGTGTCGACACAG  >  NC_002947/3646765‑3647062

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