Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 339,856 T→C 49.2% intergenic (+113/‑309) yahG → / → yahI DUF1116 family protein/carbamate kinase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913339,8560TC49.2% 11.5 / 311.2 226intergenic (+113/‑309)yahG/yahIDUF1116 family protein/carbamate kinase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (47/67);  new base C (50/61);  total (97/129)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.57e-01

Only 100 of 230 total aligned reads displayed.

GCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCG  >  NC_000913/339757‑339955
                                                                                                   |                                                                                                   
   GGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGGATTTGTCACCACTGCCGGATgcg                                                                                                  <  2:1197792‑M1/100‑1 (MQ=255)
    GATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATgcgg                                                                                                 >  2:1314003‑M1/1‑100 (MQ=255)
    GATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACgcgg                                                                                                 >  2:3935703‑M1/1‑100 (MQ=255)
    GATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACgcgg                                                                                                 >  1:3313116‑M1/1‑100 (MQ=255)
      TGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACgcggcg                                                                                               >  2:3301423‑M1/1‑100 (MQ=255)
        cggcgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTg                                                                                             <  1:2870284/100‑1 (MQ=38)
        cggcgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTg                                                                                             <  1:1409130/100‑1 (MQ=38)
        cggcgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTg                                                                                             <  1:2525880/100‑1 (MQ=255)
        cggcgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTg                                                                                             >  1:2817825/1‑100 (MQ=38)
         gtcgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGa                                                                                            >  1:2781634/3‑100 (MQ=255)
            gTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACg                                                                                         <  1:1547317/100‑1 (MQ=255)
             tGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGc                                                                                        <  1:93856/100‑1 (MQ=40)
             tGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGc                                                                                        <  1:4049021/100‑1 (MQ=255)
              gAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGcc                                                                                       <  1:460630/100‑1 (MQ=40)
              gAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGcc                                                                                       <  1:3215637/100‑1 (MQ=38)
              gAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGcc                                                                                       <  2:1621423/100‑1 (MQ=255)
              gAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGcc                                                                                       >  1:911307/1‑100 (MQ=255)
                aCGCCTTACCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGcc                                                                                       <  2:1598712/98‑1 (MQ=25)
                  gCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATACGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAAAGCCTTAt                                                                                   >  2:850741/1‑100 (MQ=255)
                  gCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCCTTTGTCACCATTGCCGGATGCGGCGTGCACGCCTTAt                                                                                   >  2:1567322/1‑100 (MQ=255)
                   ccTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGGCACCATTGCCGGACGCGGCGTGAACGCCTTATc                                                                                  >  2:1018926/1‑100 (MQ=255)
                    cttccgcGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATcc                                                                                 <  2:2574227/94‑1 (MQ=255)
                        tCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGcc                                                                             <  2:3872096/100‑1 (MQ=255)
                           ggCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGcc                                                                             <  2:3282518/97‑1 (MQ=25)
                             ccTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGccagatc                                                                        >  2:4036185/1‑95 (MQ=255)
                             ccTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGa                                                                        >  1:1517192/1‑100 (MQ=38)
                              cTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGaa                                                                       >  1:2574242/1‑100 (MQ=255)
                               tACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACTATTGCAGGATGCGGCGAGAACGAGTTATCCGGCCTACg                                                                         >  2:3591043/1‑97 (MQ=11)
                               tACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAAt                                                                      <  1:120586/100‑1 (MQ=255)
                               tACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTCCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAAt                                                                      <  2:2213884/100‑1 (MQ=255)
                                 cGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCTGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATgg                                                                    <  2:2201729/100‑1 (MQ=255)
                                  gAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGgc                                                                   >  2:1458069/1‑100 (MQ=18)
                                   aaTGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGgcg                                                                  >  2:2125182/1‑100 (MQ=40)
                                      ggCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGgcgc                                                                 <  2:2271896/98‑1 (MQ=38)
                                        cgcAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGa                                                             >  1:3623322/1‑100 (MQ=40)
                                        cgcAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGa                                                             >  1:727275/1‑100 (MQ=255)
                                         gcAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGaa                                                            <  1:1931360/100‑1 (MQ=40)
                                         gcAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGaa                                                            <  1:3559596/100‑1 (MQ=255)
                                         gcAACAATCTGTAGGCCGGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAg                                                              <  2:2748293/98‑1 (MQ=255)
                                            aGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCt                                                         >  1:2233642/1‑100 (MQ=255)
                                             gAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTg                                                        <  1:3187908/100‑1 (MQ=255)
                                               aTCTGTAGGCCTGGTAAGCGTAGCGCATCAGGCATTTGTCGCCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTa                                                      <  2:1706275/100‑1 (MQ=255)
                                               aTCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTa                                                      >  1:670464/1‑100 (MQ=255)
                                                 cTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAgg                                                    >  1:2673102/1‑100 (MQ=255)
                                                 cTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAgg                                                    >  1:3229568/1‑100 (MQ=255)
                                                   gTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGcc                                                  >  1:2390193/1‑100 (MQ=255)
                                                   gTAGGCATGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGcc                                                  <  1:1392085/100‑1 (MQ=255)
                                                    tAGGCCTGACAAGCGTAGCGCCGCAGGCCTTTCGCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGACCTACGAATGGCGCAAGAATCTGTAGGcc                                                  <  2:3249857/99‑1 (MQ=14)
                                                     aGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTg                                                <  1:3630279/100‑1 (MQ=255)
                                                        aaTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATa                                             >  1:3836392‑M1/3‑100 (MQ=255)
                                                        CCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATa                                             <  1:169865‑M1/100‑1 (MQ=255)
                                                          TGAGAAGCGTAGCGCATCAGGCATTTGTCACCAGGGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATa                                             >  1:3687243‑M1/1‑98 (MQ=255)
                                                            ATAAGCGTAGCGCATCAGGCATTTGTCATCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCg                                         >  1:3457441‑M1/1‑100 (MQ=255)
                                                              aggcctAGCGCTTCAGGCACTACGCACCATTGACCTACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTa                                       <  2:4001238/95‑1 (MQ=11)
                                                               aGCGTAGCGGATGAGGCCTTTGACATCCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGAt                                              <  2:3654668/92‑1 (MQ=14)
                                                               aGCGTAGCGCATTAGGAATTTGTCGCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAg                                           <  2:3976093/95‑1 (MQ=255)
                                                                gCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAgc                                     <  1:1106162/100‑1 (MQ=255)
                                                                gCGTAGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAgc                                     <  1:3209985/100‑1 (MQ=255)
                                                                 cGTAGCGCATCAGGCATTTGTTACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAgcg                                    <  2:647698/100‑1 (MQ=255)
                                                                 cGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAgcg                                    >  2:1084995/1‑100 (MQ=255)
                                                                  gTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCGGTAGGCCTGATAAGCGt                                        >  2:3124895/1‑95 (MQ=255)
                                                                   tAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCGAGAATCGGTAGGCCTGATAGGCGTAccgc                                   >  2:484814/1‑99 (MQ=11)
                                                                    aGCGCATCAGGCATTTGTCACCATTGCCGGATGTGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCAt                                 >  1:4011785/1‑100 (MQ=255)
                                                                    aGCGCATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAgcgc                                   >  1:3780252/1‑98 (MQ=255)
                                                                       gcATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAg                              <  2:932693/100‑1 (MQ=255)
                                                                        cATCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAgg                             <  1:3848513/100‑1 (MQ=255)
                                                                         aTCAGGCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGc                            >  1:1774295/1‑100 (MQ=255)
                                                                             ggCATTTGTCGCCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                          <  2:575782/98‑1 (MQ=255)
                                                                             ggCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                        >  1:1724807/1‑100 (MQ=255)
                                                                             ggCATTTGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                        >  1:1613922/1‑100 (MQ=255)
                                                                                aTTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGcc                     >  1:204091/1‑100 (MQ=255)
                                                                                 tttGTCGCCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGcca                    <  1:2961874/100‑1 (MQ=255)
                                                                                  ttGTCGTCATTGCCGGATGCGGCGTGAGCGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGc                      <  1:1309517/97‑1 (MQ=255)
                                                                                   tgtcgtgcatGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGccacc                  <  2:3007616/91‑1 (MQ=2)
                                                                                   tgaggtccgtGCCGGATGCGGCGTGAACGGCTGATCCGGCATACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGccacc                  <  2:1976315/91‑1 (MQ=11)
                                                                                   tGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGAAAAGCGTAGCGCATCAGGCATTTGccac                   >  1:3727239/1‑99 (MQ=255)
                                                                                   tGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGccacc                  >  1:3082111/1‑100 (MQ=255)
                                                                                   tGTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGccacc                  <  1:3046378/100‑1 (MQ=255)
                                                                                    gTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGccacca                 >  2:2386506/1‑100 (MQ=255)
                                                                                    gTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGccacca                 >  1:1280054/1‑100 (MQ=255)
                                                                                    gTCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGccacca                 >  1:648597/1‑100 (MQ=255)
                                                                                     tCACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCcaccac                <  1:3983062/100‑1 (MQ=255)
                                                                                      cACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACt               >  1:3956841/1‑100 (MQ=255)
                                                                                      cACCATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGGCTGATAAGCGTAGCACATCAGGCATTTGccacc                  >  2:62235/1‑97 (MQ=255)
                                                                                       aCCCTTGCCGGACGCGGCATGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGccac                   <  2:441812/95‑1 (MQ=255)
                                                                                       aCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTg              >  1:1926973/1‑100 (MQ=255)
                                                                                        ccATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGc             >  1:3355026/1‑100 (MQ=255)
                                                                                        ccATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGCCCTGATAAGCGTAGCGCATCAAGCATTTGACACCCCtcc             >  2:1821752/1‑98 (MQ=11)
                                                                                         cATTGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATACGCGTAGCGCATCAGGCATTTGCCACCACTGcc            >  2:1162089/1‑100 (MQ=255)
                                                                                           ttGCCGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCgg          >  1:2297611/1‑100 (MQ=255)
                                                                                            tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGa         <  1:2484418/100‑1 (MQ=38)
                                                                                              ccGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATg       <  1:78170/100‑1 (MQ=255)
                                                                                              ccGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATg       >  1:3905399/1‑100 (MQ=255)
                                                                                              ccGGACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATg       >  1:430656/1‑100 (MQ=255)
                                                                                                ggATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATgcg     <  1:3884028/100‑1 (MQ=38)
                                                                                                ggACGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATgcg     >  2:272615/1‑100 (MQ=255)
                                                                                                 gATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATgcgg    <  1:1056843/100‑1 (MQ=38)
                                                                                                   tGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATATGCGTAGCGCATCAGGCATTTGCCACCGCTGCCGGATgcgg    >  2:786019/1‑98 (MQ=255)
                                                                                                   tGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATgcggcg  <  2:3465128/100‑1 (MQ=38)
                                                                                                   tGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATgcggcg  <  1:1069136/100‑1 (MQ=38)
                                                                                                   |                                                                                                   
GCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCG  >  NC_000913/339757‑339955

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: