Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,422,895 C→G 15.7% A17G (GCC→GGC)  ulaE → L‑xylulose 5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,422,8950CG15.7% 448.7 / 71.6 166A17G (GCC→GGC) ulaEL‑xylulose 5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (79/61);  new base G (12/14);  total (91/75)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

Only 100 of 177 total aligned reads displayed.

CACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCT  >  NC_000913/4422796‑4422994
                                                                                                   |                                                                                                   
    tCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCggg                                                                                                 <  2:2475267/100‑1 (MQ=255)
    tCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCggg                                                                                                 <  2:3010546/100‑1 (MQ=255)
     cAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCgggg                                                                                                <  2:2783401/100‑1 (MQ=255)
       gTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGa                                                                                               >  1:3362255/1‑99 (MQ=255)
       gTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAg                                                                                              >  2:2111722/1‑100 (MQ=255)
       gTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAg                                                                                              >  2:47443/1‑100 (MQ=255)
       gTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAg                                                                                              >  2:3428783/1‑100 (MQ=255)
         tCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTg                                                                                            <  2:1942108/100‑1 (MQ=255)
          cAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGGCGGGGAgcg                                                                                            <  1:3758015/99‑3 (MQ=255)
          cAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTgc                                                                                           >  2:3827574/1‑100 (MQ=255)
          cAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTgc                                                                                           >  2:817913/1‑100 (MQ=255)
          cAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCACTCCCCGGCGGGGAgcgc                                                                                           <  1:3300808/100‑1 (MQ=255)
           aaaCGTTCCATCGCTGAACTGTGGGGCTAATGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTgct                                                                                          >  1:2805315/1‑100 (MQ=255)
           aaaCGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTgct                                                                                          >  1:3577704/1‑100 (MQ=255)
           aaaCGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCACTCCCCGGCGGGGAGCgct                                                                                          <  2:3710599/100‑1 (MQ=255)
            agcGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCACTCCCCGGCGGGGAGCgctt                                                                                         <  1:2428709/98‑2 (MQ=255)
            aaCGTTCGATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTCGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTgctg                                                                                         <  1:2291568/100‑1 (MQ=255)
             aCGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTgctgg                                                                                        <  1:64881/100‑1 (MQ=255)
              cGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTgctggc                                                                                       >  2:1308921/1‑100 (MQ=255)
               gTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTgctggct                                                                                      >  1:1696125/1‑100 (MQ=255)
               gTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCACTCCCCGGCGGGGAGCgcttttt                                                                                      <  2:2867993/100‑5 (MQ=38)
               gTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCACTCCCCGGCGGGGAGCgcttttt                                                                                      <  2:3105906/100‑5 (MQ=38)
               gTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCACTCCCCGGCGGGGAGCgcttttt                                                                                      <  2:2660770/100‑5 (MQ=38)
                 tCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAGACAAATCCCGCTTGGCATCTATGAAAAAGCACTCCCCGGCGGGGAGCgctttttca                                                                                    <  1:2732141/100‑7 (MQ=37)
                 tCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCACTCCCCGGCGGGGAGCgctttttca                                                                                    <  2:1836609/100‑7 (MQ=38)
                   cATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGaa                                                                                  >  1:955172/1‑100 (MQ=255)
                       gCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATTAAAAAGCGCTCCCAGCCGGGGAGTGCTGGCTGGAACGcc                                                                              >  1:689570/1‑100 (MQ=255)
                       gCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGcc                                                                              >  2:3048771/1‑100 (MQ=255)
                       gCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGcc                                                                              >  1:3124773/1‑100 (MQ=255)
                       gCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGcc                                                                              <  2:3946511/100‑1 (MQ=255)
                         tGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCTAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGTACGCCTg                                                                            <  1:2246051/100‑1 (MQ=255)
                         tGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTg                                                                            <  2:1736004/100‑1 (MQ=255)
                         tGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTg                                                                            <  2:3910962/100‑1 (MQ=255)
                          gAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGc                                                                           <  2:3674973/100‑1 (MQ=255)
                            aCTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCaa                                                                         >  2:1135039/1‑100 (MQ=255)
                            aCTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCaa                                                                         >  1:3867906/1‑100 (MQ=255)
                              tgtgGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACt                                                                       >  2:2154371/1‑100 (MQ=255)
                              tgtgGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACt                                                                       >  2:33364/1‑100 (MQ=255)
                                tgGGGCTAAGGAGCGGATATGTTGTACAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTgg                                                                     >  1:1527752/1‑100 (MQ=255)
                                 ggggCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGc                                                                    <  2:2244032/100‑1 (MQ=255)
                                  gggCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCa                                                                   >  2:858879/1‑100 (MQ=255)
                                   ggCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCaa                                                                  >  2:2019053/1‑100 (MQ=255)
                                      tAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCaaaaa                                                               >  2:3272585/1‑100 (MQ=255)
                                      tAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCaaaaa                                                               <  1:4070163/100‑1 (MQ=255)
                                       aaGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAAc                                                              >  2:1422679/1‑100 (MQ=255)
                                         ggcgCGGATATGTTGTGCAAGCAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAACACg                                                             >  1:2613875/4‑99 (MQ=255)
                                           aGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTa                                                          <  2:258609/100‑1 (MQ=255)
                                           aGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTa                                                          <  2:3391594/100‑1 (MQ=255)
                                            gCGGATATGTTGTCCAAACAAATCCCGCGTGGCATCTATGAAAAAGCGCTCCCCGCCGGAGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTa                                                          <  1:1244840/99‑1 (MQ=255)
                                              gatcttgttgtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGc                                                       <  1:1639981/95‑1 (MQ=255)
                                                 tatGTTGTCCAAACAAATCCCGCTTGGCATGTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGTTACGCCTTCTACTGGCAAAAACGTTAGGCtta                                                    <  1:2467910/100‑2 (MQ=255)
                                                 tatGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTc                                                    >  2:443469/1‑100 (MQ=255)
                                                   tgttgtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGcagatc                                                  >  2:1639985/1‑95 (MQ=255)
                                                   tgttgtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGa                                                  >  1:2782085/1‑100 (MQ=255)
                                                    gttgtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGAt                                                 >  2:1865261/1‑100 (MQ=255)
                                                       gtCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGAtttt                                              >  2:69374/1‑100 (MQ=255)
                                                         ccAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGt                                            >  1:2956949/1‑100 (MQ=255)
                                                         ccAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGt                                            >  1:837981/1‑100 (MQ=255)
                                                         ccAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCAGCCGGGGAGTGCTGGCTGGAACGCCTGGAACTGGCAAAGACGTTAGGCTTCGAtttt                                              >  1:1797257/1‑98 (MQ=255)
                                                          caaacaaaTCCCGCTTGTCATCTCTGAAAAGGCGCCCCCCGCCGGGGTGTGCTGGCTGGATCGCCTGCAGCTTGCAATAACGTTAGGCCTCGATTTTg                                             <  1:3810239/98‑1 (MQ=255)
                                                          caaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTc                                           >  2:2404953/1‑100 (MQ=255)
                                                          caaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTc                                           >  2:2516356/1‑100 (MQ=255)
                                                          caaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTc                                           <  1:2063967/100‑1 (MQ=255)
                                                           aaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCg                                          <  1:3127320/100‑1 (MQ=255)
                                                           aaacaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCg                                          >  1:1887421/1‑100 (MQ=255)
                                                             acaaaTCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGaa                                        >  1:4126490/1‑100 (MQ=255)
                                                                  tCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTc                                   <  2:4260621/100‑1 (MQ=255)
                                                                  tCCCGCTTGGCATCTATGAAAAAGCGCTACCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGACTTTGTCGAAATGTc                                   >  2:1165227/1‑100 (MQ=255)
                                                                    ccGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCgg                                 >  2:4050554/1‑100 (MQ=255)
                                                                         tGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGAt                            >  2:3143718/1‑100 (MQ=255)
                                                                         tGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGAt                            >  1:960558/1‑100 (MQ=255)
                                                                          ggCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATg                           <  2:1191842/100‑1 (MQ=255)
                                                                          ggCATCCATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATg                           <  1:494091/100‑1 (MQ=255)
                                                                            cATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGaa                         <  2:341610/100‑1 (MQ=255)
                                                                              cgttccagccaGCACTCCCCGGCGGGGAGCGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAAc                       >  1:1167245/11‑100 (MQ=37)
                                                                               cTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCTTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACt                      <  1:2482714/100‑1 (MQ=255)
                                                                               cTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACt                      <  2:2297528/100‑1 (MQ=255)
                                                                               cTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACt                      >  2:3930374/1‑100 (MQ=255)
                                                                                 aTGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGa                    <  2:511682/100‑1 (MQ=255)
                                                                                 aTGAAAAAGCGCGCCCCGCCGGGGAGTGCTGGCTGGAACACCTGCAACTGGCAAAGACTTTAGGCTTCGATTTTGTCGGAATGTCGGTAGATGAAAc                       <  1:1030112/97‑1 (MQ=255)
                                                                                  tGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGAc                   >  1:1211724/1‑100 (MQ=255)
                                                                                   gAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACg                  >  1:3263800/1‑100 (MQ=255)
                                                                                    naaaaGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGa                 >  1:1187634/2‑100 (MQ=255)
                                                                                      ccagcacTCCCCGGCGGGGAGCGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATc               >  2:707556/7‑100 (MQ=38)
                                                                                      ccagcacTCCCCGGCGGGGAGCGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATc               >  2:3988634/7‑100 (MQ=38)
                                                                                       caGCACTCCCCGGCGGGGAGCGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCg              >  1:343028/2‑100 (MQ=39)
                                                                                       caGCACTCCCCGGCGGGGAGCGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCg              >  2:3406167/2‑100 (MQ=39)
                                                                                         gcgcTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGcc            <  2:942240/100‑1 (MQ=255)
                                                                                         gcgcTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGcc            <  2:355760/100‑1 (MQ=255)
                                                                                          cgcTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCt           >  1:1032777/1‑100 (MQ=255)
                                                                                           gcTCCCCGCTGGGGAGTGTTGGCTGGTACGCCTGCAACTGGCAAAAACGGTTGGCGTTGTTTTTGTCGTAATGTCGGTAGGTGAAACTGACGATCGCCTg          >  1:2767329/1‑100 (MQ=255)
                                                                                             tCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTATGCTTCGAGTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGACTGTc        >  1:1463335/1‑100 (MQ=255)
                                                                                             tCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTc        >  2:2962398/1‑100 (MQ=255)
                                                                                              ccccGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTcg       >  2:4182936/1‑100 (MQ=255)
                                                                                              ccccGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTcg       >  1:2380576/1‑100 (MQ=255)
                                                                                               cccGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTcgc      >  2:2086353/1‑100 (MQ=255)
                                                                                               cccGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTcgc      >  2:1874306/1‑100 (MQ=255)
                                                                                                ccgccgGGGAGTGCTGGCTGGAACGCGTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTcgcg     >  1:3417355/1‑100 (MQ=255)
                                                                                                ccgccgGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATccc             >  1:349754/1‑90 (MQ=255)
                                                                                                   ccgGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCt  >  2:2462273/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGGCAAAAACGTTAGGCTTCGATTTTGTCGAAATGTCGGTAGATGAAACTGACGATCGCCTGTCGCGCCT  >  NC_000913/4422796‑4422994

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: