Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 805,960 A→C 9.5% intergenic (+196/+38) 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,9600AC 9.5% 242.5 / 17.4 84intergenic (+196/+38)ybhJ/ybhCputative hydratase YbhJ/outer membrane lipoprotein YhbC
Reads supporting (aligned to +/- strand):  ref base A (38/38);  new base C (3/5);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.38e-01

CGATGCTTGCGCATCTTATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGTTATCATCT  >  NC_000913/805812‑806109
                                                                                                                                                    |                                                                                                                                                     
cgatgcTTGCGCATCTTATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAg                                                                                                                                                      >  2:1128529/1‑150 (MQ=255)
      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:493641/1‑150 (MQ=255)
                           cGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGt                                                                                                                           <  1:70898/150‑1 (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:1370226/1‑150 (MQ=255)
                                               ggCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCg                                                                                                       >  1:326517/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:368268/1‑150 (MQ=255)
                                                                                     gTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTg                                                                 >  1:1454819/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                                        <  2:996540/150‑1 (MQ=255)
                                                                                                              aTCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGtagttag                                        <  1:281278/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                                                                       >  2:726949/1‑99 (MQ=255)
                                                                                                                                  cGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCgg                                                                       <  1:726949/99‑1 (MQ=255)
                                                                                                                                       tatttgCCGGATGCGGCGTAAACGCCTTATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTtc               >  1:575922/6‑150 (MQ=255)
                                                                                                                                        gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTtcg              <  1:408953/150‑1 (MQ=255)
                                                                                                                                        gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTtcg              <  2:493641/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)
                                                                                                                                                    |                                                                                                                                                     
CGATGCTTGCGCATCTTATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAATAGTTAATTGCTCTTACTTCTTCGCCTCTGCAACCACTTTACTACCCACGCCGCGGTTATTGTATTCCCACATGCGGTTGTAGTTAGTGTCATTCAGATTGCGCTGTATTTCGTCGTTATCATCT  >  NC_000913/805812‑806109

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: