Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 231,868 T→G 6.9% A249A (GCT→GCG yafD → endo/exonuclease/phosphatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913231,8680TG 6.9% 769.7 / 10.2 286A249A (GCT→GCGyafDendo/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base T (121/138);  new base G (2/24);  total (124/162)
Fisher's exact test for biased strand distribution p-value = 5.28e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

AACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGA  >  NC_000913/231719‑232017
                                                                                                                                                     |                                                                                                                                                     
aaCGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCt                                                                                                                                                       >  2:208560/1‑150 (MQ=255)
  cgctTTATATCGCTTTGCTCGGGTATTTTCTCTGCGCGAGGTTCTTTTTACCGATGTTCTGCGCGGTGGTGGGTTTGGTCGCCCTCTCGTTTTTGTTTTCTACCGTTGTCTTTACGTCTGTGATGCTTCTGTACTGGTTACGCGCGCTTc                                                                                                                                                     <  1:344383/150‑1 (MQ=255)
  cgcgTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTc                                                                                                                                                     >  2:2008555/1‑150 (MQ=255)
   gcgTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCGTcc                                                                                                                                                    <  1:836087/150‑1 (MQ=255)
     gTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGa                                                                                                                                                  <  2:3580899/150‑1 (MQ=255)
     gTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGa                                                                                                                                                  <  1:3397681/150‑1 (MQ=255)
      ttATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGAt                                                                                                                                                 <  1:1056257/150‑1 (MQ=255)
      ttATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCAGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCGTCCGAt                                                                                                                                                 <  1:4249674/150‑1 (MQ=255)
        atatCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATca                                                                                                                                               <  1:3644520/150‑1 (MQ=255)
           tCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGTAGCTTCTGTACTGGTTTCGCGCGCTTCCGATAACaa                                                                                                                                            >  1:2853734/1‑150 (MQ=255)
           tCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACaa                                                                                                                                            >  1:1654785/1‑150 (MQ=255)
               tttGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTATGTAATGGTTACGCGCGCTTCCGATCACAAACCg                                                                                                                                        >  1:1920053/1‑150 (MQ=255)
                    gcgGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTAca                                                                                                                                   >  1:3816122/1‑149 (MQ=255)
                    gcgGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACt                                                                                                                                   >  2:3329979/1‑150 (MQ=255)
                    gcgGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACt                                                                                                                                   >  2:2902588/1‑150 (MQ=255)
                    gcgGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACt                                                                                                                                   >  1:283533/1‑150 (MQ=255)
                    gcgGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCGCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACt                                                                                                                                   >  1:56933/1‑150 (MQ=255)
                        gAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGtt                                                                                                                               <  2:4087228/150‑1 (MQ=255)
                        gAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGtt                                                                                                                               <  1:1294537/150‑1 (MQ=255)
                          aaTGACGCTGCGCCAGGTGCGATTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGAATTTGTAAACTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGa                                                                                                                             <  2:1281342/150‑1 (MQ=255)
                           aTGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGaa                                                                                                                            <  2:158302/150‑1 (MQ=255)
                             gTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGGTTGGTCGCCCGCTGGATTTTGTTTTCTACCGGGGTCTGAACGTCAGTGAAGCGTCTGTACTGGTTACGCGCGCTTCCGCTCACACTCCGCTCCACGTTGaagt                                                                                                                          >  1:3496449/1‑148 (MQ=255)
                              tCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGGTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTc                                                                                                                         >  2:2102200/1‑150 (MQ=255)
                                gCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTTAAttaat                                                                                                                       >  1:2582153/1‑147 (MQ=255)
                                gCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAg                                                                                                                       <  2:4252582/150‑1 (MQ=255)
                                 cTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGt                                                                                                                      <  2:4142035/150‑1 (MQ=255)
                                    cgcCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCGTCCGATCACAATCCGCTACTCGTTGAATTCAGTccc                                                                                                                   <  2:2202068/150‑1 (MQ=255)
                                     gcCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCg                                                                                                                  <  1:3100272/150‑1 (MQ=255)
                                      ccAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGCTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCgg                                                                                                                 <  1:1898772/150‑1 (MQ=255)
                                        aGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCa                                                                                                               >  2:3042847/1‑150 (MQ=255)
                                         ggTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCaa                                                                                                              <  2:4203837/150‑1 (MQ=255)
                                         ggTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCaa                                                                                                              <  2:276074/150‑1 (MQ=255)
                                              gTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTg                                                                                                         >  2:499366/1‑150 (MQ=255)
                                               ttttACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGa                                                                                                        <  2:1820781/150‑1 (MQ=255)
                                                attACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTACTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGat                                                                                                       <  2:1099182/149‑1 (MQ=255)
                                                     cGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGataaata                                                                                                  <  1:2542832/150‑1 (MQ=255)
                                                      gatgatCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGataaataa                                                                                                 <  2:2042625/150‑1 (MQ=255)
                                                      gatgatCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGataaataa                                                                                                 <  2:1166670/150‑1 (MQ=255)
                                                        tgatCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAgg                                                                                               >  2:3100168/1‑150 (MQ=255)
                                                        tgatCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAgg                                                                                               <  2:886072/150‑1 (MQ=255)
                                                        tgatCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCGTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAgg                                                                                               <  2:259995/150‑1 (MQ=255)
                                                          atCAGCGCCGTCGGGCGTTTTGTCGCCCCCTCGATTTTTTTTTATAACCTGGTCTGAACGTCAGTGAATCTTATGTAATGGTTACGCGCGCTTACGATAACAAAACCCGACACGTTGAATTAAGTCCCCGCAAGCCGGATAAAAAAGGTa                                                                                             >  1:2223500/1‑150 (MQ=255)
                                                          atCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTa                                                                                             >  2:3152011/1‑150 (MQ=255)
                                                               cgcCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCa                                                                                        <  2:3702039/150‑1 (MQ=255)
                                                                      gggCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGc                                                                                 <  1:1013042/150‑1 (MQ=255)
                                                                       ggCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGcc                                                                                <  2:3588919/150‑1 (MQ=255)
                                                                            ttGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAgg                                                                           <  2:3518612/150‑1 (MQ=255)
                                                                            ttGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAgg                                                                           <  1:2351138/150‑1 (MQ=255)
                                                                               gTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCa                                                                        <  1:3310280/150‑1 (MQ=255)
                                                                                tCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAg                                                                       >  2:2843534/1‑150 (MQ=255)
                                                                                tCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAg                                                                       >  2:2490763/1‑150 (MQ=255)
                                                                                 cGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGa                                                                      <  2:1633699/150‑1 (MQ=255)
                                                                                 cGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCGTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGa                                                                      <  2:2330281/150‑1 (MQ=255)
                                                                                    ccGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCCTGTGGCAGAcct                                                                   <  1:2915658/150‑2 (MQ=255)
                                                                                     cGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCaa                                                                  >  2:3213821/1‑150 (MQ=255)
                                                                                      gCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAAc                                                                 >  1:4063186/1‑150 (MQ=255)
                                                                                          ttttttttttTCTACGTTGGTCTTAACTTCTGTTATGCTTGTGTACTGGTTACGCGCGCTTCCGATCACAATCCGATACTGTTTTATTTCTGTCCGGGCAAGCTTTATAAATATGGTATGTCTGTTCTGCCACCGGGCAGACCAACGttt                                                             <  1:2043663/148‑1 (MQ=255)
                                                                                            ttttGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTgg                                                           <  2:1410240/150‑1 (MQ=255)
                                                                                            ttttGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTgg                                                           <  2:2516736/150‑1 (MQ=255)
                                                                                             tttGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGgc                                                          >  2:7750/1‑150 (MQ=255)
                                                                                              ttGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGgcg                                                         >  1:3595992/1‑150 (MQ=255)
                                                                                              ttGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGgcg                                                         >  2:3839666/1‑150 (MQ=255)
                                                                                                gTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCt                                                       <  2:3217001/150‑1 (MQ=255)
                                                                                                 ttttATACCGTGGTCTGACCGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTTAATTCTGTCCCGGCACGCCTGATAAATAAGGTATGTCTGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTg                                                      <  1:3953196/150‑1 (MQ=255)
                                                                                                      tACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCaa                                                 >  1:3269387/1‑150 (MQ=255)
                                                                                                      tACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCaa                                                 >  1:4176581/1‑150 (MQ=255)
                                                                                                       aCCGTGGTCTGAACGTCTGTGCAGCTTCTGTACTGGTTACGCGCGCTTCCGAACACAATCCGCTACTCGTTGAATTCTGTCCCGGCACGCCTTATAAATAAGGTATGTCTGGTCTGCCACTGGGCAGACCAACGTTTGGCGCTGCGCaaa                                                <  1:981220/150‑1 (MQ=255)
                                                                                                       aCCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCaaa                                                <  2:2976255/150‑1 (MQ=255)
                                                                                                        ccGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCaaaa                                               <  2:3757896/150‑1 (MQ=255)
                                                                                                          gTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTTGCGCAGCGCCAAACg                                             <  1:3543200/150‑1 (MQ=255)
                                                                                                             gTCTGAACGTCCGTGAAGCGTCTGGACTGGGTCCGCGCGCGGCCGAGCACAAGCCGCTACTCGTTGAAGTCCGTCCCGGCAAGCCTGATAAATACGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTTGCGCAGCGCCAAACGttg                                          <  2:1223126/150‑3 (MQ=255)
                                                                                                             gTCTGAACGTCAGTGAAGCTTCTGTTCTGGTTACGCGCGCGTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATACATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGa                                          <  2:2360807/150‑1 (MQ=255)
                                                                                                              tCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAg                                         <  2:1920040/150‑1 (MQ=255)
                                                                                                              tCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAg                                         <  2:1920153/150‑1 (MQ=255)
                                                                                                               cTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAgc                                        >  1:2491657/1‑150 (MQ=255)
                                                                                                               cTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAgc                                        >  1:1142371/1‑150 (MQ=255)
                                                                                                                tGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCGTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTTGCGCAGCGCCAAACGttggtc                                       <  2:3435053/150‑6 (MQ=255)
                                                                                                                   aCGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCgg                                    <  1:2973578/150‑1 (MQ=255)
                                                                                                                      tCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCggtgg                                 <  2:2441806/150‑1 (MQ=255)
                                                                                                                       cAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCggtggt                                >  1:2747150/1‑150 (MQ=255)
                                                                                                                        aGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGgtggtg                               <  1:96601/150‑1 (MQ=255)
                                                                                                                        aGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGAACACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGgtggtg                               >  1:1197902/1‑150 (MQ=255)
                                                                                                                          tGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGc                             <  1:725520/150‑1 (MQ=255)
                                                                                                                               cgtCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATg                        <  2:3279950/148‑1 (MQ=255)
                                                                                                                               cTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATg                        <  2:1080458/150‑1 (MQ=255)
                                                                                                                                ttCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGa                       >  1:2106033/1‑150 (MQ=255)
                                                                                                                                ttCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGa                       >  1:2457800/1‑150 (MQ=255)
                                                                                                                                ttCTGTACTGGTTACGCGCGCGTCCGATCACACTCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGa                       <  2:598123/150‑1 (MQ=255)
                                                                                                                                 tCTGTACTGGTTACGCGCGCGTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGAc                      <  2:267924/150‑1 (MQ=255)
                                                                                                                                     tACTGGTTACGCGCGCTTCCGATCCCACTCCGCTACTCGTTGAATTCAGGCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTAt                  <  2:1555800/150‑1 (MQ=255)
                                                                                                                                     tACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTAt                  >  1:1049559/1‑150 (MQ=255)
                                                                                                                                     tACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTAt                  >  1:2584917/1‑150 (MQ=255)
                                                                                                                                       cTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCt                >  1:1715966/1‑150 (MQ=255)
                                                                                                                                          gTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGcc             <  2:2407758/150‑1 (MQ=255)
                                                                                                                                          gTTACGCGCGCGTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGcc             <  2:273809/150‑1 (MQ=255)
                                                                                                                                           ttACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGccc            >  1:791490/1‑150 (MQ=255)
                                                                                                                                           ttACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGccc            >  2:958741/1‑150 (MQ=255)
                                                                                                                                             aCGCGCGCGTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCcac          <  2:412471/150‑1 (MQ=255)
                                                                                                                                               gcgcgcTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACAt        >  1:1853226/1‑150 (MQ=255)
                                                                                                                                                  cgcTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCt     <  1:2283547/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAAATAAGGTATGTCAGGTCTGCCACAGGGCAGACCAACGTTTGGCGCTGCGCAAAACGTGAGCGCGGTGGTGGCGTATGACTTATCTGCCCACATGCTGGA  >  NC_000913/231719‑232017

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: