Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,616,032 T→A 5.7% intergenic (+21/‑61) PP_mr40 → / → pet TPP regulatory region/protein Pet18

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,616,0320TA 5.7% 518.2 / 14.4 157intergenic (+21/‑61)PP_mr40/petTPP regulatory region/protein Pet18
Reads supporting (aligned to +/- strand):  ref base T (78/70);  new base A (5/4);  total (83/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 161 total aligned reads displayed.

TCGGGTATACCTGAGCGCCACCAAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGTGACCCGATTGTG  >  NC_002947/3615887‑3616181
                                                                                                                                                 |                                                                                                                                                     
tCGGGTATACCTGAGCGCCACCAAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGa                                                                                                                                                   <  1:4393561/150‑1 (MQ=255)
tCGGGTATACCTGAGCGCCACCAAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGa                                                                                                                                                   <  2:2333059/150‑1 (MQ=255)
tCGGGTATACCTGAGCGCCACCAAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGa                                                                                                                                                   >  1:93505/1‑150 (MQ=255)
   ggTATACCTGAGCGCCACCAAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGc                                                                                                                                                <  1:1232089/150‑1 (MQ=255)
       taCCTGAGCGCCACCAAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGt                                                                                                                                            <  1:625123/150‑1 (MQ=255)
       taCCTGAGCGCCACCAAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGt                                                                                                                                            <  1:625092/150‑1 (MQ=255)
        aCCTGAGCGCCACCAAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGTGCGGTGACCCTTCGAACCTGATCCGGAGCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCGCAGTGAGGccggca                                                                                                                                           <  2:2735445/150‑5 (MQ=255)
        aCCTGAGCGCCACCAAGGGGAGCCCGGCAAGGGGCTGAGAAACCGCTGAACGCTACAGTGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCGCAGTGAGGccggca                                                                                                                                           <  2:2594906/150‑5 (MQ=255)
         ccTGAGCGCCACCAAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTcc                                                                                                                                          >  2:1522199/1‑150 (MQ=255)
          cTGAGCGCCACCAAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCa                                                                                                                                         >  2:171717/1‑150 (MQ=255)
           tGAGCGCCACCAAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCAt                                                                                                                                        <  1:4771155/150‑1 (MQ=255)
            gAGCGCCACCAAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATg                                                                                                                                       >  1:1263042/1‑150 (MQ=255)
                 ccaccaAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAAt                                                                                                                                  <  2:4058094/150‑1 (MQ=255)
                 ccaccaAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAAt                                                                                                                                  <  1:2168980/150‑1 (MQ=255)
                  caccaAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAAtt                                                                                                                                 >  1:1400765/1‑150 (MQ=255)
                  caccaAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAAtt                                                                                                                                 <  1:4508619/150‑1 (MQ=255)
                   accaAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTg                                                                                                                                >  1:3175255/1‑150 (MQ=255)
                   accaAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTg                                                                                                                                >  2:1657823/1‑150 (MQ=255)
                     caAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGcc                                                                                                                              >  2:4351074/1‑150 (MQ=255)
                       aGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCaa                                                                                                                            <  1:4123580/150‑1 (MQ=255)
                       aGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCaa                                                                                                                            <  2:821476/150‑1 (MQ=255)
                                cGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCt                                                                                                                   >  1:3359830/1‑150 (MQ=255)
                                cGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCt                                                                                                                   >  1:3359641/1‑150 (MQ=255)
                                  gCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGt                                                                                                                 <  2:1532578/150‑1 (MQ=255)
                                  gCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGt                                                                                                                 <  2:3310524/150‑1 (MQ=255)
                                   cAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTc                                                                                                                <  1:2660492/150‑1 (MQ=255)
                                           tGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCtcat                                                                                                        >  2:1692376/1‑150 (MQ=255)
                                               aaaCCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAg                                                                                                    <  1:256339/150‑1 (MQ=255)
                                               aaaCCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAg                                                                                                    <  1:2104584/150‑1 (MQ=255)
                                                 aCCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGa                                                                                                  <  2:4542412/150‑1 (MQ=255)
                                                  ccGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGAc                                                                                                 <  1:2842128/150‑1 (MQ=255)
                                                     cTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGc                                                                                              <  2:1633130/150‑1 (MQ=255)
                                                      tGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCa                                                                                             >  2:1224176/1‑150 (MQ=255)
                                                         aCGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCAt                                                                                                                                        <  1:3205184/104‑1 (MQ=255)
                                                         aCGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCAt                                                                                                                                        >  2:3205184/1‑104 (MQ=255)
                                                            cTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGa                                                                                       <  2:1011783/150‑1 (MQ=255)
                                                             tACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGAc                                                                                      <  1:1050936/150‑1 (MQ=255)
                                                             tACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGAc                                                                                      >  2:3647639/1‑150 (MQ=255)
                                                             tACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCAATGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGAc                                                                                      >  1:693078/1‑150 (MQ=255)
                                                             tACAGCGCGGTGACCATTCGAAACTGATCCGGATCATGCCGGCGAAGCGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAACGACAGCAACATGAc                                                                                      <  2:3348665/150‑1 (MQ=255)
                                                                aGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGa                                                                                   >  2:4964927/1‑150 (MQ=255)
                                                                   gcgGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACg                                                                                >  1:1786377/1‑150 (MQ=255)
                                                                     ggTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGtt                                                                              <  1:4610527/150‑1 (MQ=255)
                                                                         aCCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAg                                                                          >  2:2510372/1‑150 (MQ=255)
                                                                         aCCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAg                                                                          >  2:1775165/1‑150 (MQ=255)
                                                                          cccTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGc                                                                         <  2:4371453/150‑1 (MQ=255)
                                                                          cccTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGc                                                                         <  1:885431/150‑1 (MQ=255)
                                                                           ccTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGcc                                                                        >  1:1654533/1‑150 (MQ=255)
                                                                             ttCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAg                                                                      >  2:2504617/1‑150 (MQ=255)
                                                                             ttCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAg                                                                      >  2:53784/1‑150 (MQ=255)
                                                                             ttCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAg                                                                      >  2:2056648/1‑150 (MQ=255)
                                                                             ttCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAg                                                                      >  2:4481869/1‑150 (MQ=255)
                                                                                gAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTcc                                                                   >  1:2742590/1‑150 (MQ=255)
                                                                                gAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTcc                                                                   >  2:5097571/1‑150 (MQ=255)
                                                                                 aaCCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTccc                                                                  >  1:4939859/1‑150 (MQ=255)
                                                                                   ccTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTg                                                                <  2:2424132/150‑1 (MQ=255)
                                                                                   ccTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTg                                                                <  1:4909935/150‑1 (MQ=255)
                                                                                     tGATCCGGATCATGCCGGCGAAGGGAAGGGGCTTGTAACCTGCCTATTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATAAAGGACAGCAACATGACCGAAAGTTTCAGCCAGTCCCTgcg                                                              >  2:4408091/1‑150 (MQ=255)
                                                                                       aTCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGcgcc                                                            <  2:722088/150‑1 (MQ=255)
                                                                                       aTCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGcgcc                                                            >  1:3922908/1‑150 (MQ=255)
                                                                                        tCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACCGCAACATGACCGAACGTTTCAGCCAGTCCCTGcgccg                                                           >  2:1169328/1‑150 (MQ=255)
                                                                                        tCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGcgccg                                                           >  2:1169373/1‑150 (MQ=255)
                                                                                         ccGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGcgccgc                                                          >  2:2066796/1‑150 (MQ=255)
                                                                                             aTCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCaga                                                      <  1:2905921/150‑1 (MQ=255)
                                                                                              tCATGCCGGCGACGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCACGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCagag                                                     <  1:3197135/150‑1 (MQ=255)
                                                                                               cATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAgagc                                                    <  2:833211/150‑1 (MQ=255)
                                                                                                   ccGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAgagcgagc                                                <  2:1989146/150‑1 (MQ=255)
                                                                                                         aaGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCt                                          <  2:898730/150‑1 (MQ=255)
                                                                                                         aaGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCt                                          >  2:3962504/1‑150 (MQ=255)
                                                                                                           gggATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCTCTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTgg                                        <  2:27026/150‑1 (MQ=255)
                                                                                                             gATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTc                                      <  1:2656848/150‑1 (MQ=255)
                                                                                                                 gggCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTcggcg                                  >  1:1965719/1‑150 (MQ=255)
                                                                                                                     ttGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCg                              >  1:2022440/1‑150 (MQ=255)
                                                                                                                       gTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGAc                                                                                                 >  2:2247642/1‑81 (MQ=255)
                                                                                                                       gTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGtt                            <  2:1523597/150‑1 (MQ=255)
                                                                                                                        tAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTg                           <  2:2359604/150‑1 (MQ=255)
                                                                                                                        tAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTg                           <  2:674578/150‑1 (MQ=255)
                                                                                                                         aaCCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTgg                          >  1:965980/1‑150 (MQ=255)
                                                                                                                            cTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCCGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCa                       >  2:4060652/1‑150 (MQ=255)
                                                                                                                            cTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCa                       >  1:4592852/1‑150 (MQ=255)
                                                                                                                             tGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGcc                               <  2:1135200/141‑1 (MQ=255)
                                                                                                                             tGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGcc                               >  1:1135200/1‑141 (MQ=255)
                                                                                                                              gCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCAcc                     <  2:2821203/150‑1 (MQ=255)
                                                                                                                              gCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCAcc                     >  1:4869405/1‑150 (MQ=255)
                                                                                                                              gCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCAcc                     <  1:397472/150‑1 (MQ=255)
                                                                                                                              gCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCAcc                     <  1:2549635/150‑1 (MQ=255)
                                                                                                                                cTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACTGCAACATGa                                                                                       <  1:1417383/82‑1 (MQ=255)
                                                                                                                                cTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACGGCAACATGa                                                                                       >  2:1417383/1‑82 (MQ=255)
                                                                                                                                 tttttGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAGGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGtt                  >  2:3163090/1‑150 (MQ=255)
                                                                                                                                  atggactggCCTCGCAGTGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGttt                 >  2:681932/8‑150 (MQ=255)
                                                                                                                                   tttGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTc                <  1:1885839/150‑1 (MQ=255)
                                                                                                                                   tttGCCGGCCTCACTGCGAGGCCAGTCCAGGTCAATTGCCAAGGTCTCGCTGTCACCCTCCTCACGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTc                <  1:3136733/150‑1 (MQ=255)
                                                                                                                                     tGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGt              <  1:2233307/150‑1 (MQ=255)
                                                                                                                                     gactggCCTCGCAGTGAGGCCGGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGt              >  1:3130551/5‑150 (MQ=255)
                                                                                                                                       ctggCCTCGCAGTGAGGCCGGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGTGa            >  2:342233/3‑150 (MQ=255)
                                                                                                                                            cTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGTGACCCGa       >  1:2785712/1‑150 (MQ=255)
                                                                                                                                             tCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGTGACCCGAt      >  1:3409104/1‑150 (MQ=255)
                                                                                                                                             tCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGTGACCCGAt      >  1:4798643/1‑150 (MQ=255)
                                                                                                                                                cTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGTGACCCGATtgt   <  2:2400639/150‑1 (MQ=255)
                                                                                                                                                 tGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGTGACCCGAtt     <  1:1060079/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TCGGGTATACCTGAGCGCCACCAAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGTGACCCGATTGTG  >  NC_002947/3615887‑3616181

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: