Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,160,283 C→T 7.7% G68G (GGG→GGA) ‡ PP_3658 ← aromatic compound MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,160,2830CT 7.7% 713.0 / 33.8 235G68G (GGG→GGA) ‡PP_3658aromatic compound MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (117/99);  new base T (12/6);  total (129/106)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 238 total aligned reads displayed.

ATAAAGCGCATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCATG  >  NC_002947/4160134‑4160432
                                                                                                                                                     |                                                                                                                                                     
aTAAAGCGCATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTc                                                                                                                                                       <  2:5305665/150‑1 (MQ=255)
    aGCGCATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGa                                                                                                                                                   >  1:4856900/1‑150 (MQ=255)
     gcgcATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGaa                                                                                                                                                  <  2:209762/150‑1 (MQ=255)
     gcgcATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGaa                                                                                                                                                  >  1:1709801/1‑150 (MQ=255)
     gcgcATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAAGGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGaa                                                                                                                                                  <  2:1475181/150‑1 (MQ=255)
      cgcATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAAc                                                                                                                                                 >  1:1113701/1‑150 (MQ=255)
          tGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATg                                                                                                                                             >  1:3558881/1‑150 (MQ=255)
          tGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATg                                                                                                                                             <  2:2543625/150‑1 (MQ=255)
           gATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGc                                                                                                                                            >  2:5592913/1‑150 (MQ=255)
           gATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGc                                                                                                                                            >  1:213055/1‑150 (MQ=255)
            aTACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCGGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCt                                                                                                                                           <  2:5465052/150‑1 (MQ=255)
             tACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTg                                                                                                                                          >  1:4690530/1‑150 (MQ=255)
             tACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCAGTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTg                                                                                                                                          >  2:5154171/1‑150 (MQ=255)
              aCTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTTCGGGACAATGCTgg                                                                                                                                         <  2:2304969/150‑2 (MQ=255)
              aCTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGa                                                                                                                                         <  1:2563990/150‑1 (MQ=255)
              aCTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGa                                                                                                                                         <  1:3049587/150‑1 (MQ=255)
                tGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGagc                                                                                                                                       >  1:73153/1‑150 (MQ=255)
                    ggTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGagccagc                                                                                                                                   >  2:2644252/1‑150 (MQ=255)
                     gTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGagccagcc                                                                                                                                  >  1:3878541/1‑150 (MQ=255)
                      tCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAgccagcca                                                                                                                                 >  2:780488/1‑150 (MQ=255)
                         gCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTa                                                                                                                              >  2:1357405/1‑150 (MQ=255)
                         gCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTa                                                                                                                              >  1:3347795/1‑150 (MQ=255)
                             gCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCCATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTAcccc                                                                                                                          >  2:3526217/1‑150 (MQ=255)
                                    gTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTTCGGGACAATGCTgg                                                                                                                                         <  1:1519971/128‑2 (MQ=255)
                                    gTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTTCGGGACAATGCTgg                                                                                                                                         >  2:1519971/1‑127 (MQ=255)
                                    gTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTg                                                                                                                   >  2:4807864/1‑150 (MQ=255)
                                    gTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTg                                                                                                                   >  1:1262287/1‑150 (MQ=255)
                                     tGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTgg                                                                                                                  <  1:1977606/150‑1 (MQ=255)
                                     tGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTgg                                                                                                                  >  2:1726134/1‑150 (MQ=255)
                                       aaCCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGtt                                                                                                                >  2:833832/1‑150 (MQ=255)
                                       aaCCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGtt                                                                                                                <  1:816462/150‑1 (MQ=255)
                                         cccTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGc                                                                                                              <  1:96986/150‑1 (MQ=255)
                                          ccTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCg                                                                                                             <  2:4610574/150‑1 (MQ=255)
                                          ccTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCg                                                                                                             >  1:2047131/1‑150 (MQ=255)
                                              ccGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAAGTCGAACTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCaa                                                                                                         <  1:2890492/150‑1 (MQ=255)
                                                gCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGAGCTGGTTGCGTCAAtt                                                                                                       >  1:4554458/1‑150 (MQ=255)
                                                   gCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGc                                                                                                    <  1:4471899/150‑1 (MQ=255)
                                                   gCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGc                                                                                                    <  1:572313/150‑1 (MQ=255)
                                                          aaaCGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCgg                                                                                             <  2:5486008/150‑1 (MQ=255)
                                                           aaCGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGc                                                                                            >  1:1442436/1‑150 (MQ=255)
                                                               cTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCatg                                                                                        >  2:1272278/1‑150 (MQ=255)
                                                                  aaCAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCatgatg                                                                                     >  2:140041/1‑150 (MQ=255)
                                                                     atcgcTACACAGATCGCGATCATTTTCGGCCGCCCAATTCGATCTGCAAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGACATGAtggct                                                                                  <  1:3603006/148‑1 (MQ=255)
                                                                     aGCGCTACACAGATCGCGATCATTTTTGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGAtggct                                                                                  >  2:2679216/1‑150 (MQ=255)
                                                                     aGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGAtggct                                                                                  <  1:3238381/150‑1 (MQ=255)
                                                                     aGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGAtggct                                                                                  <  2:4690530/150‑1 (MQ=255)
                                                                         cTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCAGAGCCAGACATTACCCCGGCGCTGGTTGCGTCAATTTGCAAATCGGCCATGAAGGCTggca                                                                              >  2:153436/1‑150 (MQ=255)
                                                                          tACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcag                                                                             <  1:1146495/150‑1 (MQ=255)
                                                                          tACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcag                                                                             >  1:3964303/1‑150 (MQ=255)
                                                                          tACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcag                                                                             >  1:2401999/1‑150 (MQ=255)
                                                                          tACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcag                                                                             >  1:1913611/1‑150 (MQ=255)
                                                                            cacaGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcaggg                                                                           <  1:2915019/150‑1 (MQ=255)
                                                                            cacaGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcaggg                                                                           <  1:2213776/150‑1 (MQ=255)
                                                                               aGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcagggcag                                                                        >  2:4669961/1‑150 (MQ=255)
                                                                                 aTCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGcc                                                                      >  2:3010246/1‑150 (MQ=255)
                                                                                 aTCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGcc                                                                      <  1:3518834/150‑1 (MQ=255)
                                                                                    gcgATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGccccc                                                                   <  2:3964303/150‑1 (MQ=255)
                                                                                       aTCATTTTCGGCCGGCCAATTCGCTCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCccacc                                                                <  1:4951756/150‑1 (MQ=255)
                                                                                          aTTTTCGGCCGGCCCATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACt                                                             >  2:2825607/1‑150 (MQ=255)
                                                                                                  ccggccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCatgatg                                                                                     <  1:3698341/118‑1 (MQ=255)
                                                                                                   cggccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCg                                                    >  2:1340542/1‑150 (MQ=255)
                                                                                                   cggccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCg                                                    >  2:4506842/1‑150 (MQ=255)
                                                                                                     gccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGCTGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTa                                                  <  2:11263/150‑1 (MQ=255)
                                                                                                     gccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTa                                                  <  2:3875544/150‑1 (MQ=255)
                                                                                                        aaTTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGcc                                               >  2:438537/1‑150 (MQ=255)
                                                                                                        aaTTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGcc                                               >  1:5375875/1‑150 (MQ=255)
                                                                                                        aaTTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGcc                                               >  2:162681/1‑150 (MQ=255)
                                                                                                             gATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGa                                          <  1:408838/150‑1 (MQ=255)
                                                                                                             gATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCAGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGa                                          >  1:648278/1‑150 (MQ=255)
                                                                                                               tCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATa                                        >  2:5406373/1‑150 (MQ=255)
                                                                                                                 tGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAAt                                      >  1:4139330/1‑150 (MQ=255)
                                                                                                                 tGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAAt                                      >  2:416317/1‑150 (MQ=255)
                                                                                                                  gCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATc                                     <  1:3178633/150‑1 (MQ=255)
                                                                                                                  gCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATc                                     <  2:5096256/150‑1 (MQ=255)
                                                                                                                     aaTGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCACGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGa                                  <  1:4159608/150‑1 (MQ=255)
                                                                                                                      aTGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGaa                                 <  2:178661/150‑1 (MQ=255)
                                                                                                                      aTGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGaa                                 <  1:1616339/150‑1 (MQ=255)
                                                                                                                      aTGTGCCTAAAAAAATAGCGCCCAGCAATGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGaa                                 <  2:1706378/150‑1 (MQ=255)
                                                                                                                        gtgCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATc                               <  1:5170353/150‑1 (MQ=255)
                                                                                                                           ccTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAac                            <  1:5483655/150‑1 (MQ=255)
                                                                                                                           ccTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAac                            <  1:3893504/150‑1 (MQ=255)
                                                                                                                            cTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAaca                           >  1:4237178/1‑150 (MQ=255)
                                                                                                                            cTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAaca                           >  1:1780538/1‑150 (MQ=255)
                                                                                                                             tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                          >  1:346387/1‑150 (MQ=255)
                                                                                                                             tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                          >  1:4331264/1‑150 (MQ=255)
                                                                                                                             tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                          >  2:3631956/1‑150 (MQ=255)
                                                                                                                             tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                          >  1:557668/1‑150 (MQ=255)
                                                                                                                                    aTAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCag                   <  2:4856576/150‑1 (MQ=255)
                                                                                                                                     gctggctcAGCATTGTTCGGGACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCaga                  >  2:4486694/8‑150 (MQ=255)
                                                                                                                                        cgcCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaata               <  2:557668/150‑1 (MQ=255)
                                                                                                                                        cgcCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaata               <  2:2047131/150‑1 (MQ=255)
                                                                                                                                           tcAGCATTGTTCGGGACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatagaa            >  2:3503439/2‑150 (MQ=255)
                                                                                                                                           ccAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatagaa            >  2:5025026/1‑150 (MQ=255)
                                                                                                                                            cAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatagaat           >  2:3985758/1‑150 (MQ=255)
                                                                                                                                            cAGCATTGTCCCGAACAATGCTGAGCCAGCAATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGACGTAGCCATCGATAATCAGAATCAACACTCACCagaatagaat           >  2:3717450/1‑150 (MQ=255)
                                                                                                                                             aGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAgaatagaata          >  2:2536802/1‑150 (MQ=255)
                                                                                                                                              gCATTGTTCGGGACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATaa         >  2:2788298/1‑150 (MQ=255)
                                                                                                                                                   gTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCa    <  1:1357405/150‑1 (MQ=255)
                                                                                                                                                    tCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCAt   >  2:1925042/1‑150 (MQ=255)
                                                                                                                                                     cccGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCATg  >  2:4909898/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATAAAGCGCATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCATG  >  NC_002947/4160134‑4160432

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: