Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 339,886 A→G 5.6% intergenic (+143/‑279) 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,8860AG 5.6% 787.2 / 14.6 249intergenic (+143/‑279)yahG/yahIDUF1116 family protein/carbamate kinase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (92/136);  new base G (12/2);  total (111/138)
Fisher's exact test for biased strand distribution p-value = 1.27e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

GTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCAT  >  NC_000913/339738‑340034
                                                                                                                                                    |                                                                                                                                                    
gTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAAt                                                                                                                                                     <  1:4209154/150‑1 (MQ=255)
 tCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATg                                                                                                                                                    <  2:4340150/150‑1 (MQ=255)
   tGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGgc                                                                                                                                                  >  1:3737305/1‑150 (MQ=255)
   tGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGgc                                                                                                                                                  <  1:3428844/150‑1 (MQ=255)
       cTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCaa                                                                                                                                              <  2:2761575/150‑1 (MQ=255)
       cTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCaa                                                                                                                                              <  2:53294/150‑1 (MQ=255)
         aCATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGACTGGCGCgaga                                                                                                                                            >  1:4214217/1‑150 (MQ=255)
         aCATGTTCATGCCGGAAGCGGCGTGAACGCCTTATCCGGCCTACGAAAGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGAATCAGGCATTTGTCACCATTGACGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCAAGa                                                                                                                                            >  2:1374489/1‑150 (MQ=14)
          cATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGaa                                                                                                                                           <  1:3549620/150‑1 (MQ=255)
             gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTCCGAGTGGCGCGAGTAtgg                                                                                                                                        >  2:1098113/1‑148 (MQ=255)
             gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCt                                                                                                                                        <  1:3636124/150‑1 (MQ=255)
              ttCATGCCGGATGCGGCGTGAACGCCTTCTCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTg                                                                                                                                       <  2:3340065/150‑1 (MQ=255)
              gtCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTg                                                                                                                                       <  1:2338517/149‑1 (MQ=255)
                cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTa                                                                                                                                     <  2:2684418/150‑1 (MQ=255)
                cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTa                                                                                                                                     <  2:4050238/150‑1 (MQ=255)
                cATGCCGGAGGCGGCGAGAACGCCTGCTCCGGGCTCCGAATGGCGCAAGACTCTGTAGGCCGGATCAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTCCGAATGGCGCAAGAATCTGTa                                                                                                                                     <  2:2309183/150‑1 (MQ=14)
                 aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAg                                                                                                                                    <  1:514774/150‑1 (MQ=255)
                  tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTCgg                                                                                                                                   >  1:1336749/1‑150 (MQ=255)
                     cGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCt                                                                                                                                <  1:3633123/150‑1 (MQ=255)
                     cGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCt                                                                                                                                <  1:1397603/150‑1 (MQ=255)
                     cGGATGCGGCGTGAACGCCTGATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCt                                                                                                                                <  2:3849170/150‑1 (MQ=255)
                      ggATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATCCGCGTAGCGCAGCTGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGACTGGCGCGAGACTCTGTTGGCCtt                                                                                                                               >  1:1466756/1‑149 (MQ=11)
                       gttGCGGCGTGAACGCCGTCTCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGa                                                                                                                              <  2:60363/148‑1 (MQ=255)
                       gATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATCGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGa                                                                                                                              <  2:3767721/150‑1 (MQ=255)
                            ggcgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGc                                                                                                                         >  1:3971457/1‑150 (MQ=255)
                               gTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTa                                                                                                                      <  2:1381254/150‑1 (MQ=255)
                                 gAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAgc                                                                                                                    <  1:311129/150‑1 (MQ=255)
                                 gAACGCCTTAGCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAgc                                                                                                                    <  1:2251190/150‑1 (MQ=255)
                                  aaCGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAgcg                                                                                                                   >  1:1187701/1‑150 (MQ=255)
                                  aaCGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAgcg                                                                                                                   <  1:2259901/150‑1 (MQ=255)
                                   aCGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCGAGAATCGGTAGGCCTGATGAGCGTAgcgc                                                                                                                  >  2:1190982/1‑150 (MQ=255)
                                    cGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCa                                                                                                                 >  1:2674071/1‑150 (MQ=255)
                                    cGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCa                                                                                                                 >  2:3439098/1‑150 (MQ=255)
                                      ccTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATc                                                                                                               <  1:2106751/150‑1 (MQ=255)
                                       cTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCa                                                                                                              <  1:2586381/150‑1 (MQ=255)
                                        ttATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCTAGAATCTTTGGGCTTGTTAGGCGAGGCGCaacgg                                                                                                             >  1:1999938/1‑146 (MQ=11)
                                        ttATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                                                             >  1:712940/1‑150 (MQ=255)
                                         tATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                                            >  1:3479067/1‑150 (MQ=255)
                                         tATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                                                            <  2:3589767/150‑1 (MQ=255)
                                           tCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCACGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                                                          >  2:397399/1‑150 (MQ=255)
                                             cGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                                                                                                        <  1:1516923/150‑1 (MQ=255)
                                              ggCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCAAGAATCGGTAGGCCGGTTAAGCGTAGCGCATCGGCCAttt                                                                                                       >  1:2649924/1‑150 (MQ=11)
                                              cgCCAACGAATGGCGCAACAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTGACACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                                                                                                       <  2:1872776/149‑1 (MQ=255)
                                                ccTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGc                                                                                                     >  1:3153426/1‑150 (MQ=255)
                                                  tACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGcca                                                                                                   <  2:899196/150‑1 (MQ=255)
                                                         ggCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGc                                                                                            <  2:1530556/150‑1 (MQ=255)
                                                           cgcAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGACACCAGAGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGACTGGCGCAAGAATCTGTGGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCg                                                                                          >  1:499635/1‑150 (MQ=14)
                                                             cAAGAATCTGTCGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGTATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGa                                                                                        <  1:2585627/150‑1 (MQ=255)
                                                             cAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGa                                                                                        >  2:3123873/1‑150 (MQ=255)
                                                             cAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGa                                                                                        <  1:3112972/150‑1 (MQ=255)
                                                                gAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATgc                                                                                     <  1:1590517/150‑1 (MQ=255)
                                                                    cTGTTGGCCTGAAAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATgcggcg                                                                                 <  2:1051060/150‑1 (MQ=255)
                                                                    cTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATgcggcg                                                                                 <  1:3325038/150‑1 (MQ=255)
                                                                    cTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATgcggcg                                                                                 <  2:3036113/150‑1 (MQ=255)
                                                                     tGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGt                                                                                >  1:2695783/1‑150 (MQ=255)
                                                                      gTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTg                                                                               <  2:334884/150‑1 (MQ=255)
                                                                        tggCCTGATTAGCGTAGCGCATCAGGCATTTGTCCCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGa                                                                             <  2:255520/149‑1 (MQ=255)
                                                                         ggCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGAc                                                                            <  2:4152719/150‑1 (MQ=255)
                                                                         ggCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGAc                                                                            <  2:1711715/150‑1 (MQ=255)
                                                                         ggCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGAc                                                                            <  1:2767314/150‑1 (MQ=255)
                                                                          gCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACg                                                                           >  2:3154285/1‑150 (MQ=255)
                                                                           ccTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGc                                                                          <  1:4296237/150‑1 (MQ=255)
                                                                              gATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCtt                                                                       >  1:2990224/1‑150 (MQ=255)
                                                                                   gCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCg                                                                  <  1:925258/150‑1 (MQ=255)
                                                                                   gCGGAGCGCATCAGGCATATGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCg                                                                  <  1:3795556/150‑1 (MQ=255)
                                                                                     gTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGc                                                                <  2:2278445/150‑1 (MQ=255)
                                                                                        gcgcATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCtc                                                             >  1:362360/1‑149 (MQ=255)
                                                                                        gcgcATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTa                                                             >  1:103953/1‑150 (MQ=255)
                                                                                        gcgcATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACAGCCGGATGCGGCGTGGACGCCTTATCCGGCCTa                                                             >  1:1906246/1‑150 (MQ=255)
                                                                                         cgcATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTAc                                                            >  2:1823550/1‑150 (MQ=255)
                                                                                           cATCCGGCATTGGGCCCCCTGGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGCATGGCGCAAGAATCTGTAGGCCTGCGAAGCGTAGCGCATCAGGCATTTGCCACCCCTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGa                                                          <  2:3939673/150‑1 (MQ=14)
                                                                                           cATCAGGCAATTGTCACCATTGCCGGATGCGGCGTGAACGCCTGATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGa                                                          <  1:4290013/150‑1 (MQ=255)
                                                                                            aTCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGTAAGCGGCGAGCACGCCTTATCCCGCCTACGac                                                         >  2:2114489/1‑149 (MQ=14)
                                                                                             tCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGt                                                        >  1:293321/1‑150 (MQ=255)
                                                                                               aGGCATTTGTCACCCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTgg                                                      <  1:1771441/150‑1 (MQ=255)
                                                                                               aGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTgg                                                      >  1:619515/1‑150 (MQ=255)
                                                                                               aGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCGACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTGCGAGTgg                                                      >  1:4118631/1‑150 (MQ=255)
                                                                                                  cacgtGTCACCCTTGCCGGATGCGGCGTGACCGCCTTCTCCGGCCTACGAAGGGCGCACGAATCTGTAGGCCTGCTAAGCGTAGCGCATCAGGCATTTGCCCCCCCTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGgcgc                                                   <  2:2180348/146‑1 (MQ=14)
                                                                                                        tCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAAt                                             >  1:2099283/1‑150 (MQ=255)
                                                                                                           ccATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTg                                          <  2:4198872/150‑1 (MQ=255)
                                                                                                           ccATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTg                                          <  2:3654818/150‑1 (MQ=255)
                                                                                                             aTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCTAGAATATGTAGCCCTGAAAAGCGTACCGAATCAGGCAATAGCCACCACTGCCTGATGCGGAGTGGACGACTTAACCGGCCTACGAGAGGCGCGAGAATGTGTa                                        >  2:4124241/1‑150 (MQ=14)
                                                                                                               tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCt                                                                        <  1:385278/116‑1 (MQ=39)
                                                                                                               tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAgg                                      <  2:781985/150‑1 (MQ=255)
                                                                                                                 ccGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGcc                                    <  1:2123465/150‑1 (MQ=255)
                                                                                                                  cGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCt                                   <  2:2519417/150‑1 (MQ=255)
                                                                                                                    gATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGa                                 >  2:2216172/1‑150 (MQ=255)
                                                                                                                       gcggcgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATaa                              <  1:1823550/150‑1 (MQ=255)
                                                                                                                         ggcgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGc                            <  2:3890979/150‑1 (MQ=255)
                                                                                                                         ggcgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGc                            >  1:1574185/1‑150 (MQ=255)
                                                                                                                          gcgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCg                           >  1:1849924/1‑150 (MQ=255)
                                                                                                                          gcgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATATGTAGGCCTGATAAGCg                           >  2:1639069/1‑150 (MQ=255)
                                                                                                                           cgTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGt                          >  1:3670027/1‑150 (MQ=255)
                                                                                                                               aaCGCCTTCTCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTAGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAgcg                      <  2:3000893/150‑1 (MQ=255)
                                                                                                                                aCGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAgcgc                     >  2:3240505/1‑150 (MQ=255)
                                                                                                                                    cTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCa                 >  2:1587162/1‑150 (MQ=255)
                                                                                                                                    cTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCCGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATcg                 >  2:3630117/1‑149 (MQ=255)
                                                                                                                                     ttATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATGTGTTGGCATGATATGCTTTGCGAAtatg                >  1:2490965/1‑147 (MQ=14)
                                                                                                                                               tACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCa      >  2:4294549/1‑150 (MQ=255)
                                                                                                                                               tACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCa      <  2:1956602/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTCTGACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGACGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCAT  >  NC_000913/339738‑340034

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: