Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,994,290 (T)5→4 100% pseudogene (352/444 nt) pbl → pseudogene, peptidoglycan‑binding enzyme family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,994,2860T.76.1% 454.5 / 216.7 205pseudogene (348/444 nt)pblpseudogene, peptidoglycan‑binding enzyme family
Reads supporting (aligned to +/- strand):  ref base T (18/30);  new base . (74/82);  total (93/112)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 227 total aligned reads displayed.

CCCCTGTTAAAAAAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATAT  >  NC_000913/2994188‑2994385
                                                                                                  |                                                                                                   
cgcttGTTCAAAAAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTGTTCTCAACAATAAAATAGGAACCGAGATTTTATCTAACCAtt                                                                                                    >  2:1181471/5‑100 (MQ=255)
ccccTGTTAAAAAAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCCtt                                                                                                    >  2:1736084/1‑100 (MQ=255)
 cccTGTTAAAAAAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCAttt                                                                                                   <  2:1218627/100‑1 (MQ=255)
  ccTGTTAAAAAAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCAtttt                                                                                                  >  1:1534161/1‑100 (MQ=255)
    tGTTAAAAAAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCAttttca                                                                                                >  1:2804001/1‑98 (MQ=255)
    tGTTAAAAAAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTc                                                                                                <  1:1580826/100‑1 (MQ=255)
     gTTAAAAAAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCAttttcac                                                                                               >  2:182803/1‑97 (MQ=255)
     gTTAAAAAAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCAttttcac                                                                                               >  1:3160477/1‑97 (MQ=255)
     gTTAAAAAAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCAttttcac                                                                                               <  2:3046789/100‑4 (MQ=255)
     gTTAAAAAAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCAttttcac                                                                                               >  1:1248328/1‑97 (MQ=255)
     gTTAAAAAAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCAttttc                                                                                                 >  2:2528806/1‑97 (MQ=255)
      ttAAAAAAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTc                                                                                                <  2:3955465/98‑1 (MQ=255)
       tAAAAAAACTTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTTAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACg                                                                                             <  1:344516/100‑1 (MQ=255)
       tAAAAAAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCAttttcac                                                                                               >  2:3978380/1‑95 (MQ=255)
        aaaaaaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGtt                                                                                           <  2:2339384/100‑1 (MQ=255)
         aaaaaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTtg                                                                                          >  1:456897/1‑100 (MQ=255)
         aaaaaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTtg                                                                                          >  2:1640435/1‑100 (MQ=255)
          aaaaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTtg                                                                                          >  2:3179496/1‑100 (MQ=255)
          aaaaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTtgt                                                                                         <  1:2044493/100‑1 (MQ=255)
           aaaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTtgtg                                                                                        >  2:378708/1‑100 (MQ=255)
           aaaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTtgtg                                                                                        <  2:1907694/100‑1 (MQ=255)
            aaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTgg                                                                                       <  2:2335097/100‑1 (MQ=255)
            aaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTgg                                                                                       <  1:838996/100‑1 (MQ=255)
                  ggAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTaa                                                                                  >  2:130007/1‑100 (MQ=255)
                  ggAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAAc                                                                                 >  2:2384192/1‑100 (MQ=255)
                   gAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACg                                                                                <  1:3240713/100‑1 (MQ=255)
                    aataatAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGt                                                                               >  1:3660829/1‑100 (MQ=255)
                     ataataAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTg                                                                              >  1:4172106/1‑100 (MQ=255)
                     ataataAAAGATAAACGAGATGTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTa                                                                                   <  2:2146317/95‑1 (MQ=255)
                       ttccgatcctGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTgg                                                                             <  1:3066949/91‑1 (MQ=255)
                       aataaAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGc                                                                            >  1:3150448/1‑100 (MQ=255)
                       aataaAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGc                                                                            <  1:37511/100‑1 (MQ=255)
                       aataaAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGc                                                                            >  2:797509/1‑100 (MQ=255)
                          aaaaGATGAACGAGATTTGTTAGATAACCTTTGTCTGACCATTAAAATAAGAACCGAAATTCTATTTAACTATTTTAAAGTTGTGGAGGAAcaacacaat                                                                         >  2:4043052/1‑92 (MQ=255)
                           aaaGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATg                                                                        >  1:3350284/1‑100 (MQ=255)
                                tGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTggagatcggaa                                                                    >  2:3066931/1‑91 (MQ=255)
                                 gAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTgg                                                                   <  1:2129360/100‑1 (MQ=255)
                                   aCGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGa                                                                <  1:2969083/100‑1 (MQ=255)
                                      agaTTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACtt                                                             <  2:2752201/100‑1 (MQ=255)
                                      agaTTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACt                                                              >  2:3446145/1‑99 (MQ=255)
                                      agaTTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCNGGGGACt                                                              <  2:250558/99‑1 (MQ=255)
                                       gaTTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTta                                                            <  1:1887227/100‑1 (MQ=255)
                                       gaTTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTta                                                            >  1:690960/1‑100 (MQ=255)
                                         tttGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACTGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTTACAATGCCTGGGGACTtata                                                          >  2:2871618/1‑100 (MQ=255)
                                          ttGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATaa                                                         <  2:3660802/100‑1 (MQ=255)
                                            gTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATg                                                       >  1:2140417/1‑100 (MQ=255)
                                               aGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGAGTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCg                                                     >  1:224168/1‑99 (MQ=255)
                                                gATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCggg                                                   <  1:1165437/100‑1 (MQ=255)
                                                gATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCggg                                                   >  1:3191333/1‑100 (MQ=255)
                                                 aTAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCggg                                                   <  2:2117779/100‑1 (MQ=255)
                                                  tAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGc                                                      >  2:2565182/1‑95 (MQ=255)
                                                  tAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGt                                                  >  1:567695/1‑99 (MQ=255)
                                                     cccTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGtttt                                               <  2:3501305/99‑1 (MQ=255)
                                                      ccTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGc                                             <  1:4126914/100‑1 (MQ=255)
                                                      ccTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGc                                             <  1:3523083/100‑1 (MQ=255)
                                                       cTTGTCTGAACATTAACATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTg                                              <  1:3610347/98‑1 (MQ=255)
                                                       cTTGTCTGAACATTAAAATAGGATCCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCa                                            <  1:1740064/100‑1 (MQ=255)
                                                       cTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGc                                             <  2:3215918/100‑1 (MQ=255)
                                                          gTCTGAACATTAAAATAGGAACCGAGATTTTATATATCCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATg                                         <  1:2477179/100‑1 (MQ=255)
                                                          gTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATg                                         >  1:1466225/1‑100 (MQ=255)
                                                          gTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATg                                         <  1:2469000/100‑1 (MQ=255)
                                                          gTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATg                                         >  1:4072763/1‑100 (MQ=255)
                                                          gTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATg                                         <  2:1022442/100‑1 (MQ=255)
                                                           tCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATgg                                        <  2:1356751/100‑1 (MQ=255)
                                                           tCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATgg                                        <  1:331129/100‑1 (MQ=255)
                                                            cTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTTGGGACTTATAATGCGGGATTTGAAACGGa                                       >  2:690342/1‑100 (MQ=255)
                                                            cTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGa                                       >  1:1693905/1‑100 (MQ=255)
                                                             tGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGa                                       <  1:3514085/100‑1 (MQ=255)
                                                             tGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGa                                       >  2:2715990/1‑100 (MQ=255)
                                                             tGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGCCTTATAATGCGGGTTTTGCAATGGa                                       <  1:1040963/99‑1 (MQ=255)
                                                             tGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGAt                                      >  2:2193666/1‑100 (MQ=255)
                                                             tGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGAt                                      >  1:2698101/1‑100 (MQ=255)
                                                             tGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTATGAAGTAACATGGCAATGCCTGGGGCTTTATAATGCGGGTTTTTCAATGGa                                       >  2:184030/1‑99 (MQ=255)
                                                               aaCATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATaa                                    >  2:4062349/1‑100 (MQ=255)
                                                                 cATTAAAGTAGGAACCGAGATTTTATATAACGATTTTTCACGTTGTGGATTGACGTGGCAATGCCTGTTGACTTATAATGCGGGTTTTGCAATGGATaa                                    >  2:1857716/1‑99 (MQ=255)
                                                                 cATTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAAcc                                  <  1:340367/100‑1 (MQ=255)
                                                                  aTTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCa                                 <  2:830173/100‑1 (MQ=255)
                                                                  aTTAAAATAGGAACCGAGATTTTATATAACCATTTTCACGATGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCa                                 >  2:2470341/1‑100 (MQ=255)
                                                                     aaaaTAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCaaaa                              <  2:713481/100‑1 (MQ=255)
                                                                         tAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAag                          <  1:3812230/100‑1 (MQ=255)
                                                                         tAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAag                          <  1:1467681/100‑1 (MQ=255)
                                                                         tAGGAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAag                          <  1:1761689/100‑1 (MQ=255)
                                                                           ggAACCGAGATTTTATATAACCATTTTTTTTTTTGTTGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAaga                         >  2:2066323/1‑100 (MQ=255)
                                                                           ggAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAaga                         <  1:1694986/100‑1 (MQ=255)
                                                                           ggAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAagag                        <  1:2267162/100‑1 (MQ=255)
                                                                           ggAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAagag                        >  2:2240966/1‑100 (MQ=255)
                                                                           ggAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAaga                         <  2:754338/99‑1 (MQ=255)
                                                                            gAACCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAgaga                       <  1:1944468/100‑1 (MQ=255)
                                                                              aCCGAGATTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACa                     <  1:784472/100‑1 (MQ=255)
                                                                               ccGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACa                     >  2:165001/1‑100 (MQ=255)
                                                                                 gagaTTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCa                  <  2:1721294/100‑1 (MQ=255)
                                                                                   gaTTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCaa                 >  1:594261/1‑100 (MQ=255)
                                                                                    aTTTTATATAACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAata               >  1:1362880/1‑100 (MQ=255)
                                                                                         atataACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGccc          <  1:2528436/100‑1 (MQ=255)
                                                                                         atataACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGcc           >  2:555405/1‑99 (MQ=255)
                                                                                          tataACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGATTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGcccc         >  2:2593042/1‑100 (MQ=255)
                                                                                            taACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCaa       >  2:3125130/1‑100 (MQ=255)
                                                                                            taACCATTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCaa       >  1:721585/1‑100 (MQ=255)
                                                                                                 aTTTTCCACCTCATGGAGTAACGAGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAata   <  1:3902057/100‑1 (MQ=255)
                                                                                                  tttttCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATAGATAAACAAAAAAAGAGACAGCAATATGCCCCAAAatat  <  2:3346103/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CCCCTGTTAAAAAAACGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCAAAATAT  >  NC_000913/2994188‑2994385

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