Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,418,628 A→G A286A (GCA→GCG b3269 → putative transport system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,418,6280AG99.1% 656.1 / ‑3.0 213A286A (GCA→GCGb3269putative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (109/102);  total (110/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

Only 100 of 217 total aligned reads displayed.

TTTAATTCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCAGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCATTGCG  >  NC_000913/3418538‑3418717
                                                                                          |                                                                                         
tttAATTCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGt                                                                                         >  2:568404/1‑93 (MQ=255)
 ttAATTCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTg                                                                                        >  2:2173223/1‑93 (MQ=255)
 ttAATTCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTg                                                                                        <  2:3624530/93‑1 (MQ=255)
    aTTCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCt                                                                                     <  1:44061/93‑1 (MQ=255)
    aTTCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATACGCGCTGGGATCCAGGCGGTGCCt                                                                                     >  2:4070962/1‑93 (MQ=255)
     ttCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCtt                                                                                    <  2:959725/93‑1 (MQ=255)
      tCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTa                                                                                   >  1:4280835/1‑93 (MQ=255)
      tCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTa                                                                                   <  1:4080698/93‑1 (MQ=255)
        cTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATg                                                                                 >  2:1788033/1‑93 (MQ=255)
         tGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATgg                                                                                >  1:3064445/1‑93 (MQ=255)
         tGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATgg                                                                                <  2:3199586/93‑1 (MQ=255)
           aaCTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTc                                                                              >  2:3753885/1‑93 (MQ=255)
           aaCTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTc                                                                              >  1:2351978/1‑93 (MQ=255)
             cTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCaa                                                                            >  1:3209198/1‑93 (MQ=255)
             cTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCaa                                                                            >  1:4155837/1‑93 (MQ=255)
              tGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAAc                                                                           <  1:1010592/93‑1 (MQ=255)
               ggCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACa                                                                          <  1:830764/93‑1 (MQ=255)
                gCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACAt                                                                         <  2:1669873/93‑1 (MQ=255)
                gCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACAt                                                                         <  2:1223523/93‑1 (MQ=255)
                 cAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATg                                                                        >  2:4147940/1‑93 (MQ=255)
                  aGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGa                                                                       <  2:2408168/93‑1 (MQ=255)
                   gCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGaa                                                                      >  2:1958224/1‑93 (MQ=255)
                   gCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGaa                                                                      >  1:3228748/1‑93 (MQ=255)
                   gCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGaa                                                                      >  1:3027421/1‑93 (MQ=255)
                   gCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGaa                                                                      <  1:1585714/93‑1 (MQ=255)
                     cTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGc                                                                    <  2:351310/93‑1 (MQ=255)
                     cTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGc                                                                    <  2:4320716/93‑1 (MQ=255)
                       tAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCgg                                                                  <  2:4413317/93‑1 (MQ=255)
                          cGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTc                                                               <  1:3875178/93‑1 (MQ=255)
                            cTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCgg                                                             >  2:1185411/1‑93 (MQ=255)
                            cTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCgg                                                             >  2:4101677/1‑93 (MQ=255)
                             tGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGt                                                            <  1:1356210/93‑1 (MQ=255)
                              ggCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTc                                                           >  1:361453/1‑93 (MQ=255)
                              ggCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTc                                                           <  2:3246187/93‑1 (MQ=255)
                              ggCACTTTCGGTTTATACATCTGCATTTAGCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTc                                                           <  1:2705019/93‑1 (MQ=255)
                                cACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCAc                                                         >  1:1283077/1‑93 (MQ=255)
                                  cTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTccacg                                                       >  2:3999667/1‑89 (MQ=255)
                                  cTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTg                                                       >  1:1629537/1‑93 (MQ=255)
                                   tttCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTgg                                                      >  1:2900476/1‑93 (MQ=255)
                                    ttCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTggg                                                     <  2:4465371/93‑1 (MQ=255)
                                     tCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGa                                                    >  2:2150445/1‑93 (MQ=255)
                                     tCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGa                                                    <  1:559142/93‑1 (MQ=255)
                                        gTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTa                                                 >  1:4609266/1‑93 (MQ=255)
                                        gTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTa                                                 <  2:798520/93‑1 (MQ=255)
                                          ttATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTAcc                                               <  2:1543029/93‑1 (MQ=255)
                                          ttATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTAcc                                               <  2:1029442/93‑1 (MQ=255)
                                            ataCATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCa                                             <  1:917284/93‑1 (MQ=255)
                                              aCATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAAc                                           <  1:1045844/93‑1 (MQ=255)
                                               cATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAAcc                                          <  1:1929162/93‑1 (MQ=255)
                                               cATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAAcc                                          >  1:1195462/1‑93 (MQ=255)
                                               cATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACACAAcc                                          >  2:1469258/1‑93 (MQ=255)
                                                aTCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCGAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAAccc                                         <  2:1996361/93‑1 (MQ=255)
                                                 tCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCg                                        >  2:250094/1‑93 (MQ=255)
                                                 tCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCg                                        >  2:256200/1‑93 (MQ=255)
                                                 tCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCg                                        >  1:4040673/1‑93 (MQ=255)
                                                  cTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCgg                                       <  1:4133062/93‑1 (MQ=255)
                                                  cTGCATTTATCGCCGAGATCATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCgg                                       <  1:2930865/93‑1 (MQ=255)
                                                    gCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGtt                                     <  2:3031865/93‑1 (MQ=255)
                                                      aTTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAATATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtac                                   <  2:1562846/93‑1 (MQ=255)
                                                      aTTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAAACCGGTtac                                   >  1:1880111/1‑93 (MQ=255)
                                                       tttATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacg                                  >  2:1463532/1‑93 (MQ=255)
                                                        ttATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacgc                                 >  1:2855476/1‑93 (MQ=255)
                                                        ttATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGCGCCTTATGGTCAACATGAAGCGGCTCGGCCACTGGGATTACCCAACCCGGGtacgc                                 >  1:818045/1‑93 (MQ=255)
                                                         tATCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacgct                                >  2:1204568/1‑93 (MQ=255)
                                                          aTCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacgcta                               <  2:1254118/93‑1 (MQ=255)
                                                           tCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacgctac                              <  2:3894802/93‑1 (MQ=255)
                                                           tCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacgctac                              >  1:4673278/1‑93 (MQ=255)
                                                           tCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacgctac                              >  1:1659961/1‑93 (MQ=255)
                                                           tCGCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacgctac                              <  1:744963/93‑1 (MQ=255)
                                                             gCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacgctacgc                            >  2:2892841/1‑93 (MQ=255)
                                                             gCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacgctacgc                            >  1:3480250/1‑93 (MQ=255)
                                                             gCCGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTtacgctacgc                            <  1:1430007/93‑1 (MQ=255)
                                                              ccGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGcc                           <  2:2280272/93‑1 (MQ=255)
                                                               cGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCa                          <  1:2907139/93‑1 (MQ=255)
                                                               cGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCa                          <  1:2489413/93‑1 (MQ=255)
                                                               cGAGATTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGGTACGCTACGCCa                          >  2:4041258/1‑93 (MQ=255)
                                                               cGAGATTATCCGCGCTGGGATCCAGGCGGAGCCTTATGGTCAACATGAAGCGGCTCGGTTACTGGGATTACCCAACCCGGTTACGCTACGCCa                          >  1:6414/1‑93 (MQ=255)
                                                                gagaTTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAg                         >  2:692129/1‑93 (MQ=255)
                                                                gagaTTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAg                         >  2:3549479/1‑93 (MQ=255)
                                                                gagaTTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAg                         >  1:2645971/1‑93 (MQ=255)
                                                                 agaTTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAgg                        <  1:3267727/93‑1 (MQ=255)
                                                                  gaTTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGt                       >  2:4397661/1‑93 (MQ=255)
                                                                   aTTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTc                      <  1:2063662/93‑1 (MQ=255)
                                                                   aTTATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTc                      <  1:1254865/93‑1 (MQ=255)
                                                                    ttATCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCa                     >  2:2057949/1‑93 (MQ=255)
                                                                       tCCGCGCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCatta                  <  2:1732965/93‑1 (MQ=255)
                                                                         cgcgCTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCattatt                <  1:3800453/93‑1 (MQ=255)
                                                                          gcgcTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTc               <  1:4546426/93‑1 (MQ=255)
                                                                           cgcTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGGTACGCTACGCCAGGTCATTATTcc              >  2:4560853/1‑93 (MQ=255)
                                                                            gcTGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTccc             >  2:2810921/1‑93 (MQ=255)
                                                                              tGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCa           >  2:3568421/1‑93 (MQ=255)
                                                                              tGGGATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCa           >  2:3150641/1‑93 (MQ=255)
                                                                              tGGGATCCAGGCGGTCCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCa           <  1:1480056/93‑1 (MQ=255)
                                                                               gggATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAg          >  2:3326038/1‑93 (MQ=255)
                                                                                ggATCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCCGGCCATTATTCCCCAgg         >  2:870356/1‑93 (MQ=255)
                                                                                   tCCAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCAt      >  1:221506/1‑93 (MQ=255)
                                                                                    ccAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTCCCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCAtt     <  2:4618889/93‑1 (MQ=255)
                                                                                    ccAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCAtt     <  1:4559039/93‑1 (MQ=255)
                                                                                     cAGGCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCATTg    >  2:1345296/1‑93 (MQ=255)
                                                                                       ggCGGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGGCATTATTCCCCACGCATTTCg  >  1:1139297/1‑93 (MQ=255)
                                                                                          |                                                                                         
TTTAATTCCTGAACTGGCAGCCTTAACGCTGGCACTTTCGGTTTATACATCTGCATTTATCGCCGAGATTATCCGCGCTGGGATCCAGGCAGTGCCTTATGGTCAACATGAAGCGGCTCGGTCACTGGGATTACCCAACCCGGTTACGCTACGCCAGGTCATTATTCCCCAGGCATTGCG  >  NC_000913/3418538‑3418717

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