Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,613,561 C→T R27R (CGC→CGT)  yjjV → putative DNase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,613,5610CT97.7% 726.2 / 3.0 219R27R (CGC→CGT) yjjVputative DNase
Reads supporting (aligned to +/- strand):  ref base C (1/4);  new base T (118/94);  total (119/100)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

Only 100 of 226 total aligned reads displayed.

AACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTAT  >  NC_000913/4613464‑4613660
                                                                                                 |                                                                                                   
aaCAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACCCCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAAcgtgc                                                                                                   <  2:4688635/100‑1 (MQ=255)
 aCAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAAcgtgcg                                                                                                  <  2:2867621/100‑1 (MQ=255)
  cAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAAcgtgcgc                                                                                                 >  2:1153977/1‑99 (MQ=255)
  cAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCgg                                                                                                 >  2:4098501/1‑100 (MQ=255)
  cAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGGTTAGTGGCGATGAAGAGGCCCGCCTGCAACGTGCgg                                                                                                 <  2:5351401/100‑1 (MQ=255)
   aaTGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGc                                                                                                <  2:610729/100‑1 (MQ=255)
   aaTGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGc                                                                                                >  2:2779340/1‑100 (MQ=255)
    aTGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGca                                                                                               <  2:1368383/100‑1 (MQ=255)
    aTGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGca                                                                                               >  1:3018000/1‑100 (MQ=255)
     tGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGcac                                                                                              >  1:3147031/1‑100 (MQ=255)
      gAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGcaca                                                                                             <  2:2092706/100‑1 (MQ=255)
      gAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACGGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGcaca                                                                                             <  1:485666/100‑1 (MQ=255)
         gATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGc                                                                                          >  1:3043226/1‑100 (MQ=255)
            cTGGTTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCggg                                                                                       <  2:3611652/100‑1 (MQ=255)
            cTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCggg                                                                                       <  2:332514/100‑1 (MQ=255)
             tGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGc                                                                                      <  1:391813/100‑1 (MQ=255)
             tGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGc                                                                                      <  2:1504894/100‑1 (MQ=255)
              ggCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGTGGCACAAGCGGGCg                                                                                     <  2:2012701/100‑1 (MQ=255)
              ggCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCg                                                                                     >  2:747574/1‑100 (MQ=255)
              ggCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCg                                                                                     >  1:600963/1‑100 (MQ=255)
                cTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTa                                                                                   <  1:2032073/100‑1 (MQ=255)
                cTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTGTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTa                                                                                   <  2:2348204/100‑1 (MQ=255)
                  tGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAgg                                                                                 <  2:5308949/100‑1 (MQ=255)
                  tGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAgg                                                                                 <  1:4037969/100‑1 (MQ=255)
                    aTTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCa                                                                               >  1:1158052/1‑100 (MQ=255)
                       tGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGa                                                                            >  1:3348399/1‑100 (MQ=255)
                         ccGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATc                                                                          >  1:5244612/1‑100 (MQ=255)
                          cGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCa                                                                         >  2:3254977/1‑100 (MQ=255)
                          cGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCa                                                                         >  2:3970708/1‑100 (MQ=255)
                             tttATCTACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCACGATCATTg                                                                      <  2:3734132/100‑1 (MQ=255)
                             tttATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTg                                                                      >  1:785800/1‑100 (MQ=255)
                             tttATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTg                                                                      <  1:5541552/100‑1 (MQ=255)
                             tttATCGACACCCACTGCCATTTTGATTTCACGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTg                                                                      >  1:2063480/1‑100 (MQ=255)
                              ttATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGt                                                                     <  2:3895402/100‑1 (MQ=255)
                                aTCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTc                                                                   >  1:281768/1‑100 (MQ=255)
                                aTCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTc                                                                   >  1:4628760/1‑100 (MQ=255)
                                aTCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTc                                                                   >  1:1926059/1‑100 (MQ=255)
                                 tCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTcc                                                                  >  2:3483028/1‑100 (MQ=255)
                                 tCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTcc                                                                  >  1:5454689/1‑100 (MQ=255)
                                  cGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATCAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCAGTGTTCCg                                                                 <  2:4534697/100‑1 (MQ=255)
                                   gACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCgg                                                                <  2:1581692/100‑1 (MQ=255)
                                   gACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCgg                                                                >  2:1935773/1‑100 (MQ=255)
                                   gACACCCACTACCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCgg                                                                >  1:4109191/1‑100 (MQ=255)
                                    acacCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGc                                                               <  1:1229903/100‑1 (MQ=255)
                                     cacCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCa                                                              <  1:4541104/100‑1 (MQ=255)
                                     cacCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCa                                                              <  1:1537138/100‑1 (MQ=255)
                                       cccACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAAc                                                            >  2:5013927/1‑100 (MQ=255)
                                       cccACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAAc                                                            >  2:4466036/1‑100 (MQ=255)
                                       cccACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAAc                                                            >  1:2361331/1‑100 (MQ=255)
                                        ccACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGAGCATTGTTCCGGCAACa                                                           >  1:1626493/1‑100 (MQ=255)
                                         cACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACag                                                          >  1:3423759/1‑100 (MQ=255)
                                         cACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATCAAGTGGCCAGCCTGCAAGGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACag                                                          >  2:5164813/1‑100 (MQ=255)
                                             gCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGc                                                      >  1:3196878/1‑100 (MQ=255)
                                             gCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGc                                                      <  1:1861589/100‑1 (MQ=255)
                                             gCCATTTTGATTTCCCGCCGTTTAATGGCAATGAAGAGGCCAGCCTGCAACATGCGGCACAAACGGGCGTAGGCAAGAGCATTGTTCCGGCAACAGAGGc                                                      >  1:5190352/1‑100 (MQ=255)
                                                aTTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGa                                                   >  1:5048875/1‑100 (MQ=255)
                                                 ttttGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGaa                                                  >  1:2216146/1‑100 (MQ=255)
                                                 ttttGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGaa                                                  >  2:1036977/1‑100 (MQ=255)
                                                   ttGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGaaaa                                                >  1:2240390/1‑100 (MQ=255)
                                                    tGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGTAACAGAGGCGGAAAAt                                               <  2:5322660/100‑1 (MQ=255)
                                                      aTTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAAttt                                             <  2:3085246/100‑1 (MQ=255)
                                                      aTTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAAttt                                             >  2:2514960/1‑100 (MQ=255)
                                                        ttCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTg                                           >  2:49188/1‑100 (MQ=255)
                                                        ttCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTg                                           <  1:233462/100‑1 (MQ=255)
                                                          cccGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCAGTGTTCCGGCAACAGAGGCGGAAAATTTTGcc                                         >  2:4833920/1‑100 (MQ=255)
                                                           ccgccgTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGccc                                        >  1:2183392/1‑100 (MQ=255)
                                                                gTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCgtgtg                                   <  2:3199477/100‑1 (MQ=255)
                                                                    aGTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTgg                               >  1:4062800/1‑100 (MQ=255)
                                                                     gTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGc                              >  1:2634946/1‑100 (MQ=255)
                                                                     gTGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGc                              <  2:4256472/100‑1 (MQ=255)
                                                                      tGGCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTATGCCCGTGTGTTGGCa                             <  2:5327169/100‑1 (MQ=255)
                                                                       ggCGATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCAt                            >  2:555804/1‑100 (MQ=255)
                                                                          gATGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAg                         >  1:3599238/1‑100 (MQ=255)
                                                                           aTGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGc                        <  2:2999792/100‑1 (MQ=255)
                                                                           aTGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGc                        <  2:1983020/100‑1 (MQ=255)
                                                                            tGAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCg                       >  1:4970891/1‑100 (MQ=255)
                                                                             gAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCgg                      <  2:4805351/100‑1 (MQ=255)
                                                                             gAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCgg                      >  2:2543356/1‑100 (MQ=255)
                                                                             gAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCgg                      <  1:3988611/100‑1 (MQ=255)
                                                                             gAAGAGGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCgg                      <  1:303147/100‑1 (MQ=255)
                                                                               agagGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGaa                    <  1:1040757/100‑1 (MQ=255)
                                                                                gagGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGaaa                   >  1:5383130/1‑100 (MQ=255)
                                                                                gagGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGaaa                   >  1:3689206/1‑100 (MQ=255)
                                                                                 agGCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGaaaa                  <  1:4720870/100‑1 (MQ=255)
                                                                                  ggCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAAt                 <  2:4897744/100‑1 (MQ=255)
                                                                                  ggCCAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAAt                 <  2:80975/100‑1 (MQ=255)
                                                                                    ccAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTa               >  1:3839702/1‑100 (MQ=255)
                                                                                    ccAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTa               <  1:291674/100‑1 (MQ=255)
                                                                                     cAGCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTAt              <  1:717910/100‑1 (MQ=255)
                                                                                       gCCTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCa            <  2:4649798/100‑1 (MQ=255)
                                                                                        ccTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCaa           <  2:1891836/100‑1 (MQ=255)
                                                                                        ccTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCaa           >  2:1256614/1‑100 (MQ=255)
                                                                                         cTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAAc          <  1:4385708/100‑1 (MQ=255)
                                                                                         cTGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAAc          <  1:2104992/100‑1 (MQ=255)
                                                                                          tGCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCGTTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAAcc         >  1:1149314/1‑100 (MQ=255)
                                                                                           gCAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCg        >  2:5232599/1‑100 (MQ=255)
                                                                                            cAACGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGc       <  1:3515599/100‑1 (MQ=255)
                                                                                             aaCGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCCGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCt      <  1:2744559/100‑1 (MQ=255)
                                                                                              aCGTGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTg     >  1:2648431/1‑100 (MQ=255)
                                                                                              aCGGGCGGCACAAGCGGGCGTAGGCAAGATCAGTGTTAAGGCAACAGAGGCGGAAAATTTTGCACGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTg     <  1:344814/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
AACAATGAGGATCTGGCTTGATTTGCCGTTTTATCGACACCCACTGCCATTTTGATTTCCCGCCGTTTAGTGGCGATGAAGAGGCCAGCCTGCAACGCGCGGCACAAGCGGGCGTAGGCAAGATCATTGTTCCGGCAACAGAGGCGGAAAATTTTGCCCGTGTGTTGGCATTAGCGGAAAATTATCAACCGCTGTAT  >  NC_000913/4613464‑4613660

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