Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,160,285 C→G 5.3% G68R (GGG→CGG) ‡ 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,2850CG 5.3% 1027.4 / 33.6 302G68R (GGG→CGG) ‡PP_3658aromatic compound MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (154/132);  new base G (7/9);  total (161/141)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 310 total aligned reads displayed.

TAAAGCGCATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCATGAA  >  NC_002947/4160135‑4160434
                                                                                                                                                      |                                                                                                                                                     
    gcgcATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGaa                                                                                                                                                    <  1:846512/150‑1 (MQ=255)
      gcATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACa                                                                                                                                                  >  2:4784869/1‑150 (MQ=255)
        aTGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTTCGGGACAAt                                                                                                                                                <  2:6336918/150‑1 (MQ=255)
         tGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATg                                                                                                                                               >  2:4935439/1‑150 (MQ=255)
         tGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATg                                                                                                                                               >  1:3626996/1‑150 (MQ=255)
          gATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCTCAGCATTGTTCGGGACAATGc                                                                                                                                              >  1:4935048/1‑150 (MQ=255)
          gATACTGAAGCTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGc                                                                                                                                              <  1:5122671/150‑1 (MQ=255)
           aTACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCt                                                                                                                                             <  2:3696530/150‑1 (MQ=255)
            tACTGAACGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTTCGGGACAATGCTg                                                                                                                                            <  1:4739989/150‑1 (MQ=255)
              cTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTTCGGGACAATGCTggg                                                                                                                                          <  2:1955637/150‑3 (MQ=255)
               tGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGagc                                                                                                                                         >  1:6702475/1‑150 (MQ=255)
               tGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGagc                                                                                                                                         >  1:5961808/1‑150 (MQ=255)
                    gTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCTCAGCATTGTTCGGGACAATGCTGggcgctat                                                                                                                                    <  1:2001491/150‑6 (MQ=255)
                    gTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGagccagcc                                                                                                                                    <  2:6423789/150‑1 (MQ=255)
                     tCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAgccagcca                                                                                                                                   >  1:391320/1‑150 (MQ=255)
                         cGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTAc                                                                                                                               <  2:740478/150‑1 (MQ=255)
                              tGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTGTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCgg                                                                                                                          <  2:6599026/150‑1 (MQ=255)
                                  aGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCt                                                                                                                      >  1:1851018/1‑150 (MQ=255)
                                   gTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTg                                                                                                                     >  2:6500591/1‑150 (MQ=255)
                                   gTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTg                                                                                                                     <  2:300269/150‑1 (MQ=255)
                                   gTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTg                                                                                                                     >  1:5227028/1‑150 (MQ=255)
                                      aaCCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGtt                                                                                                                  >  1:916440/1‑150 (MQ=255)
                                      aaCCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGtt                                                                                                                  <  1:568094/150‑1 (MQ=255)
                                        cccTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGc                                                                                                                >  1:4010761/1‑150 (MQ=255)
                                        cccTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGc                                                                                                                >  1:591635/1‑150 (MQ=255)
                                          cTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGt                                                                                                              <  1:3820257/150‑1 (MQ=255)
                                                cAGCGGTGAAAACGCTGAACAGCGCTACACCGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAAttt                                                                                                        >  2:53665/1‑150 (MQ=255)
                                                   cGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCa                                                                                                     >  1:3014968/1‑150 (MQ=255)
                                                       gAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATc                                                                                                 >  2:6074930/1‑150 (MQ=255)
                                                        aaaaCGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCg                                                                                                >  1:6785138/1‑150 (MQ=255)
                                                           aCGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGcc                                                                                             <  1:3733528/150‑1 (MQ=255)
                                                                  aCAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGAtgg                                                                                      <  1:2805222/150‑1 (MQ=255)
                                                                   cAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGAtggc                                                                                     <  1:2648312/150‑1 (MQ=255)
                                                                    aGCGCTACACCGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGAtggct                                                                                    >  2:3100013/1‑150 (MQ=255)
                                                                          acacAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcagg                                                                              <  1:4983127/150‑1 (MQ=255)
                                                                           cacaGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcaggg                                                                             <  2:1483045/150‑1 (MQ=255)
                                                                           cacaGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcaggg                                                                             <  1:197948/150‑1 (MQ=255)
                                                                           cacaGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcaggg                                                                             <  1:3365679/150‑1 (MQ=255)
                                                                           cacaGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcaggg                                                                             >  1:1507999/1‑150 (MQ=255)
                                                                             caGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcagggca                                                                           <  2:1145662/150‑1 (MQ=255)
                                                                              aGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcagggcag                                                                          <  1:5372982/150‑1 (MQ=255)
                                                                              aGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTggcagggcag                                                                          >  1:6492303/1‑150 (MQ=255)
                                                                               gATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGc                                                                         >  2:6222061/1‑150 (MQ=255)
                                                                                 tCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGccc                                                                       >  2:4444317/1‑150 (MQ=255)
                                                                                   gcgATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGccccc                                                                     >  2:677964/1‑150 (MQ=255)
                                                                                     gATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCccac                                                                   >  2:2840300/1‑150 (MQ=255)
                                                                                     gATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCccac                                                                   >  2:5864146/1‑150 (MQ=255)
                                                                                      aTCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCccacc                                                                  >  1:1397554/1‑150 (MQ=255)
                                                                                       tCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCccacca                                                                 >  1:425672/1‑150 (MQ=255)
                                                                                         aTTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACt                                                               >  2:437066/1‑150 (MQ=255)
                                                                                           tttCGGCCGGCCAATTCGATCTGCCAAGGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGc                                                             <  1:4881612/150‑1 (MQ=255)
                                                                                                 ccggccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTc                                                       <  2:5645087/150‑1 (MQ=255)
                                                                                                 ccggccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTc                                                       <  2:1365497/150‑1 (MQ=255)
                                                                                                 ccggccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTc                                                       <  2:3251533/150‑1 (MQ=255)
                                                                                                    gccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTa                                                    <  2:1303692/150‑1 (MQ=255)
                                                                                                    gccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTa                                                    <  2:4011722/150‑1 (MQ=255)
                                                                                                     ccAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAg                                                   >  1:4925096/1‑150 (MQ=255)
                                                                                                       aaTTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGcc                                                 <  1:4068308/150‑1 (MQ=255)
                                                                                                         ttCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCCCCACTGCCAAGTCGTAGCCAt                                               <  2:2456411/150‑1 (MQ=255)
                                                                                                            gATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGa                                            <  1:2394753/150‑1 (MQ=255)
                                                                                                             aTCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGAt                                           >  2:2062714/1‑150 (MQ=255)
                                                                                                                  ccAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCa                                      <  2:5804651/150‑1 (MQ=255)
                                                                                                                  ccAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCa                                      <  2:2177953/150‑1 (MQ=255)
                                                                                                                    aaTGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGa                                    <  1:2159152/150‑1 (MQ=255)
                                                                                                                    aaTGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGa                                    >  2:3578070/1‑150 (MQ=255)
                                                                                                                         gCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCaa                               <  1:6276554/150‑1 (MQ=255)
                                                                                                                          ccTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAac                              >  1:3653349/1‑150 (MQ=255)
                                                                                                                          ccTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAac                              >  2:1229231/1‑150 (MQ=255)
                                                                                                                           cTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAaca                             >  1:4587694/1‑150 (MQ=255)
                                                                                                                           cTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAaca                             >  2:1788903/1‑150 (MQ=255)
                                                                                                                           cTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAaca                             >  2:4169298/1‑150 (MQ=255)
                                                                                                                           cTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAaca                             <  1:1067713/150‑1 (MQ=255)
                                                                                                                            tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTTGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                            >  2:5456274/1‑150 (MQ=255)
                                                                                                                            tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                            >  1:24353/1‑150 (MQ=255)
                                                                                                                            tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                            >  1:2510622/1‑150 (MQ=255)
                                                                                                                            tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                            >  2:4048731/1‑150 (MQ=255)
                                                                                                                            tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAacac                            >  1:1596888/1‑150 (MQ=255)
                                                                                                                            tAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAAGCAGAATCAacac                            >  2:5993375/1‑150 (MQ=255)
                                                                                                                              aaaaaaTAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGc                          >  2:183732/1‑150 (MQ=255)
                                                                                                                               aaaaaTAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCa                         <  1:1269568/150‑1 (MQ=255)
                                                                                                                               aaaaaTAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCa                         <  1:6539479/150‑1 (MQ=255)
                                                                                                                               aaaaaTAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCa                         >  2:5315357/1‑150 (MQ=255)
                                                                                                                                         ctcAGCATTGTTCGGGACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGACGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaataga               >  2:1031680/3‑150 (MQ=255)
                                                                                                                                          tcAGCATTGTTCGGGACAATGCTGGGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatagaa              >  2:5087117/2‑150 (MQ=255)
                                                                                                                                          tcAGCATTGTTCGGGACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatagaa              >  1:448959/2‑150 (MQ=255)
                                                                                                                                          tcAGCATTGTTCGGGACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatagaa              >  1:6095602/2‑150 (MQ=255)
                                                                                                                                          ccAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGTTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatagaa              <  2:6856573/150‑1 (MQ=255)
                                                                                                                                          ccAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatagaa              <  2:1707443/150‑1 (MQ=255)
                                                                                                                                          ccAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatagaa              <  1:3213419/150‑1 (MQ=255)
                                                                                                                                          ccAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatagaa              <  2:3405335/150‑1 (MQ=255)
                                                                                                                                           cAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCGCCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCagaatagaat             >  1:613344/1‑150 (MQ=255)
                                                                                                                                            aGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCcgaatagaaga            >  2:4595981/1‑148 (MQ=255)
                                                                                                                                             gCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATaa           <  1:2298896/150‑1 (MQ=255)
                                                                                                                                              cATTGTTCGGGACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAg          >  1:1798671/1‑150 (MQ=255)
                                                                                                                                              cATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAg          >  2:390906/1‑150 (MQ=255)
                                                                                                                                                 tGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTcc       <  1:751479/150‑1 (MQ=255)
                                                                                                                                                   tCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGTTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCAt     >  1:124428/1‑150 (MQ=255)
                                                                                                                                                   tCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCAt     >  1:3490101/1‑150 (MQ=255)
                                                                                                                                                   tCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCAt     <  1:2143761/150‑1 (MQ=255)
                                                                                                                                                      cGAACCATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCAGATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCAGCGATCATCAGAATCAACACGCACCCGAATAGAGTCAGTcccgggc  >  2:1605424/1‑145 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TAAAGCGCATGATACTGAAGGTCAGCGGGCTGTCAGTGAACCCTGCCGCAGCGGTGAAAACGCTGAACAGCGCTACACAGATCGCGATCATTTTCGGCCGGCCAATTCGATCTGCCAATGTGCCTAAAAAAATAGCGCCCAGCATTGTCCCGAACAATGCTGAGCCAGCCATTACCCCGGCGCTGGTTGCGTCAATTTGCATATCGGCCATGATGGCTGGCAGGGCAGCCCCCACCACTGCCAAGTCGTAGCCATCGATAATCAGAATCAACACGCACCAGAATAGAATAAGTCCATGAA  >  NC_002947/4160135‑4160434

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: