Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,851,049 G→T 12.4% D17Y (GAT→TAT)  hypB → hydrogenase isoenzymes nickel incorporation protein HypB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,851,0490GT12.4% 256.2 / 15.6 85D17Y (GAT→TAT) hypBhydrogenase isoenzymes nickel incorporation protein HypB
Reads supporting (aligned to +/- strand):  ref base G (41/33);  new base T (2/9);  total (43/42)
Fisher's exact test for biased strand distribution p-value = 2.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGC  >  NC_000913/2850900‑2851198
                                                                                                                                                     |                                                                                                                                                     
aGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGGAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGtg                                                                                                                                                       >  1:759803/1‑150 (MQ=255)
aGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGtg                                                                                                                                                       >  2:134979/1‑150 (MQ=255)
aGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGtg                                                                                                                                                       <  1:637362/150‑1 (MQ=255)
   gcgTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGtgatg                                                                                                                                                    <  1:985822/150‑1 (MQ=255)
    cgTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGCCGCCGGGTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGGGAAGGCACCCGGTATATAGAGGGtgatga                                                                                                                                                   >  2:60547/1‑150 (MQ=255)
       ccgccgCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGttatgttca                                                                                                                                                <  2:812319/150‑5 (MQ=255)
         gccgcTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATa                                                                                                                                              >  2:1123504/1‑150 (MQ=255)
          ccgcTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGCGATGAACATaa                                                                                                                                             >  1:941610/1‑150 (MQ=255)
               gTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTTATGTTCATCACCCTc                                                                                                                                        <  2:978757/150‑1 (MQ=255)
               gTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTTATGTTCATCACCCTc                                                                                                                                        <  2:1441299/150‑1 (MQ=255)
                 ccACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTTATGTTCATCACCCTcga                                                                                                                                      <  2:1284907/150‑3 (MQ=255)
                 ccACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTTATGTTCATCACCCTcga                                                                                                                                      <  2:952491/150‑3 (MQ=255)
                 ccACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCAt                                                                                                                                      >  2:536735/1‑150 (MQ=255)
                  cacaGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTTATGTTCATCACCCTcgat                                                                                                                                     <  1:823667/150‑4 (MQ=255)
                  cacaGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCAtt                                                                                                                                     >  1:425906/1‑150 (MQ=255)
                  cacaGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCAtt                                                                                                                                     >  1:431398/1‑150 (MQ=255)
                  cacaGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCAtt                                                                                                                                     >  2:852508/1‑150 (MQ=255)
                   acaGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTTATGTTCATCACCCTcgata                                                                                                                                    <  1:76218/150‑5 (MQ=255)
                    caGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTTATGTTCATCACCCTcgatat                                                                                                                                   <  1:1450482/150‑6 (MQ=255)
                      gtgtCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTTATGTTCATCACCCTcgatatac                                                                                                                                 <  2:1257118/150‑8 (MQ=255)
                                gACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAgc                                                                                                                       >  2:1414900/1‑150 (MQ=255)
                                   aTGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAgcgcg                                                                                                                    >  1:738940/1‑150 (MQ=255)
                                    tgctgcAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAgcgcgc                                                                                                                   >  2:547943/1‑150 (MQ=255)
                                     gctgcAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGcc                                                                                                                  >  2:500157/1‑150 (MQ=255)
                                     gctgcAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGcc                                                                                                                  >  1:926772/1‑150 (MQ=255)
                                       tgcAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCAt                                                                                                                <  1:1123504/150‑1 (MQ=255)
                                       tgcAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCAt                                                                                                                >  2:421368/1‑150 (MQ=255)
                                       tgcAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCAt                                                                                                                >  2:278139/1‑150 (MQ=255)
                                          aGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTg                                                                                                             <  2:1314501/150‑1 (MQ=255)
                                               gtgGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCcg                                                                                                        >  1:1148146/1‑150 (MQ=255)
                                               gtgGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCcg                                                                                                        >  2:765042/1‑150 (MQ=255)
                                                 ggCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCcggc                                                                                                      >  1:579124/1‑150 (MQ=255)
                                                  gCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCcggcg                                                                                                     >  2:158417/1‑150 (MQ=255)
                                                   cAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCcggcgg                                                                                                    >  1:1054343/1‑150 (MQ=255)
                                                       cgacgGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACg                                                                                                >  1:662162/1‑150 (MQ=255)
                                                         acgGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGcc                                                                                              >  1:140448/1‑150 (MQ=255)
                                                          cgGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGccc                                                                                             <  2:1235513/150‑1 (MQ=255)
                                                          cgGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGccc                                                                                             <  1:352948/150‑1 (MQ=255)
                                                             tttACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGaa                                                                                          >  2:25833/1‑150 (MQ=255)
                                                              ttACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAg                                                                                         <  1:852508/150‑1 (MQ=255)
                                                              ttACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAg                                                                                         <  2:1415897/150‑1 (MQ=255)
                                                                aCAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGAt                                                                                       >  1:245391/1‑150 (MQ=255)
                                                                   gATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGaa                                                                                    >  2:809754/1‑150 (MQ=255)
                                                                     ttCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGaaaa                                                                                  <  1:959166/150‑1 (MQ=255)
                                                                     ttCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGaaaa                                                                                  <  1:1453092/150‑1 (MQ=255)
                                                                       cggcggATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATc                                                                                <  1:1178077/150‑1 (MQ=255)
                                                                           ggATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCg                                                                            <  1:1247964/150‑1 (MQ=255)
                                                                             aTAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGc                                                                          <  2:426780/150‑1 (MQ=255)
                                                                                    tAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAg                                                                   <  2:738940/150‑1 (MQ=255)
                                                                                     aGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATACCCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAgc                                                                  >  1:1257979/1‑150 (MQ=255)
                                                                                     aGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAgc                                                                  >  1:1263876/1‑150 (MQ=255)
                                                                                      gACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAgcg                                                                 >  2:1377742/1‑150 (MQ=255)
                                                                                              gTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAAtt                                                         <  2:1391483/150‑1 (MQ=255)
                                                                                               tGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAAttt                                                        <  2:140448/150‑1 (MQ=255)
                                                                                                gAGCGATGTGTACAACATGCGGTTGCTTTGAAGGCAACCTTTATATCGAGGGTGATTAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGTCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTa                                                       <  1:1484747/150‑1 (MQ=255)
                                                                                                gAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTa                                                       >  1:1074297/1‑150 (MQ=255)
                                                                                                gAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTa                                                       >  1:1086942/1‑150 (MQ=255)
                                                                                                   cGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTAccc                                                    <  1:328460/150‑1 (MQ=255)
                                                                                                       gtgtACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCccag                                                >  1:1404463/1‑150 (MQ=255)
                                                                                                            caacaTGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGa                                           >  2:1225165/1‑150 (MQ=255)
                                                                                                            caacaTGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGa                                           >  2:917371/1‑150 (MQ=255)
                                                                                                              acaTGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACt                                         <  2:175972/150‑1 (MQ=255)
                                                                                                                 tgcggttgcggtGAAGGCAACCTGTATATAGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTgaa                                      >  1:858760/1‑150 (MQ=255)
                                                                                                                  gcggttgcggtGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTgaag                                     <  2:1086942/150‑1 (MQ=255)
                                                                                                                   cggttgcggtGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTgaaga                                    <  2:729166/150‑1 (MQ=255)
                                                                                                                    ggttgcggtGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTgaagaa                                   <  2:1054343/150‑1 (MQ=255)
                                                                                                                        gcggtGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCg                               <  1:1189758/150‑1 (MQ=255)
                                                                                                                          ggtGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGATGAAGGCGAc                             <  2:245391/150‑1 (MQ=255)
                                                                                                                          ggtGAAGGCAACCTGTATATCGAGGGTGATGAACATACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCccag                                                >  1:365246/1‑130 (MQ=255)
                                                                                                                          ggtGAAGGCAACCTGTATATCGAGGGTGATGAACATACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCccag                                                <  2:365246/130‑1 (MQ=255)
                                                                                                                           gtGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGAcc                            <  1:1500450/150‑1 (MQ=255)
                                                                                                                             gAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTg                          >  2:642078/1‑150 (MQ=255)
                                                                                                                               aGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTGGTAGCGCGCCATTTGCCCCGGCGGCACGCCGGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCa                        <  1:158417/150‑1 (MQ=255)
                                                                                                                               aGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCa                        <  1:321087/150‑1 (MQ=255)
                                                                                                                               aGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCa                        >  1:1103334/1‑150 (MQ=255)
                                                                                                                               aGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCa                        <  1:809754/150‑1 (MQ=255)
                                                                                                                               aGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCa                        >  1:466891/1‑150 (MQ=255)
                                                                                                                                  cAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTa                     <  2:1394451/150‑1 (MQ=255)
                                                                                                                                   aaCCTGTATATCGAGGGTGATGAAGATAACCCTCATTCCGCGTTTCGTAGCGCGCCTTTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTAc                    <  2:1267743/150‑1 (MQ=255)
                                                                                                                                     ccTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACgg                  <  2:926772/150‑1 (MQ=255)
                                                                                                                                       tGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTc                >  1:1205914/1‑150 (MQ=255)
                                                                                                                                          ggaatgAGGGTTATGTTCATCACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATg             >  2:853121/6‑150 (MQ=255)
                                                                                                                                          atatCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATg             <  2:662162/150‑1 (MQ=255)
                                                                                                                                           tatCGAGGGTTATGATCATTACCCTCATTCCGCGTTTCGTAGCGCGCGATTTGCCGCGGCGGCACGCCCGAAGATGAATATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTAGGGTCTTgg            <  2:1089409/150‑1 (MQ=255)
                                                                                                                                           tatCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATgg            <  1:1414900/150‑1 (MQ=255)
                                                                                                                                            atCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGc           >  2:53210/1‑150 (MQ=255)
                                                                                                                                            atCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGc           <  1:1377742/150‑1 (MQ=255)
                                                                                                                                            atCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGc           >  1:1591604/1‑150 (MQ=255)
                                                                                                                                              ctaGGGTGATGAACATAACCCTCATTCAGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGa         >  2:900657/3‑150 (MQ=255)
                                                                                                                                               gAGGGTTATGTTCATCACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGaa        >  1:1164136/1‑150 (MQ=255)
                                                                                                                                                  ggTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCg     >  1:899342/1‑150 (MQ=255)
                                                                                                                                                    tgatgaACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCggg   >  2:113505/1‑150 (MQ=255)
                                                                                                                                                     gatgaACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGc  <  2:497300/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCGCGTCCGCCGCTGTCCACAGTGTCATGGTGACATGCTGCAGATTGTGGCAGACGACGGTTTACAGATTCGGCGGATAGAAATAGACCAGGAGTGAGCGATGTGTACAACATGCGGTTGCGGTGAAGGCAACCTGTATATCGAGGGTGATGAACATAACCCTCATTCCGCGTTTCGTAGCGCGCCATTTGCCCCGGCGGCACGCCCGAAGATGAAAATCACCGGCATTAAAGCGCCTGAATTTACCCCCAGCCAGACTGAAGAAGGCGACCTGCATTACGGTCATGGCGAAGCGGGC  >  NC_000913/2850900‑2851198

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