Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,934,461 A→G P186P (CCA→CCG rbsB → D‑ribose periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,934,4610AG98.9% 1461.7 / 1.9 458P186P (CCA→CCGrbsBD‑ribose periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base G (230/222);  total (234/224)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

Only 100 of 557 total aligned reads displayed.

CATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCAGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGTTCAGGCTGTATTCGCGCAGAATGA  >  NC_000913/3934371‑3934550
                                                                                          |                                                                                         
 aTTGCTGGTACATCCGCAGGCCGTGAACCTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCcggca                                                                                        <  2:5115954/93‑1 (MQ=255)
 aTTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCcggca                                                                                        <  1:2609985/93‑1 (MQ=255)
  ttGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCCGCcggcag                                                                                       >  1:8528472/1‑93 (MQ=255)
  ttGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCcggcag                                                                                       >  2:4758307/1‑93 (MQ=255)
   tGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGa                                                                                      <  2:6048694/93‑1 (MQ=255)
      tGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGAttt                                                                                   <  2:3610171/93‑1 (MQ=255)
       ggTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGAtttt                                                                                  <  2:5276831/93‑1 (MQ=255)
       ggTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGAtttt                                                                                  <  2:2768410/93‑1 (MQ=255)
         tACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGa                                                                                >  2:9868354/1‑93 (MQ=255)
          aCATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGAt                                                                               <  1:7987025/93‑1 (MQ=255)
           cATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATc                                                                              <  1:8715553/93‑1 (MQ=255)
           cATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATc                                                                              >  2:191945/1‑93 (MQ=255)
                gCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCAtt                                                                         <  1:5301229/93‑1 (MQ=255)
                gCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCAtt                                                                         <  2:1403113/93‑1 (MQ=255)
                gCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCAtt                                                                         >  2:3247802/1‑93 (MQ=255)
                gCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCAtt                                                                         >  2:3841730/1‑93 (MQ=255)
                   gcgcGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTaaa                                                                      <  1:2861921/90‑1 (MQ=255)
                   gCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTaaa                                                                      <  1:8295945/93‑1 (MQ=255)
                   gCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTaaa                                                                      >  2:8864797/1‑93 (MQ=255)
                   gCCCGTGAACGTGGCGAAGGCTTCCAGCACGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTaaa                                                                      >  2:2650187/1‑93 (MQ=255)
                    cccGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAg                                                                     <  2:8763490/93‑1 (MQ=255)
                      cGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGt                                                                   <  1:7838629/93‑1 (MQ=255)
                       gTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGtt                                                                  >  1:9364986/1‑93 (MQ=255)
                        tGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGttt                                                                 <  1:9580762/93‑1 (MQ=255)
                        tGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGttt                                                                 >  1:1837829/1‑93 (MQ=255)
                        tGAACGCGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGttt                                                                 >  1:6820015/1‑93 (MQ=255)
                          aaCGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGa                                                               <  1:4415402/93‑1 (MQ=255)
                               ggCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTa                                                          <  2:768805/93‑1 (MQ=255)
                               ggCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTa                                                          <  2:2466591/93‑1 (MQ=255)
                                gCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTaa                                                         <  2:669116/93‑1 (MQ=255)
                                gCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTaa                                                         <  1:4002192/93‑1 (MQ=255)
                                gCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTaa                                                         <  1:6198230/93‑1 (MQ=255)
                                   aaGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGc                                                      >  1:6654750/1‑93 (MQ=255)
                                   aaGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGc                                                      <  1:3348811/93‑1 (MQ=255)
                                   aaGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGc                                                      >  1:9240245/1‑93 (MQ=255)
                                   aaGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGc                                                      >  1:7593327/1‑93 (MQ=255)
                                    aGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCa                                                     >  1:2113891/1‑93 (MQ=255)
                                    aGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCa                                                     <  2:113504/93‑1 (MQ=255)
                                       cTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGaa                                                  >  1:4773475/1‑93 (MQ=255)
                                       cTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGaa                                                  <  2:7179515/93‑1 (MQ=255)
                                        ttCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAAc                                                 <  1:3036701/93‑1 (MQ=255)
                                         tCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAAcc                                                >  1:9683705/1‑93 (MQ=255)
                                          ccAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCt                                               >  1:4094257/1‑93 (MQ=255)
                                           cagcagGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTg                                              >  2:4358833/1‑93 (MQ=255)
                                             gcagGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGCACCTGtt                                            >  1:6440780/1‑93 (MQ=255)
                                               agGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGa                                          <  2:4988799/93‑1 (MQ=255)
                                               agGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGa                                          <  1:9195526/93‑1 (MQ=255)
                                                ggCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGAc                                         <  1:1359512/93‑1 (MQ=255)
                                                ggCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGAc                                         <  1:355413/93‑1 (MQ=255)
                                                ggCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGCTTGAACGTAATGCAGAACCTGTTGAc                                         >  1:1459550/1‑93 (MQ=255)
                                                 gCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGAcc                                        >  2:5964811/1‑93 (MQ=255)
                                                  ccGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCg                                       >  2:9498412/1‑93 (MQ=255)
                                                   cGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCAGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGc                                      >  1:2408927/1‑93 (MQ=255)
                                                    gTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCt                                     >  2:8613618/1‑93 (MQ=255)
                                                     ttGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTc                                    >  1:10017636/1‑93 (MQ=255)
                                                     ttGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTc                                    >  1:9083244/1‑93 (MQ=255)
                                                     ttGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTc                                    >  1:5809566/1‑93 (MQ=255)
                                                     ttGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTc                                    >  1:5034979/1‑93 (MQ=255)
                                                     ttGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTc                                    >  2:3107021/1‑93 (MQ=255)
                                                      tgctgcTCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCa                                   >  2:1039279/1‑93 (MQ=255)
                                                        ctgctCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATc                                 <  2:3426866/93‑1 (MQ=255)
                                                        ctgctCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCACATTTTGATCGCATTAAATGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATc                                 <  2:7535425/93‑1 (MQ=255)
                                                         tgctCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATcc                                <  1:4866155/93‑1 (MQ=255)
                                                            tCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGa                             <  2:5749141/93‑1 (MQ=255)
                                                            tCACAAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGa                             <  1:3661646/93‑1 (MQ=255)
                                                             cacaAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGAt                            >  1:7189728/1‑93 (MQ=255)
                                                              acaAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATg                           <  1:3213589/93‑1 (MQ=255)
                                                               caAGTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATgc                          >  2:2049449/1‑92 (MQ=255)
                                                                  gTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCa                       >  2:6872940/1‑93 (MQ=255)
                                                                  gTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCa                       <  1:430199/93‑1 (MQ=255)
                                                                  gTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCa                       >  2:8930187/1‑93 (MQ=255)
                                                                  gTTTAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCa                       >  2:2245525/1‑93 (MQ=255)
                                                                   tttAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAg                      >  2:806565/1‑93 (MQ=255)
                                                                   tttAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAg                      >  1:6960904/1‑93 (MQ=255)
                                                                    ttAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAgg                     <  1:1639209/93‑1 (MQ=255)
                                                                     tAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGc                    <  2:5416553/93‑1 (MQ=255)
                                                                     tAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGc                    >  1:8076362/1‑93 (MQ=255)
                                                                     tAATGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGc                    <  1:5942039/93‑1 (MQ=255)
                                                                      aaTGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCt                   >  2:1652970/1‑93 (MQ=255)
                                                                      aaTGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCt                   >  1:5432012/1‑93 (MQ=255)
                                                                        tGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGt                 >  2:7788711/1‑93 (MQ=255)
                                                                        tGTTCTTGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGt                 >  2:4030226/1‑93 (MQ=255)
                                                                             ttGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTcg            >  2:5629237/1‑93 (MQ=255)
                                                                             ttGCCAGCCAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTcg            <  2:4247265/93‑1 (MQ=255)
                                                                              tGCGAGCGAGCCGGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTcgc           <  1:1755419/93‑1 (MQ=255)
                                                                                ccagccagccggCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCgcgc         <  1:3897574/93‑1 (MQ=255)
                                                                                 cagccagccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCa        <  2:8405912/93‑1 (MQ=255)
                                                                                 cagccagccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCa        <  1:5602221/93‑1 (MQ=255)
                                                                                 cagccagccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCa        >  2:4637809/1‑93 (MQ=255)
                                                                                   gccagccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGa      >  1:7325831/1‑93 (MQ=255)
                                                                                   gccagccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGa      >  1:315852/1‑93 (MQ=255)
                                                                                    ccagccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGaa     <  1:8817613/93‑1 (MQ=255)
                                                                                     cagccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAat    <  2:6402733/93‑1 (MQ=255)
                                                                                      agccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGaacg   <  1:8977729/93‑3 (MQ=255)
                                                                                      agccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAatg   <  2:9054321/93‑1 (MQ=255)
                                                                                      agccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAatg   >  1:1181222/1‑93 (MQ=255)
                                                                                       gccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAatga  >  1:273452/1‑93 (MQ=255)
                                                                                       gccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAatga  >  2:6453380/1‑93 (MQ=255)
                                                                                       gccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAatga  >  2:925915/1‑93 (MQ=255)
                                                                                       gccggcAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGCTCAGGCTGTATTCGCGCAGAatga  <  2:9671879/93‑1 (MQ=255)
                                                                                          |                                                                                         
CATTGCTGGTACATCCGCAGCCCGTGAACGTGGCGAAGGCTTCCAGCAGGCCGTTGCTGCTCACAAGTTTAATGTTCTTGCCAGCCAGCCAGCAGATTTTGATCGCATTAAAGGTTTGAACGTAATGCAGAACCTGTTGACCGCTCATCCGGATGTTCAGGCTGTATTCGCGCAGAATGA  >  NC_000913/3934371‑3934550

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