Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,186,374 T→C R1149R (CGT→CGC rpoC → RNA polymerase beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,186,3740TC81.0% 544.2 / 132.3 252R1149R (CGT→CGCrpoCRNA polymerase beta prime subunit
Reads supporting (aligned to +/- strand):  ref base T (24/24);  new base C (89/115);  total (113/139)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-01

Only 100 of 270 total aligned reads displayed.

GACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGG  >  NC_000913/4186282‑4186465
                                                                                            |                                                                                           
    cccTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGa                                                                                         >  2:3178471/1‑93 (MQ=255)
      cTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGaaa                                                                                       >  1:4668457/1‑93 (MQ=255)
         gcgcGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAagag                                                                                    <  1:307978/93‑1 (MQ=255)
         gcgcGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAagag                                                                                    <  2:1175616/93‑1 (MQ=255)
         gcgcGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAagag                                                                                    >  1:1420688/1‑93 (MQ=255)
           gcgTATTCCGCAGTAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGcc                                                                                  <  1:6625539/93‑1 (MQ=255)
             gTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCgg                                                                                <  1:1872702/93‑1 (MQ=255)
             gTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCgg                                                                                >  2:3385169/1‑93 (MQ=255)
                   cGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGACGCACGTCGCCCGAAAGAGCCGGCAATcc                                                                          >  2:291589/1‑93 (MQ=255)
                   cGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATcc                                                                          >  2:4991264/1‑93 (MQ=255)
                   cGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATcc                                                                          >  2:4937820/1‑93 (MQ=255)
                   cGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATcc                                                                          <  1:2477188/93‑1 (MQ=255)
                   cGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATcc                                                                          <  2:2005026/93‑1 (MQ=255)
                    gcagcaGTCCGGCGGTCCCAAGGACATCACCGGTCGTCTGCCGCGCGTTGCGGCCCTGTTCGAAGCACGTCGCACGAAAGAGCCGGCAATCCt                                                                         <  1:3329144/93‑1 (MQ=255)
                    gCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCt                                                                         <  1:1373430/93‑1 (MQ=255)
                     cAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTg                                                                        <  1:5491156/93‑1 (MQ=255)
                     cAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTg                                                                        >  1:2769646/1‑93 (MQ=255)
                     cAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTg                                                                        <  1:3014658/93‑1 (MQ=255)
                     cAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTg                                                                        >  1:3282927/1‑93 (MQ=255)
                      aGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTgg                                                                       <  2:5694900/93‑1 (MQ=255)
                       ggAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGc                                                                      >  2:252189/1‑93 (MQ=255)
                        gAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCt                                                                     >  1:2193665/1‑93 (MQ=255)
                            ccGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGaaa                                                                 <  2:1648551/93‑1 (MQ=255)
                            ccGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGaaa                                                                 >  1:5670590/1‑93 (MQ=255)
                            ccGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGaaa                                                                 >  1:4152524/1‑93 (MQ=255)
                             cggcggTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAAt                                                                >  2:1464857/1‑93 (MQ=255)
                              ggcggTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATc                                                               >  2:642837/1‑93 (MQ=255)
                               gcggTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCa                                                              <  2:5105779/93‑1 (MQ=255)
                                 ggTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGc                                                            >  1:6054127/1‑93 (MQ=255)
                                  gTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCg                                                           >  1:651926/1‑93 (MQ=255)
                                   tACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCgg                                                          <  2:800843/93‑1 (MQ=255)
                                   gaCCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCgg                                                          <  1:3159836/92‑1 (MQ=255)
                                    aCCAAGGACATCACCGGTGGTCTGCCGCGAGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGt                                                         >  2:4861336/1‑93 (MQ=255)
                                     ccAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTa                                                        >  1:4985016/1‑93 (MQ=255)
                                     ccAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTa                                                        >  2:3694463/1‑93 (MQ=255)
                                      cAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTAt                                                       <  2:3868725/93‑1 (MQ=255)
                                      cAAGGACATCACCGGTCGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTAt                                                       <  1:1502885/93‑1 (MQ=255)
                                       aaGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATc                                                      >  2:5441210/1‑93 (MQ=255)
                                        aGGACATCACCGGTGGTCTGCCGCGCGGTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCg                                                     >  2:4069138/1‑93 (MQ=255)
                                         ggACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGt                                                    <  2:4817285/93‑1 (MQ=255)
                                           ccATCACCGGTGGTCGGCCGCGCGTTGCGGACCTGTCCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGttt                                                  <  1:3838192/92‑1 (MQ=255)
                                            cATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTc                                                 <  2:5857503/93‑1 (MQ=255)
                                             aTCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTcc                                                >  2:6007701/1‑93 (MQ=255)
                                              tCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCCGGCTGAAATCAGCGGTATCGTTTCCt                                               >  1:5646740/1‑93 (MQ=255)
                                               cACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCtt                                              <  2:2692142/93‑1 (MQ=255)
                                               cACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCtt                                              <  1:5588996/93‑1 (MQ=255)
                                               cACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCtt                                              >  2:349323/1‑93 (MQ=255)
                                               cACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCtt                                              <  2:6148605/93‑1 (MQ=255)
                                                cccGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTc                                             <  2:1311168/92‑1 (MQ=255)
                                                aCCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTc                                             <  1:4311319/93‑1 (MQ=255)
                                                aCCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTc                                             <  1:516670/93‑1 (MQ=255)
                                                 ccGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCg                                            >  1:5733542/1‑93 (MQ=255)
                                                    gtggtCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAGAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTa                                         >  2:2269437/1‑93 (MQ=255)
                                                    gtggtCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCggga                                         >  2:1686209/1‑91 (MQ=255)
                                                     tggcCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTaa                                        <  2:2561491/93‑1 (MQ=255)
                                                      ggtCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTaaa                                       <  1:2809973/93‑1 (MQ=255)
                                                          tGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGaaa                                   >  2:4200565/1‑93 (MQ=255)
                                                          tGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGaaa                                   <  1:5096393/93‑1 (MQ=255)
                                                           gCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAAc                                  >  1:2720143/1‑93 (MQ=255)
                                                           gCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAAc                                  <  2:4720251/93‑1 (MQ=255)
                                                            ccGCGCGTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAAcc                                 >  2:6506760/1‑93 (MQ=255)
                                                              gcgcgTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCaa                               >  1:6203824/1‑93 (MQ=255)
                                                              gcgcgTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCaa                               >  1:4163151/1‑93 (MQ=255)
                                                              gcgcgTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCaa                               >  1:3814641/1‑93 (MQ=255)
                                                              gcgcgTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCaa                               >  2:5727311/1‑93 (MQ=255)
                                                              gcgcgTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCaa                               <  2:2419762/93‑1 (MQ=255)
                                                              gcgcgTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCaa                               <  1:4012633/93‑1 (MQ=255)
                                                               cgcgTTGCGGACCTGTTCGAAGCACGTCTCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCaaa                              <  1:758284/93‑1 (MQ=255)
                                                                  gTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGt                           <  2:3069349/93‑1 (MQ=255)
                                                                  gTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGt                           <  1:3784738/93‑1 (MQ=255)
                                                                  gTTGCGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTACTTCTGTAAAGAAACCAAAGGt                           >  1:3535364/1‑93 (MQ=255)
                                                                    tGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTaa                         >  2:5256876/1‑93 (MQ=255)
                                                                      cGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAAc                       <  1:903587/93‑1 (MQ=255)
                                                                      cGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAAc                       >  1:3940293/1‑93 (MQ=255)
                                                                      cGGACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAAc                       <  2:4570029/93‑1 (MQ=255)
                                                                       ggACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAAcg                      <  1:529657/93‑1 (MQ=255)
                                                                       ggACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAAcg                      <  1:6163430/93‑1 (MQ=255)
                                                                       ggACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAAcg                      <  2:5270554/93‑1 (MQ=255)
                                                                        gACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAAcgt                     >  2:4279451/1‑93 (MQ=255)
                                                                        gACCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAAcgt                     <  1:1625770/93‑1 (MQ=255)
                                                                         gccTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAAcgtc                    <  2:1777988/92‑1 (MQ=255)
                                                                         aCCTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAAcgtc                    <  1:1069544/93‑1 (MQ=255)
                                                                          ccTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAAcgtcg                   >  1:2901345/1‑93 (MQ=255)
                                                                          ccTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAAcgtcg                   >  1:2055113/1‑93 (MQ=255)
                                                                           cTGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGAATCGTTTCCTTCCGTAAAGAAACCAAAGGAAAAcgtcgt                  >  1:4908144/1‑93 (MQ=255)
                                                                            tGTTCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACgtcgtc                 >  1:4207323/1‑93 (MQ=255)
                                                                              ttCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTg               >  2:3317271/1‑93 (MQ=255)
                                                                               tCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTgg              <  1:5814360/93‑1 (MQ=255)
                                                                               tCGAAGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTgg              <  1:1019743/93‑1 (MQ=255)
                                                                                cGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGt             <  2:2168728/93‑1 (MQ=255)
                                                                                  aaGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTa           <  2:2974903/93‑1 (MQ=255)
                                                                                  aaGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTa           >  2:6110677/1‑93 (MQ=255)
                                                                                   aGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTAt          >  1:52026/1‑93 (MQ=255)
                                                                                   aGCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTAt          <  2:5340387/93‑1 (MQ=255)
                                                                                    gCACGTCGCCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATc         <  1:6424960/93‑1 (MQ=255)
                                                                                     cACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCa        >  1:1274460/1‑93 (MQ=255)
                                                                                        gtcgccCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCAccc     <  2:4406473/93‑1 (MQ=255)
                                                                                         tcgtcCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCAcccc    >  1:2618909/1‑93 (MQ=255)
                                                                                         tcgccCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCAcccc    <  2:5623069/93‑1 (MQ=255)
                                                                                           gtcCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCgg  >  1:1562101/1‑93 (MQ=255)
                                                                                            |                                                                                           
GACACCCTGGCGCGTATTCCGCAGGAATCCGGCGGTACCAAGGACATCACCGGTGGTCTGCCGCGCGTTGCGGACCTGTTCGAAGCACGTCGTCCGAAAGAGCCGGCAATCCTGGCTGAAATCAGCGGTATCGTTTCCTTCGGTAAAGAAACCAAAGGTAAACGTCGTCTGGTTATCACCCCGG  >  NC_000913/4186282‑4186465

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