Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,995,4570TG41.9% 109.3 / 12.8 63G232G (GGT→GGGdapFdiaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base T (10/26);  new base G (26/0);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 1.30e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGA  >  NC_000913/3995309‑3995599
                                                                                                                                                    |                                                                                                                                                
aGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTTTGCAAGGgtt                                                                                                                                                >  2:633107/1‑151 (MQ=255)
   gTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGGTGCAGAAGGGAGTCAGCAAGGGTTgt                                                                                                                                              >  2:402485/1‑149 (MQ=255)
     tttCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTGCTgg                                                                                                                                           >  2:1197533/1‑151 (MQ=255)
     tttCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTg                                                                                                                                            <  2:922941/150‑1 (MQ=255)
      ttCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTgg                                                                                                                                           <  2:877252/150‑1 (MQ=255)
       tCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTCAGCAAGGGTTTCTTGcc                                                                                                                                         >  1:449755/1‑151 (MQ=255)
         cGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCga                                                                                                                                       <  2:1223672/151‑1 (MQ=255)
           gAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGGAGGAAGGAATCAGCAAGGGTTGCTgg                                                                                                                                           >  1:998551/1‑145 (MQ=255)
                cgcCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGTAGGGATTCAGCAAGGTTTGCTGGCCgaagaa                                                                                                                                   >  2:560184/1‑148 (MQ=255)
                     atatCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGGTTCAGCAAGGGTTGCTGGCCgaagaa                                                                                                                                   >  2:1041358/1‑143 (MQ=255)
                       atCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTgg                                                                                                                         <  1:140608/151‑1 (MQ=255)
                         cGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCGGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGaa                                                                                                                       <  1:819872/151‑1 (MQ=255)
                            ttttATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTTCTGGCCGAAGAAGGACGCGTGGAACTc                                                                                                                    >  2:67608/1‑151 (MQ=255)
                               tATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGGCGTAGGGATTCAGCAAGGGTTTGTGGCCGAAGAAGGACGCGTGGGACTcccc                                                                                                                 >  1:1042497/1‑151 (MQ=255)
                                aTGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCcg                                                                                                                <  2:38008/151‑1 (MQ=255)
                                 tGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGGTGGAGGAGGGGTTTAACAAGGGTTTCTGGCCCGAGAAGGACGCGCGGGACTTCCcgg                                                                                                               >  2:676746/1‑151 (MQ=255)
                                  gCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGTAGGGGGGGGTCTGCAAGGGTTGGTTGGCGGAGAAATACGCGGGGGAATcccccg                                                                                                               >  1:549819/1‑148 (MQ=255)
                                  gCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCcgg                                                                                                               <  2:464725/150‑1 (MQ=255)
                                   cAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGGAGGGAATAAGCCAGGGTTGCTGGCCGAAAGAGGACACGGGGAACTcccc                                                                                                                 >  1:43023/1‑147 (MQ=255)
                                   cAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCATGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTGCTGGCCGAAGAAAGACGCGTGGAACTcccc                                                                                                                 >  1:781533/1‑147 (MQ=255)
                                    aaGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAAATCCCccgcgg                                                                                                            >  2:258166/1‑151 (MQ=255)
                                        ggTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGgtcg                                                                                                         <  2:171726/150‑1 (MQ=255)
                                        ggTTAAGCGCGAGCAGATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGgtcg                                                                                                         <  1:218153/150‑1 (MQ=255)
                                         gTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTCAGCAAGGGTTGCTGGCCGAAGAAGTAAGCGGGGGACTCCCCCGCGggcg                                                                                                         >  2:644306/1‑149 (MQ=255)
                                          ttAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAAGGTTGCTGGCCGAAGAAGGACGCGGGGGACTCCCCGGGGGGCGTCt                                                                                                      >  2:822643/1‑151 (MQ=255)
                                           tAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCt                                                                                                      <  1:505036/150‑1 (MQ=255)
                                            aaGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCtt                                                                                                     <  2:1016636/150‑1 (MQ=255)
                                            aaGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTg                                                                                                    <  2:1161861/151‑1 (MQ=255)
                                            aaGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTTAGCAAGGGTTTGTGGGCGGAGAAGGACGCGGGGGAATCCCCGGCGGGCGTCtt                                                                                                     >  2:177512/1‑150 (MQ=255)
                                              gcgcGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGGTTTAGCAAGGGTTTCTGGCCGGAGAAAAACGCGGGGGACTCCCCGGCGGGCGTCt                                                                                                      >  2:724283/1‑147 (MQ=255)
                                               cgcgAGGATATTCGTTTACGCGTTTGTTGGCGTGGGGGAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCTGTCGTCTTGat                                                                                                  <  1:903102/150‑1 (MQ=255)
                                                   aGCATATTCGTTTACGCGTTTATGAGCGTGGGGCCGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGCTGCCGGCGAGAGTCATCACGGAATGCTGGCCGGAGACGCAAGCGCGGTAGTTCCCGGAGGCAGCCTTGATATTGc                                                                                             >  1:686396/1‑151 (MQ=255)
                                                       tatTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAACAAGGGTTGCTGGCCCAAGAAGTACGCGGGGAGCTCCCCCGCGGGCGTCCTGatattg                                                                                              >  1:681754/1‑144 (MQ=255)
                                                                aCGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTGCTGGCCGGAAAAGGACGCGGGGAACTCCCCGGCGGGCGGCTTGATATTGCCTGGGAAAggcc                                                                                 >  2:808423/1‑148 (MQ=255)
                                                                 cgcgTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGTGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCgg                                                                              <  1:1172968/151‑1 (MQ=255)
                                                                     tttATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAAGGTTGGTCGGCGGAGAAGTAAGCGTGGAACTCCCCCGCGGGCGGCTTTATATCGCCCGGGAAAGGCCggg                                                                             >  2:445678/1‑148 (MQ=255)
                                                                        aTGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCCGTAGGGGTTCCGCAAGGGTTGCTGGCCGAAGAAATAAGCGTGGAAATCCCCGGCGGGCGTCTTGATATCGCCTGGGAAGGGCCGGGGGAccc                                                                      >  1:600838/1‑151 (MQ=255)
                                                                        aTGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCAcc                                                                       <  2:1225410/150‑1 (MQ=255)
                                                                            gCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTGCTGGCCGGAAAAGTAAGCGTGGGACTCCCCGGCGGGCGGCTTTGTATCGCCCGGGAAAGGCCGGGGCAccccgt                                                                   >  2:584396/1‑147 (MQ=255)
                                                                             cGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGTATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAGCTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCCCCCGtt                                                                   <  1:1158322/149‑1 (MQ=255)
                                                                                  ggCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATg                                                            <  1:364671/151‑1 (MQ=255)
                                                                                      ggAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTg                                                        <  2:743534/151‑1 (MQ=255)
                                                                                       gagaAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGGAGCAAGGGGTCAGCCAGGGGTGCTGGCCGGAGACGGAAGCGGGGGAACCCCCGGCGGGCGTCCTGATATCACCTGCAAAGGTGCGGGGCCCCCGTtatattttaatg                                                        >  1:89736/1‑143 (MQ=255)
                                                                                        agaAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGGAGGGATTCAGCAAGGGTTTCTGGCCGAAGAAGGACGCGGGGAAATCCCCGGCGGGCGGCTTGGTATCGCCCGGGAAGGGCCGGGGGACCCGGtatat                                                               >  1:927395/1‑142 (MQ=255)
                                                                                         gaAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGcc                                                     <  1:1197178/151‑1 (MQ=255)
                                                                                         gaAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGcc                                                     <  1:812942/151‑1 (MQ=255)
                                                                                              ccAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCcggcg                                                <  1:787256/151‑1 (MQ=255)
                                                                                                  gCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAAGGTTTCTGGGCGAAGAAATACGCGTGGAACTCCCCGGCGGGCGTCTTGATATTGCCTGGGAAGGTTCGGGGCACCCCTTATATATGACTTGCCCCGCGGGAc                                            >  2:774912/1‑151 (MQ=255)
                                                                                                    cTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGGTTTCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGGCCGGGGCACCCGTTATATATTACTGGCCCGGCGGGACa                                           >  2:788424/1‑150 (MQ=255)
                                                                                                       cGGCAGCGGCGCGTGTGCGGCGGGTGCAGTAGGGGTTCAGCAAAGGTTTCGCGCCCAGGAAGGACGCGCGGGGATCCCCCGCGGGCGGCTTGGTATCGCCCGGGCAGGGCCGGGGGACCCCGTTAGTATGACTGGgccccgggggaa                                           >  1:179108/1‑138 (MQ=255)
                                                                                                            gcggcgCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTacga                                  <  2:628048/151‑1 (MQ=255)
                                                                                                             cggcgCGTGTGCGGCGGTTGCAGTAGGGGTTCAGCAAGGGTTGCTGGCCGGAGAAGGACGCGTGGAACTCCCCGGCGGGCGTCTTGGTATTGGCCGGGAAGGGCCGGGGCACCCGCTGTATATGACTGGCCCGGGGGGGCATGTCTac                                    >  1:383678/1‑148 (MQ=255)
                                                                                                                 gcgTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGAt                             <  1:86598/151‑1 (MQ=255)
                                                                                                                         ggcggTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGGTTTATTCAt                     >  2:763036/1‑151 (MQ=255)
                                                                                                                           cggTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCt                   >  1:165601/1‑151 (MQ=255)
                                                                                                                            ggTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGGCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCt                   >  2:780243/1‑150 (MQ=255)
                                                                                                                             gTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTAt                 <  1:232460/151‑1 (MQ=255)
                                                                                                                                  aGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGc            >  2:793873/1‑151 (MQ=255)
                                                                                                                                     aGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAAc         >  2:769982/1‑151 (MQ=255)
                                                                                                                                     aGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAAc         <  1:67596/151‑1 (MQ=255)
                                                                                                                                       ggATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAAcc        >  2:52194/1‑150 (MQ=255)
                                                                                                                                        gATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAg      <  2:512471/151‑1 (MQ=255)
                                                                                                                                              gCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGga  >  2:406990/1‑149 (MQ=255)
                                                                                                                                               cAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGCCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAg      >  1:887655/1‑144 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGCGGCAGCGGCGCGTGTGCGGCGGTTGCAGTAGGGATTCAGCAAGGTTTGCTGGCCGAAGAAGTACGCGTGGAACTCCCCGGCGGTCGTCTTGATATCGCCTGGAAAGGTCCGGGTCACCCGTTATATATGACTGGCCCGGCGGTACATGTCTACGACGGATTTATTCATCTATGAAGCAACCAGGGGA  >  NC_000913/3995309‑3995599

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