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)
        aaaaaaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGtt                                                                                           <  2:2339384/100‑1 (MQ=255)
         aaaaaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTtg                                                                                          >  1:456897/1‑100 (MQ=255)
         aaaaaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTtg                                                                                          >  2:1640435/1‑100 (MQ=255)
          aaaaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTtg                                                                                          >  2:3179496/1‑100 (MQ=255)
          aaaaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTtgt                                                                                         <  1:2044493/100‑1 (MQ=255)
           aaaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTtgtg                                                                                        >  2:378708/1‑100 (MQ=255)
           aaaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTtgtg                                                                                        <  2:1907694/100‑1 (MQ=255)
            aaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTgg                                                                                       <  2:2335097/100‑1 (MQ=255)
            aaaCGTGGAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTgg                                                                                       <  1:838996/100‑1 (MQ=255)
                  ggAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTaa                                                                                  >  2:130007/1‑100 (MQ=255)
                  ggAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAAc                                                                                 >  2:2384192/1‑100 (MQ=255)
                   gAATAATAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACg                                                                                <  1:3240713/100‑1 (MQ=255)
                    aataatAAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGt                                                                               >  1:3660829/1‑100 (MQ=255)
                     ataataAAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTg                                                                              >  1:4172106/1‑100 (MQ=255)
                     ataataAAAGATAAACGAGATGTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTa                                                                                   <  2:2146317/95‑1 (MQ=255)
                       ttccgatcctGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTgg                                                                             <  1:3066949/91‑1 (MQ=255)
                       aataaAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGc                                                                            >  1:3150448/1‑100 (MQ=255)
                       aataaAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGc                                                                            <  1:37511/100‑1 (MQ=255)
                       aataaAAGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGc                                                                            >  2:797509/1‑100 (MQ=255)
                          aaaaGATGAACGAGATTTGTTAGATAACCTTTGTCTGACCATTAAAATAAGAACCGAAATTCTATTTAACTA‑TTTTAAAGTTGTGGAGGAAcaacacaat                                                                         >  2:4043052/1‑92 (MQ=255)
                           aaaGATGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATg                                                                        >  1:3350284/1‑100 (MQ=255)
                                tGAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTggagatcggaa                                                                    >  2:3066931/1‑91 (MQ=255)
                                 gAACGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTgg                                                                   <  1:2129360/100‑1 (MQ=255)
                                   aCGAGATTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGa                                                                <  1:2969083/100‑1 (MQ=255)
                                      agaTTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACtt                                                             <  2:2752201/100‑1 (MQ=255)
                                      agaTTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACt                                                              >  2:3446145/1‑99 (MQ=255)
                                      agaTTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCNGGGGACt                                                              <  2:250558/99‑1 (MQ=255)
                                       gaTTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTta                                                            <  1:1887227/100‑1 (MQ=255)
                                       gaTTTGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTta                                                            >  1:690960/1‑100 (MQ=255)
                                         tttGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACTGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTTACAATGCCTGGGGACTtata                                                          >  2:2871618/1‑100 (MQ=255)
                                          ttGTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATaa                                                         <  2:3660802/100‑1 (MQ=255)
                                            gTTAGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATg                                                       >  1:2140417/1‑100 (MQ=255)
                                               aGATAACCCTTGTCTGAACATTAAAATAGGAACCGAGAGTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCg                                                     >  1:224168/1‑99 (MQ=255)
                                                gATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCggg                                                   <  1:1165437/100‑1 (MQ=255)
                                                gATAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCggg                                                   >  1:3191333/1‑100 (MQ=255)
                                                 aTAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCggg                                                   <  2:2117779/100‑1 (MQ=255)
                                                  tAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGc                                                      >  2:2565182/1‑95 (MQ=255)
                                                  tAACCCTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGt                                                  >  1:567695/1‑99 (MQ=255)
                                                     cccTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGtttt                                               <  2:3501305/99‑1 (MQ=255)
                                                      ccTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGc                                             <  1:4126914/100‑1 (MQ=255)
                                                      ccTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGc                                             <  1:3523083/100‑1 (MQ=255)
                                                       cTTGTCTGAACATTAACATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTg                                              <  1:3610347/98‑1 (MQ=255)
                                                       cTTGTCTGAACATTAAAATAGGATCCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCa                                            <  1:1740064/100‑1 (MQ=255)
                                                       cTTGTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGc                                             <  2:3215918/100‑1 (MQ=255)
                                                          gTCTGAACATTAAAATAGGAACCGAGATTTTATATATCCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATg                                         <  1:2477179/100‑1 (MQ=255)
                                                          gTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATg                                         >  1:1466225/1‑100 (MQ=255)
                                                          gTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATg                                         <  1:2469000/100‑1 (MQ=255)
                                                          gTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATg                                         >  1:4072763/1‑100 (MQ=255)
                                                          gTCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATg                                         <  2:1022442/100‑1 (MQ=255)
                                                           tCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATgg                                        <  2:1356751/100‑1 (MQ=255)
                                                           tCTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATgg                                        <  1:331129/100‑1 (MQ=255)
                                                            cTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTTGGGACTTATAATGCGGGATTTGAAACGGa                                       >  2:690342/1‑100 (MQ=255)
                                                            cTGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGa                                       >  1:1693905/1‑100 (MQ=255)
                                                             tGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGa                                       <  1:3514085/100‑1 (MQ=255)
                                                             tGAACATTAAAATAGGAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGa                                       >  2:2715990/1‑100 (MQ=255)
                                                             tGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGCCTTATAATGCGGGTTTTGCAATGGa                                       <  1:1040963/99‑1 (MQ=255)
                                                             tGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGAt                                      >  2:2193666/1‑100 (MQ=255)
                                                             tGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGAt                                      >  1:2698101/1‑100 (MQ=255)
                                                             tGAACATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTATGAAGTAACATGGCAATGCCTGGGGCTTTATAATGCGGGTTTTTCAATGGa                                       >  2:184030/1‑99 (MQ=255)
                                                               aaCATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATaa                                    >  2:4062349/1‑100 (MQ=255)
                                                                 cATTAAAGTAGGAACCGAGATTTTATATAACGATTTTTCACGTTGTGGATTGACGTGGCAATGCCTGTTGACTTATAATGCGGGTTTTGCAATGGATaa                                    >  2:1857716/1‑99 (MQ=255)
                                                                 cATTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAAcc                                  <  1:340367/100‑1 (MQ=255)
                                                                  aTTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCa                                 <  2:830173/100‑1 (MQ=255)
                                                                  aTTAAAATAGGAACCGAGATTTTATATAACCA‑TTTTCACGATGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCa                                 >  2:2470341/1‑100 (MQ=255)
                                                                     aaaaTAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCaaaa                              <  2:713481/100‑1 (MQ=255)
                                                                         tAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAag                          <  1:3812230/100‑1 (MQ=255)
                                                                         tAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAag                          <  1:1467681/100‑1 (MQ=255)
                                                                         tAGGAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAag                          <  1:1761689/100‑1 (MQ=255)
                                                                           ggAACCGAGATTTTATATAACCATTTTTTTTTTTGTTGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAaga                         >  2:2066323/1‑100 (MQ=255)
                                                                           ggAACCGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAaga                         <  1:1694986/100‑1 (MQ=255)
                                                                           ggAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAagag                        <  1:2267162/100‑1 (MQ=255)
                                                                           ggAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAagag                        >  2:2240966/1‑100 (MQ=255)
                                                                           ggAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAaga                         <  2:754338/99‑1 (MQ=255)
                                                                            gAACCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAgaga                       <  1:1944468/100‑1 (MQ=255)
                                                                              aCCGAGATTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACa                     <  1:784472/100‑1 (MQ=255)
                                                                               ccGAGATTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACa                     >  2:165001/1‑100 (MQ=255)
                                                                                 gagaTTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCa                  <  2:1721294/100‑1 (MQ=255)
                                                                                   gaTTTTATATAACCATTTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCaa                 >  1:594261/1‑100 (MQ=255)
                                                                                    aTTTTATATAACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAata               >  1:1362880/1‑100 (MQ=255)
                                                                                         atataACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGccc          <  1:2528436/100‑1 (MQ=255)
                                                                                         atataACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGcc           >  2:555405/1‑99 (MQ=255)
                                                                                          tataACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGATTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGcccc         >  2:2593042/1‑100 (MQ=255)
                                                                                            taACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCaa       >  2:3125130/1‑100 (MQ=255)
                                                                                            taACCA‑TTTTCACGTTGTGGAGTAACGTGGCAATGCCTGGGGACTTATAATGCGGGTTTTGCAATGGATAACCAAAAAAAGAGACAGCAATATGCCCCaa       >  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: ‑