Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,253,616 A→G 5.0% I200T (ATT→ACT)  yhaJ ← LysR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,253,6160AG 5.0% 253.3 / 2.9 79I200T (ATT→ACT) yhaJLysR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (36/39);  new base G (2/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

ACCATAGGATACGGCATCGTCGCCACGCCAAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGCTTGCGCGAGTTGATCTCCGACGAGGAACGAAAATGCATATCCGGCGC  >  NC_000913/3253469‑3253764
                                                                                                                                                   |                                                                                                                                                    
aCCATAGGATACGGCATCGTCGCCACGCCAAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCCCCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAAt                                                                                                                                                     >  8:358383/1‑149 (MQ=255)
 ccATAGGATACGGCATCGTCGCCACGCCAAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAAtt                                                                                                                                                    <  3:71461/149‑1 (MQ=255)
   aTAGGATACGGCATCGTCGCCACGCCAAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTc                                                                                                                                                   >  7:304320/1‑148 (MQ=255)
    tAGGATACGGCATCGTCGCCACGCCAAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCa                                                                                                                                                 <  7:369007/149‑1 (MQ=255)
    tAGGATACGGCATCGTCGCCACGCCAAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCa                                                                                                                                                 <  1:234329/149‑1 (MQ=255)
             gCATCGTCGCCACGCCAAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTc                                                                                                                                        <  8:13085/149‑1 (MQ=255)
             gCATCGTCGCCACGCCAAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTc                                                                                                                                        <  4:154844/149‑1 (MQ=255)
              cATCGTCGCCACGCCAAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTca                                                                                                                                       >  2:194915/1‑149 (MQ=255)
               aTCGTCGCCACGCCAAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTcac                                                                                                                                      >  6:208399/1‑149 (MQ=255)
                     gccacgccaAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCTCGATATTTCACACGCGt                                                                                                                                <  8:304320/149‑1 (MQ=255)
                        acgccaAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCAc                                                                                                                             >  2:197855/1‑149 (MQ=255)
                            caAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACtt                                                                                                                           >  7:315066/1‑147 (MQ=255)
                                   cGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACg                                                                                                                  <  8:315066/149‑1 (MQ=255)
                                   cGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACg                                                                                                                  >  4:69445/1‑149 (MQ=255)
                                     ccAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGt                                                                                                                >  2:19640/1‑149 (MQ=255)
                                       aGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTc                                                                                                              <  7:68895/149‑1 (MQ=255)
                                       aGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTc                                                                                                              <  4:216355/149‑1 (MQ=255)
                                           aTGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGc                                                                                                          <  1:194915/149‑1 (MQ=255)
                                             gCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTc                                                                                                        >  5:206788/1‑149 (MQ=255)
                                              ccTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTcc                                                                                                       <  1:325295/149‑1 (MQ=255)
                                               cTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCt                                                                                                      >  7:385597/1‑149 (MQ=255)
                                               cTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCt                                                                                                      >  5:183439/1‑149 (MQ=255)
                                                  tcGTTTATCTTCATTCGTGCTGACCGTTTAGCGGGGGTGTTTGTCCAGCATCTTTACGGTCAACACGGGGCTCTCTCTAGCGGTATCCGCCACCGCAATTACACGATATTTCACTCGCGTCTCTTTAGATAACGGTTCGGGCTcgttat                                                                                                   <  6:168319/148‑6 (MQ=255)
                                                  aCGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCtgat                                                                                                   <  6:206788/149‑1 (MQ=255)
                                                  aCGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCtga                                                                                                    <  4:257151/148‑1 (MQ=255)
                                                         tCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCgg                                                                                            >  4:88707/1‑149 (MQ=255)
                                                          cTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCggg                                                                                           >  8:39080/1‑149 (MQ=255)
                                                           ttCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTc                                                                                                                               <  2:240719/112‑1 (MQ=255)
                                                           ttCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTc                                                                                                                               >  1:240719/1‑112 (MQ=255)
                                                             cAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGt                                                                                          <  8:159113/147‑1 (MQ=255)
                                                             cAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGt                                                                                          <  3:69445/147‑1 (MQ=255)
                                                              aaTCGTGCTCACCGTTAAGCGCGGCTGTTGGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCCCACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGAt                                                                                       >  2:254099/1‑149 (MQ=255)
                                                                tCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGAt                                                                                       >  5:254161/1‑147 (MQ=255)
                                                                       cACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCgg                                                                              >  4:152356/1‑149 (MQ=255)
                                                                        aCCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGgcgc                                                                                <  3:354220/146‑1 (MQ=255)
                                                                        aCCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGgcgc                                                                                >  4:354220/1‑146 (MQ=255)
                                                                        aCCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGcgcg                                                                               <  1:19640/147‑1 (MQ=255)
                                                                        aCCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCgg                                                                              <  7:130263/148‑1 (MQ=255)
                                                                          cGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCaa                                                                           >  6:264963/1‑149 (MQ=255)
                                                                            ttAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACg                                                                         >  7:357581/1‑149 (MQ=255)
                                                                               aGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTa                                                                       <  3:88707/148‑1 (MQ=255)
                                                                               aGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAg                                                                      <  4:379195/149‑1 (MQ=255)
                                                                                gcgcGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGa                                                                     <  3:152356/149‑1 (MQ=255)
                                                                                 cgcgGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGAc                                                                    <  1:197855/149‑1 (MQ=255)
                                                                                 cgcgGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGAc                                                                    <  5:208399/149‑1 (MQ=255)
                                                                                    ggCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGtt                                                                 >  2:223548/1‑149 (MQ=255)
                                                                                     gCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGtt                                                                 >  7:12332/1‑148 (MQ=255)
                                                                                     gCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTc                                                                <  7:358383/149‑1 (MQ=255)
                                                                                        gTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAGTTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCAt                                                              >  3:33647/1‑148 (MQ=255)
                                                                                        gTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAGTTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCAt                                                              >  3:28505/1‑148 (MQ=255)
                                                                                            gTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAAc                                                          <  1:223548/148‑1 (MQ=255)
                                                                                             tCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAGTTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGt                                                        <  4:33647/149‑1 (MQ=255)
                                                                                             tCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAGTTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGt                                                        <  4:28505/149‑1 (MQ=255)
                                                                                               cagcagCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCTCGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGta                                                       <  6:183439/148‑1 (MQ=255)
                                                                                                     cTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACAcgcg                                                                                                                                 <  8:266220/68‑1 (MQ=255)
                                                                                                     cTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACAcgcg                                                                                                                                 >  7:266220/1‑68 (MQ=255)
                                                                                                       gTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTcc                                                                                                             <  4:244313/86‑1 (MQ=255)
                                                                                                       gTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTcc                                                                                                             >  3:244313/1‑86 (MQ=255)
                                                                                                        tACGGTCAGCACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCTACGTAGACGTTCTTTAAGGTATCGAGCTTgc                                             >  4:75507/1‑149 (MQ=255)
                                                                                                        tACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGCTTgc                                             >  7:402430/1‑149 (MQ=255)
                                                                                                        tACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAAACGTTCATTAACGTATAGAGCTTgc                                             >  4:272351/1‑149 (MQ=255)
                                                                                                        tACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTGCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTTTAGAGCTTgc                                             >  5:201063/1‑149 (MQ=255)
                                                                                                          cGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGCTTgcgc                                           >  2:167841/1‑149 (MQ=255)
                                                                                                           ggTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGCTTgcgc                                           >  2:221588/1‑148 (MQ=255)
                                                                                                             tCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGCTTGCGCGAg                                        <  8:385597/149‑1 (MQ=255)
                                                                                                             tCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGCTTGCGCGAg                                        <  5:264963/149‑1 (MQ=255)
                                                                                                             tCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAAACGTTCATTAACGTATAGAGCTTGCGCGAg                                        <  3:272351/149‑1 (MQ=255)
                                                                                                                 cacCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGc                                                 <  1:211548/136‑1 (MQ=255)
                                                                                                                 cacCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGc                                                 >  2:211548/1‑136 (MQ=255)
                                                                                                                 cacCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGCTTGCGCGAGTTGa                                    >  8:337247/1‑149 (MQ=255)
                                                                                                                        cgcTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGCTTGCGCGAGTTGATCTCcga                             <  6:254161/149‑1 (MQ=255)
                                                                                                                             aCGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGCTTGCGCGAGTTGATCTCCGACGAg                         >  4:245197/1‑148 (MQ=255)
                                                                                                                             aCGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTACGTATAGAGCTTGCGCGAGTTGATCTCCGACGAg                         <  8:357581/147‑1 (MQ=255)
                                                                                                                                     tATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGCTTGCGCGAGTTGATCTCcga                             <  8:255888/136‑1 (MQ=255)
                                                                                                                                     tATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGCTTGCGCGAGTTGATCTCcga                             >  7:255888/1‑136 (MQ=255)
                                                                                                                                         cGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGCTTGCGCGAGTTGATCTCCGACGAGGAACGAAAATGc             <  1:221588/148‑1 (MQ=255)
                                                                                                                                         cGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAAACGGGTGATCAGGCGCGGCAACGTAGACGGTCATTAACGTATAGAGCTTGCGCGAGTTGATATCCGACGAGGAACGAAAATGcc            >  8:303318/1‑148 (MQ=255)
                                                                                                                                              ccGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGCTTGCGCGAGTTGATCTCCGACGAGGAACGAAAATGCATATcc       >  7:66568/1‑149 (MQ=255)
                                                                                                                                                gCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGCTTGCGCGAGTTGATCTCCGACGAGGAACGAAAATGCATATCCgg     <  3:245197/149‑1 (MQ=255)
                                                                                                                                                   aTTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGCTCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGCTTGCGCGAGTTGATCTCCGACGAGGAACGAAAATGCATATCCGgcgc  >  6:92899/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACCATAGGATACGGCATCGTCGCCACGCCAAGCCCCGCCAGTAATGCCTGACGTTTATCTTCAATCGTGCTCACCGTTAAGCGCGGCTGTTTGTCCAGCAGCTGTACGGTCAACACCGGGCGCTCACGAGCGGTATCCGCCACCGCAATTCCACGATATTTCACACGCGTCACTTCAGATAACGGTTCCGGCTCCTGATGAATCGGGTGATCAGGCGCGGCAACGTAGACGTTCATTAACGTATAGAGCTTGCGCGAGTTGATCTCCGACGAGGAACGAAAATGCATATCCGGCGC  >  NC_000913/3253469‑3253764

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