Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,143,579 T→G 5.6% S29A (TCG→GCG) ‡ PP_3643 → oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,143,5790TG 5.6% 982.4 / 32.8 286S29A (TCG→GCG) ‡PP_3643oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (145/125);  new base G (7/9);  total (152/134)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

Only 100 of 293 total aligned reads displayed.

GTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGAACAGCAG  >  NC_002947/4143430‑4143728
                                                                                                                                                     |                                                                                                                                                     
 tataCCATCAGACACCTCCACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTc                                                                                                                                                      <  2:5633768/150‑1 (MQ=255)
   taCCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCgg                                                                                                                                                    <  2:5022398/150‑1 (MQ=255)
     ccATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGgtc                                                                                                                                                  <  1:3730090/150‑1 (MQ=255)
     ccATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGgtc                                                                                                                                                  <  1:1806909/150‑1 (MQ=255)
      cATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGgtcg                                                                                                                                                 >  1:3233451/1‑150 (MQ=255)
      cATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGgtcg                                                                                                                                                 >  1:2050022/1‑150 (MQ=255)
       aTCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGgtcgt                                                                                                                                                >  2:5268498/1‑150 (MQ=255)
       aTCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGgtcgt                                                                                                                                                >  1:288554/1‑150 (MQ=255)
        tCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGgtcgtc                                                                                                                                               <  2:4126606/150‑1 (MQ=255)
           gACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGc                                                                                                                                            <  1:93239/150‑1 (MQ=255)
           gACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGc                                                                                                                                            >  2:3210460/1‑150 (MQ=255)
           gACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGc                                                                                                                                            <  1:5407163/150‑1 (MQ=255)
            acacATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGcc                                                                                                                                           >  1:1379989/1‑150 (MQ=255)
            acacATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACCGAGCCGTCGTCGTcc                                                                                                                                           <  1:5229456/150‑1 (MQ=255)
              acaTCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTg                                                                                                                                         >  2:5104662/1‑150 (MQ=255)
               caTCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGa                                                                                                                                        <  2:2669747/150‑1 (MQ=255)
                 tCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACCAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGAcc                                                                                                                                      >  1:195740/1‑150 (MQ=255)
                 tCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGgt                                                                                                                                                   <  2:3937321/137‑1 (MQ=255)
                 tCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGAcc                                                                                                                                      >  2:5837109/1‑150 (MQ=255)
                   aaCGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGc                                                                                                                                    >  2:250511/1‑150 (MQ=255)
                    aCGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGcc                                                                                                                                   <  1:3503840/150‑1 (MQ=255)
                     cGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCGACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGCCGGCTCGGTCGTCGGCCTGACCGccc                                                                                                                                  >  1:3440441/1‑149 (MQ=255)
                           aaaCCGCGGACCACCGCAGCCCTTCATTGGTACTCCCTATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCt                                                                                                                            <  1:1885714/150‑1 (MQ=255)
                            aaCCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTATACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTAATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCtt                                                                                                                           <  1:5586439/150‑1 (MQ=255)
                               cgcgGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTtctc                                                                                                                        >  2:6382904/1‑150 (MQ=255)
                                 cgGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTctcct                                                                                                                      <  1:467169/150‑1 (MQ=255)
                                   gACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCtcctcc                                                                                                                    <  2:6791286/150‑1 (MQ=255)
                                    accaccGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTccc                                                                                                                   >  2:4351753/1‑150 (MQ=255)
                                    accaccGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTccc                                                                                                                   <  2:698898/150‑1 (MQ=255)
                                     ccaccGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCt                                                                                                                  >  1:7222712/1‑150 (MQ=255)
                                         cGCAGCCCTTCAGAGGTCCTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCg                                                                                                              >  2:7342991/1‑150 (MQ=255)
                                          gCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGa                                                                                                             <  2:4616627/150‑1 (MQ=255)
                                          gCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGa                                                                                                             <  1:6819830/150‑1 (MQ=255)
                                           cAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGAt                                                                                                            <  2:2937941/150‑1 (MQ=255)
                                           cAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGAt                                                                                                            <  1:40108/150‑1 (MQ=255)
                                            aGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATg                                                                                                           >  1:2644098/1‑150 (MQ=255)
                                            aGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATg                                                                                                           <  2:3532568/150‑1 (MQ=255)
                                              cccTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGa                                                                                                         <  1:2213436/150‑1 (MQ=255)
                                                 ttCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGAccc                                                                                                      >  2:4952469/1‑150 (MQ=255)
                                                   cAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACccgc                                                                                                    <  1:1542736/150‑1 (MQ=255)
                                                        gTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTg                                                                                               >  1:304565/1‑150 (MQ=255)
                                                          aCTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGt                                                                                             <  2:6159324/150‑1 (MQ=255)
                                                           cTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTg                                                                                            <  1:160564/150‑1 (MQ=255)
                                                           cTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTAATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTg                                                                                            <  2:1465061/150‑1 (MQ=255)
                                                            tCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGc                                                                                           >  1:6247052/1‑150 (MQ=255)
                                                              ccGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTg                                                                                         >  2:5857866/1‑150 (MQ=255)
                                                              ccGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTg                                                                                         >  2:375154/1‑150 (MQ=255)
                                                                gatgatCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGtt                                                                                       <  2:2778983/150‑1 (MQ=255)
                                                                  tgatCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCTCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCt                                                                                     <  2:1352543/150‑1 (MQ=255)
                                                                    atCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGc                                                                                   <  1:814392/150‑1 (MQ=255)
                                                                       gACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGcccc                                                                                <  2:304565/150‑1 (MQ=255)
                                                                         cGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGTCGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCaa                                                                              <  2:2170421/150‑1 (MQ=255)
                                                                          gCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAAc                                                                             <  2:1813218/150‑1 (MQ=255)
                                                                            aaCTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTa                                                                           >  2:7155689/1‑150 (MQ=255)
                                                                            aaCTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTa                                                                           >  2:4188189/1‑150 (MQ=255)
                                                                            aaCTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTa                                                                           <  1:1827983/150‑1 (MQ=255)
                                                                               tGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTacac                                                                        <  2:2726393/150‑1 (MQ=255)
                                                                                 tCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCt                                                                      >  1:2121538/1‑150 (MQ=255)
                                                                                   tACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTcc                                                                    >  1:6974722/1‑150 (MQ=255)
                                                                                   tACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCCGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTcc                                                                    >  2:6273900/1‑150 (MQ=255)
                                                                                    aCAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCg                                                                   >  2:1845625/1‑150 (MQ=255)
                                                                                      aaGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACAGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGAc                                                                 >  1:885986/1‑150 (MQ=255)
                                                                                       agcaCGTTCGGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACt                                                                <  1:6723270/147‑1 (MQ=255)
                                                                                        gAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTc                                                               <  2:2114589/150‑1 (MQ=255)
                                                                                        gAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTc                                                               >  2:7111799/1‑150 (MQ=255)
                                                                                               cTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCg                                                        <  1:2355806/150‑1 (MQ=255)
                                                                                                  ggCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTa                                                     >  2:3407880/1‑150 (MQ=255)
                                                                                                  ggCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTa                                                     <  2:7246502/150‑1 (MQ=255)
                                                                                                    cTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCCACTACACCTCCGACTCCTACCCGGTACt                                                   >  2:4081449/1‑150 (MQ=255)
                                                                                                     tCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTg                                                  >  2:2957048/1‑150 (MQ=255)
                                                                                                     tCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTg                                                  >  2:1604686/1‑150 (MQ=255)
                                                                                                     tCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTg                                                  >  2:819261/1‑150 (MQ=255)
                                                                                                       cTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGAt                                                >  2:192128/1‑150 (MQ=255)
                                                                                                       cTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGAt                                                >  1:640805/1‑150 (MQ=255)
                                                                                                       cTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGAt                                                <  1:2830153/150‑1 (MQ=255)
                                                                                                        ttCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCaa                                                                              <  2:2290576/119‑1 (MQ=255)
                                                                                                        ttCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTCCCCGGTCCTGATc                                               >  2:1124215/1‑150 (MQ=255)
                                                                                                         tCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCa                                              <  1:4395757/150‑1 (MQ=255)
                                                                                                          cccGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCaa                                             <  2:934512/150‑1 (MQ=255)
                                                                                                          cccGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCaa                                             <  1:6628189/150‑1 (MQ=255)
                                                                                                           ccGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAg                                            >  1:2563188/1‑150 (MQ=255)
                                                                                                            cGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGc                                           <  1:731155/150‑1 (MQ=255)
                                                                                                             gTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCa                                          >  1:948186/1‑150 (MQ=255)
                                                                                                             gTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCa                                          <  1:510023/150‑1 (MQ=255)
                                                                                                              tCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAg                                         >  2:6680384/1‑150 (MQ=255)
                                                                                                                       aCCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCttt                                <  1:2505867/150‑1 (MQ=255)
                                                                                                                        ccGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTg                               <  1:4952469/150‑1 (MQ=255)
                                                                                                                          gTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGcc                             >  1:1111353/1‑150 (MQ=255)
                                                                                                                           tcatcaCCACCAAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCACCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCa                            >  2:2586677/1‑150 (MQ=255)
                                                                                                                             atcaCCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCAtt                          <  2:3887238/150‑1 (MQ=255)
                                                                                                                               caccacCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCa                        >  2:3119042/1‑150 (MQ=255)
                                                                                                                                  caccCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTg                                                                                            <  2:5993764/79‑1 (MQ=255)
                                                                                                                                  caccCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTg                                                                                            >  1:5993764/1‑79 (MQ=255)
                                                                                                                                         gccgacgaccGAGCCGTCGTCGTCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTAc                                                    <  2:208127/110‑1 (MQ=25)
                                                                                                                                         gccgacgaccGAGCCGTCGTCGTCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTcc              >  1:1073112/3‑150 (MQ=37)
                                                                                                                                         gccgacgaccGAGCCGTCGTCGTCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACACGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTcc              >  2:4830922/3‑150 (MQ=37)
                                                                                                                                          acgacgacgGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCg             <  1:5705140/150‑1 (MQ=255)
                                                                                                                                           cgacgacgGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCgg            <  1:5837109/150‑1 (MQ=255)
                                                                                                                                             acgacgGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTg          <  2:2458525/150‑1 (MQ=255)
                                                                                                                                              cgacgGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTCTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGa         <  2:4466408/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTATACCATCAGACACATCAACGATAAAAACCGCGGACCACCGCAGCCCTTCAGAGGTACTCCCGATGATCGACGCAACTGTCTACAAGAACGTTCTGGGCTCCTTCCCGTCCGGCGTTACCGTCATCACCACCCAGGACGACGACGGCTCGGTCGTCGGCCTGACCGCCAGCGCCTTCTCCTCCCTGTCGATGGACCCGCCGCTGGTGCTGTTCTGCCCCAACTACACCTCCGACTCCTACCCGGTACTGATCAAGCAGAAGCGCTTTGCCATTCACCTGCTGTCCGGTGAACAGCAG  >  NC_002947/4143430‑4143728

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: ‑