Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,160,287 A→G 7.5% F67S (TTC→TCC)  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,2870AG 7.5% 745.9 / 32.7 240F67S (TTC→TCC) PP_3658aromatic compound MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (122/99);  new base G (12/6);  total (134/106)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 248 total aligned reads displayed.

AGCGCATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCATGAAAG  >  NC_002947/4160138‑4160436
                                                                                                                                                     |                                                                                                                                                     
aGCGCATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGa                                                                                                                                                       >  1:4856900/1‑150 (MQ=255)
 gcgcATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGaa                                                                                                                                                      <  2:209762/150‑1 (MQ=255)
     aTGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAAt                                                                                                                                                  >  2:1666738/1‑150 (MQ=255)
      tGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATg                                                                                                                                                 >  1:3558881/1‑150 (MQ=255)
      tGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATg                                                                                                                                                 <  1:2376316/150‑1 (MQ=255)
      tGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATg                                                                                                                                                 <  2:2543600/150‑1 (MQ=255)
       gATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGc                                                                                                                                                >  2:5592913/1‑150 (MQ=255)
        aTACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCGGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCt                                                                                                                                               <  2:5465052/150‑1 (MQ=255)
          aCTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTTCGGGACAATGCTgg                                                                                                                                             <  1:1252042/150‑2 (MQ=255)
           cTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCTCAGCATTGTTCGGGACAATGCTggg                                                                                                                                            <  2:3932277/150‑3 (MQ=255)
            tGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGagc                                                                                                                                           >  1:73153/1‑150 (MQ=255)
             gAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGagcc                                                                                                                                          >  2:5483655/1‑150 (MQ=255)
                 gTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGagccagcc                                                                                                                                      >  1:3878541/1‑150 (MQ=255)
                  tCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAgccagcca                                                                                                                                     >  2:3602390/1‑150 (MQ=255)
                  tCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAgccagcca                                                                                                                                     >  2:780488/1‑150 (MQ=255)
                         gCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCCATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTAcccc                                                                                                                              >  2:3526217/1‑150 (MQ=255)
                               aGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCt                                                                                                                        >  1:4081648/1‑150 (MQ=255)
                                gTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTTCGGGACAATGCTgg                                                                                                                                             <  1:1519971/128‑2 (MQ=255)
                                gTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTTCGGGACAATGCTgg                                                                                                                                             >  2:1519971/1‑127 (MQ=255)
                                 tGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTgg                                                                                                                      >  2:1725864/1‑150 (MQ=255)
                                 tGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTgg                                                                                                                      <  1:3356148/150‑1 (MQ=255)
                                   aaCCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGtt                                                                                                                    <  1:816462/150‑1 (MQ=255)
                                   aaCCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGtt                                                                                                                    >  2:833832/1‑150 (MQ=255)
                                    aCCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTg                                                                                                                   <  2:122878/150‑1 (MQ=255)
                                     cccTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGc                                                                                                                  <  1:96986/150‑1 (MQ=255)
                                      ccTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCg                                                                                                                 >  2:3374513/1‑150 (MQ=255)
                                          ccGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCaa                                                                                                             <  1:1347656/150‑1 (MQ=255)
                                            gCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGAGCTGGTTGCGTCAAtt                                                                                                           >  1:4554458/1‑150 (MQ=255)
                                             cAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAAttt                                                                                                          <  1:4156824/150‑1 (MQ=255)
                                              aGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTg                                                                                                         <  2:2587871/150‑1 (MQ=255)
                                               gCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGc                                                                                                        <  2:4620426/150‑1 (MQ=255)
                                               gCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGc                                                                                                        <  1:1942153/150‑1 (MQ=255)
                                                    gAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATc                                                                                                   >  1:3304258/1‑150 (MQ=255)
                                                     aaaaCGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCg                                                                                                  >  2:4670381/1‑150 (MQ=255)
                                                       aaCGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGc                                                                                                <  2:1948343/150‑1 (MQ=255)
                                                       aaCGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGc                                                                                                >  1:1442269/1‑150 (MQ=255)
                                                       aaCGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGc                                                                                                >  2:2659135/1‑150 (MQ=255)
                                                          gCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGACCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCat                                                                                             >  1:4612294/1‑150 (MQ=255)
                                                             gAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCatgat                                                                                          >  1:245625/1‑150 (MQ=255)
                                                              aaCAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGa                                                                                                                                             <  1:1521450/98‑1 (MQ=255)
                                                              aaCAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGa                                                                                                                                             >  2:1521450/1‑98 (MQ=255)
                                                                 atcgcTACACAGATCGCGATCATTTTCGGCCGCCCAATTCGATCTGCAAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGACATGAtggct                                                                                      <  1:3603006/148‑1 (MQ=255)
                                                                     cTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggca                                                                                  <  1:4228633/150‑1 (MQ=255)
                                                                      tACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcag                                                                                 >  1:2126728/1‑150 (MQ=255)
                                                                      tACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcag                                                                                 >  2:4813216/1‑150 (MQ=255)
                                                                        cacaGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcaggg                                                                               <  1:2914908/150‑1 (MQ=255)
                                                                        cacaGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcaggg                                                                               <  1:2213776/150‑1 (MQ=255)
                                                                        cacaGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcaggg                                                                               <  2:1423182/150‑1 (MQ=255)
                                                                             aTCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGcc                                                                          >  2:3010246/1‑150 (MQ=255)
                                                                             aTCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGcc                                                                          <  1:3518834/150‑1 (MQ=255)
                                                                                  gATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCccac                                                                     >  1:2041479/1‑150 (MQ=255)
                                                                                   aTCATTTTCGGCCGGCCAATTCGCTCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCccacc                                                                    <  1:4951756/150‑1 (MQ=255)
                                                                                        tttCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGc                                                               <  1:4383506/150‑1 (MQ=255)
                                                                                           cggcctgcCAAGTCGATCTGCCAATGTGCATATAAAATTAGCGCCCAGCATTGTCCCGATCAATGCTGAGCCAGCCATTACCCCTGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCaa                                                            <  2:5398461/150‑1 (MQ=255)
                                                                                              ccggccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCatgatg                                                                                         <  1:3698341/118‑1 (MQ=255)
                                                                                              cccccaaaTTCTAACTCCCAATGTGCCAAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTc                                                         <  2:5078199/144‑1 (MQ=255)
                                                                                               cggccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCg                                                        >  2:1340632/1‑150 (MQ=255)
                                                                                               cggccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCg                                                        >  2:4506842/1‑150 (MQ=255)
                                                                                                 gccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTa                                                      <  2:3875544/150‑1 (MQ=255)
                                                                                                 gccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTa                                                      <  2:512459/150‑1 (MQ=255)
                                                                                                   cAATTCGATCTGCCTATGTGCCTAAAAAAATAGCGCCCAGCATTGTGCGGAACTATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTTGCAGGGCAGCCCCCACCACTGGCTAGTCGTatc                                                    <  2:4317703/150‑3 (MQ=255)
                                                                                                    aaTTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGcc                                                   >  2:438537/1‑150 (MQ=255)
                                                                                                    aaTTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGcc                                                   >  1:5375875/1‑150 (MQ=255)
                                                                                                      ttCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCAt                                                 <  2:1277673/150‑1 (MQ=255)
                                                                                                         gATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCAGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGa                                              >  1:648278/1‑150 (MQ=255)
                                                                                                              gCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATc                                         <  1:3178633/150‑1 (MQ=255)
                                                                                                                 aaTGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGa                                      >  2:1746518/1‑150 (MQ=255)
                                                                                                                  aTGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGaa                                     <  2:1389628/150‑1 (MQ=255)
                                                                                                                    gtgCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATc                                   <  1:5170353/150‑1 (MQ=255)
                                                                                                                       ccTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAac                                >  1:640643/1‑150 (MQ=255)
                                                                                                                       ccTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAac                                <  1:3893504/150‑1 (MQ=255)
                                                                                                                        cTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGGTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAaca                               >  1:4090581/1‑150 (MQ=255)
                                                                                                                         tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                              >  1:346387/1‑150 (MQ=255)
                                                                                                                         tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                              >  2:4696055/1‑150 (MQ=255)
                                                                                                                         tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                              >  2:4036411/1‑150 (MQ=255)
                                                                                                                         tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                              >  1:4331264/1‑150 (MQ=255)
                                                                                                                         tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                              >  1:557668/1‑150 (MQ=255)
                                                                                                                         tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACGGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                              >  2:4998785/1‑150 (MQ=255)
                                                                                                                             aaaaTAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCAc                          >  1:1126319/1‑150 (MQ=255)
                                                                                                                               aaTAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCa                        <  1:1666738/150‑1 (MQ=255)
                                                                                                                                aTAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCag                       <  2:4856900/150‑1 (MQ=255)
                                                                                                                                 gctggctcAGCATTGTTCGGGACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCaga                      >  2:4486697/8‑150 (MQ=255)
                                                                                                                                    cgcCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaata                   <  1:1725864/150‑1 (MQ=255)
                                                                                                                                    cgcCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaata                   <  2:557668/150‑1 (MQ=255)
                                                                                                                                     gcCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCAACGATAATCAGAATCAACACGCACCagaatag                  >  1:1884674/1‑150 (MQ=255)
                                                                                                                                      cccAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaataga                 >  1:2628421/1‑150 (MQ=255)
                                                                                                                                        cAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatagaat               >  2:3985758/1‑150 (MQ=255)
                                                                                                                                        cAGCATTGTCCCGAACAATGCTGAGCCAGCAATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGACGTAGCCATCGATAATCAGAATCAACACTCACCagaatagaat               >  2:3717450/1‑150 (MQ=255)
                                                                                                                                         aGCATTGTTCGGGACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAgaatagaata              >  1:5495733/1‑150 (MQ=255)
                                                                                                                                         aGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAgaatagaata              >  2:2536802/1‑150 (MQ=255)
                                                                                                                                          gCATTGTTCGGGACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATaa             >  2:2788283/1‑150 (MQ=255)
                                                                                                                                          gCATTGTTCGGGACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATaa             >  2:2788329/1‑150 (MQ=255)
                                                                                                                                          gCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATaa             <  2:1889932/150‑1 (MQ=255)
                                                                                                                                             ttGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTc          >  1:3069968/1‑150 (MQ=255)
                                                                                                                                              tGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTcc         <  2:2300128/150‑1 (MQ=255)
                                                                                                                                                  ccGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCATGa     >  1:1100855/1‑150 (MQ=255)
                                                                                                                                                  ccGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCATGa     <  2:1113701/150‑1 (MQ=255)
                                                                                                                                                  ccGAACAATGCTGAGACAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACAACTGCCAAGTCGTAGCCATCTATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCATGa     >  2:1607244/1‑150 (MQ=255)
                                                                                                                                                   cGAACAATGCTGAGCCAGCCATTACCCCGGTGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCATGaa    >  1:2929783/1‑150 (MQ=255)
                                                                                                                                                   cGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCATGaa    >  2:52501/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCGCATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCATGAAAG  >  NC_002947/4160138‑4160436

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