Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 805,967 T→A 9.2% intergenic (+203/+31) ybhJ → / ← ybhC putative hydratase YbhJ/outer membrane lipoprotein YhbC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913805,9670TA 9.2% 257.1 / 10.1 85intergenic (+203/+31)ybhJ/ybhCputative hydratase YbhJ/outer membrane lipoprotein YhbC
Reads supporting (aligned to +/- strand):  ref base T (40/37);  new base A (3/5);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TTGCGCATCTTATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGTTATCATCTACGCTGC  >  NC_000913/805818‑806116
                                                                                                                                                     |                                                                                                                                                     
ttGCGCATCTTATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGttt                                                                                                                                                       <  1:261108/150‑1 (MQ=255)
 tGCGCATCTTATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTa                                                                                                                                                      >  1:1435037/1‑150 (MQ=255)
        cTTATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCa                                                                                                                                               >  2:281278/1‑150 (MQ=255)
           aTCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                            <  1:448009/150‑1 (MQ=37)
             ccGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                          <  1:1128529/150‑1 (MQ=255)
             ccGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                          <  1:75516/150‑1 (MQ=255)
              cGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                         <  1:399201/150‑1 (MQ=37)
               gACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                        <  1:938866/150‑2 (MQ=37)
                 ccTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCaaa                                                                                                                                      >  2:382355/1‑150 (MQ=255)
                 ccTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaaa                                                                                                                                      <  1:66675/150‑1 (MQ=255)
                   tACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATa                                                                                                                                    <  2:1085198/150‑1 (MQ=255)
                    aCGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcctacgg                                                                                                                                   <  2:245503/150‑7 (MQ=25)
                    aCGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAg                                                                                                                                   <  1:1452295/150‑1 (MQ=255)
                     cGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGt                                                                                                                                  <  1:70898/150‑1 (MQ=255)
                     cGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGt                                                                                                                                  >  1:493641/1‑150 (MQ=255)
                       aaTCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTa                                                                                                                                <  1:1122358/150‑1 (MQ=255)
                              tCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGctc                                                                                                                         >  1:1179090/1‑150 (MQ=255)
                                gAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCtt                                                                                                                       >  2:332662/1‑150 (MQ=255)
                                    cTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTActt                                                                                                                   >  2:971156/1‑150 (MQ=255)
                                     tGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTActtg                                                                                                                  >  2:719234/1‑149 (MQ=255)
                                     tGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTActtc                                                                                                                  >  1:862698/1‑150 (MQ=255)
                                     tGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTActtc                                                                                                                  >  2:584711/1‑150 (MQ=255)
                                     tGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTActta                                                                                                                  >  1:1289388/1‑149 (MQ=255)
                                     tGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTActta                                                                                                                  >  1:1483650/1‑149 (MQ=255)
                                      gTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTActtct                                                                                                                 <  1:335863/150‑1 (MQ=255)
                                         ggCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCg                                                                                                              >  1:326517/1‑150 (MQ=255)
                                         ggCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCg                                                                                                              >  1:1370226/1‑150 (MQ=255)
                                           ccAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGcc                                                                                                            >  1:560057/1‑150 (MQ=255)
                                            cAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCCGATGCGATTCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTACGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGAct                                                                                                           >  2:746787/1‑150 (MQ=255)
                                              gATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCctct                                                                                                         >  2:303815/1‑150 (MQ=255)
                                               aTAAGGCATTTTCGCATCATCCTTCACTTTTTGTCGGATGAGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGAGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTg                                                                                                        <  2:1370226/150‑1 (MQ=37)
                                                 aaGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCa                                                                                                      >  1:349882/1‑150 (MQ=255)
                                                           cGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTtac                                                                                            <  2:1496160/150‑1 (MQ=255)
                                                            gcagcaTCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTtact                                                                                           <  2:1435037/150‑1 (MQ=255)
                                                            atcgcaTCCGACAATAAGTGCCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTtact                                                                                           >  1:985791/4‑150 (MQ=37)
                                                               gcaTCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTAcc                                                                                        >  1:691500/1‑150 (MQ=255)
                                                                   ccGGCACTTATTGTCGGATGCGATGCTTTCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCAcg                                                                                    >  2:551438/1‑150 (MQ=255)
                                                                   ccGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCAcg                                                                                    >  2:478289/1‑150 (MQ=255)
                                                                     ggCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCAcgcc                                                                                  >  2:349050/1‑150 (MQ=255)
                                                                           tATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGtt                                                                            >  1:866926/1‑150 (MQ=255)
                                                                               gTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAACTCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTg                                                                        <  2:976198/150‑1 (MQ=255)
                                                                               gTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTg                                                                                                        >  2:898987/1‑118 (MQ=255)
                                                                               gTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTg                                                                                                        <  1:898987/118‑1 (MQ=255)
                                                                               gTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTg                                                                        <  1:63291/150‑1 (MQ=255)
                                                                               gTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTg                                                                        >  1:1454819/1‑150 (MQ=255)
                                                                               gTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTg                                                                        >  1:368268/1‑150 (MQ=255)
                                                                                    atgcgatgcTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTc                                                                   <  1:660390/150‑1 (MQ=255)
                                                                                     tgcgatgcTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGt                                                                       <  2:900418/145‑1 (MQ=255)
                                                                                     tgcgatgcTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGt                                                                       >  1:900418/1‑145 (MQ=255)
                                                                                        gatgcTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCcac                                                               >  1:924147/1‑150 (MQ=255)
                                                                                             ttGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCg                                                          <  1:332662/150‑1 (MQ=255)
                                                                                               gcgcATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGt                                                        <  1:1434449/150‑1 (MQ=255)
                                                                                                cgcATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGtt                                                       >  1:631125/1‑150 (MQ=255)
                                                                                                  cATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGt                                                     <  2:691500/150‑1 (MQ=255)
                                                                                                   aTCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGta                                                    <  1:437245/150‑1 (MQ=255)
                                                                                                    tCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGtag                                                   >  1:231653/1‑150 (MQ=255)
                                                                                                        aTCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGtagttag                                               <  1:281278/150‑1 (MQ=255)
                                                                                                        aTCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGtagttag                                               <  2:996540/150‑1 (MQ=255)
                                                                                                        aTCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGtagttag                                               >  2:111328/1‑150 (MQ=255)
                                                                                                          ccGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAgtg                                             <  2:1289388/150‑1 (MQ=255)
                                                                                                          ccGACCTACAAATCGCATCGAACCGTACGCCGGATAAGGCGGTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTGGTAGTTAgtg                                             <  2:886146/150‑1 (MQ=255)
                                                                                                           cGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCTTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAgtgt                                            <  1:551438/150‑1 (MQ=255)
                                                                                                               cTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCAtt                                        <  1:448567/150‑1 (MQ=255)
                                                                                                               cTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCAtt                                        >  2:1065405/1‑150 (MQ=255)
                                                                                                                tACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTc                                       <  2:1506521/150‑1 (MQ=255)
                                                                                                                tACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTc                                       >  1:1196854/1‑150 (MQ=255)
                                                                                                                        gCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTgcg                               >  1:309640/1‑150 (MQ=255)
                                                                                                                          aTCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCt                             <  1:1509962/150‑1 (MQ=255)
                                                                                                                            cGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCgg                                                                              <  1:726949/99‑1 (MQ=255)
                                                                                                                            cGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCgg                                                                              >  2:726949/1‑99 (MQ=255)
                                                                                                                                 tatttgCCGGATGCGGCGTAAACGCCTTATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTtc                      >  1:575922/6‑150 (MQ=255)
                                                                                                                                  gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTtcg                     <  2:493641/150‑1 (MQ=255)
                                                                                                                                  gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTtcg                     <  1:408953/150‑1 (MQ=255)
                                                                                                                                     tgCCGGATGCGGCGTAAACGCCTTATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTg                                                                        >  1:345556/2‑96 (MQ=37)
                                                                                                                                     tgCCGGATGCGGCGTAAACGCCTTATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTg                                                                        <  2:345556/95‑1 (MQ=37)
                                                                                                                                     ggCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTtcgtcg                  <  2:368268/150‑1 (MQ=255)
                                                                                                                                     ggCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTtcgtcg                  <  2:1080659/150‑1 (MQ=255)
                                                                                                                                      gCCGGATGCGGCGTAAACGCCTTATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTcgtcgt                 >  1:1257016/1‑150 (MQ=255)
                                                                                                                                       cctgATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGtt                <  2:560057/147‑1 (MQ=255)
                                                                                                                                       ccGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGtt                >  1:323935/1‑150 (MQ=255)
                                                                                                                                        cGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGTTa               <  2:231653/150‑1 (MQ=255)
                                                                                                                                          gATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGTTatc             <  1:971156/150‑1 (MQ=255)
                                                                                                                                           aTAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGTTatca            <  1:349050/150‑1 (MQ=255)
                                                                                                                                              aGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGTTATCATCt         >  1:323260/1‑150 (MQ=255)
                                                                                                                                                gCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGTTATCATCTAc       >  1:1472161/1‑150 (MQ=255)
                                                                                                                                                 cGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGTTATCATCTACg      >  1:923824/1‑150 (MQ=255)
                                                                                                                                                   tttACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTATGCAACCACTTTACTACCCACGCCGCGGTTATTGTAATCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGTTATCATCTACGCt    >  2:1191955/1‑150 (MQ=255)
                                                                                                                                                     tACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGTTATCATCTACGCTGc  <  1:111328/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTGCGCATCTTATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGTTATCATCTACGCTGC  >  NC_000913/805818‑806116

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