Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,972,795 G→T 6.5% intergenic (‑104/+41) tar ← / ← cheW methyl‑accepting chemotaxis protein II/purine‑binding chemotaxis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,972,7950GT 6.5% 263.5 / 5.6 92intergenic (‑104/+41)tar/cheWmethyl‑accepting chemotaxis protein II/purine‑binding chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base G (43/43);  new base T (2/4);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

CAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCCCAGTCC  >  NC_000913/1972647‑1972943
                                                                                                                                                    |                                                                                                                                                    
cAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATgcgc                                                                                                                                                     >  2:922403/1‑150 (MQ=255)
 aGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCa                                                                                                                                                    >  2:837108/1‑150 (MQ=255)
  gCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAg                                                                                                                                                   >  2:1022450/1‑150 (MQ=255)
     ccATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGatc                                                                                                                                                >  2:704451/1‑150 (MQ=255)
         tACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCg                                                                                                                                            <  1:464466/150‑1 (MQ=255)
          aCCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCgg                                                                                                                                           >  1:512583/1‑150 (MQ=255)
          aCCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGGTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCgg                                                                                                                                           >  1:1013917/1‑150 (MQ=255)
             aaCAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGtt                                                                                                                                        <  1:1587678/150‑1 (MQ=255)
                 gCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCAtt                                                                                                                                    >  1:1337667/1‑150 (MQ=255)
                 gCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGATCTGCGCATCGGGtttttt                                                                                                                                    <  1:646898/150‑5 (MQ=255)
                   gTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTc                                                                                                                                  <  1:704451/150‑1 (MQ=255)
                    tGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCa                                                                                                                                 <  2:551007/150‑1 (MQ=255)
                    tGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCa                                                                                                                                 <  2:544369/150‑1 (MQ=255)
                    tGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGATCTGCGCATCGGGttttttatt                                                                                                                                 <  2:1183950/150‑8 (MQ=255)
                       cTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCaagt                                                                                                                              >  1:839715/1‑148 (MQ=255)
                       cTACGCGGATACGGTTAATCATAAGGCACATTCATGAAAACAAGTTGATCTCGTTATCGCCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAAtt                                                                                                                              >  1:759535/1‑150 (MQ=255)
                           gctgATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAgg                                                                                                                          <  1:1337587/147‑1 (MQ=255)
                                tACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATc                                                                                                                     <  1:1516838/150‑1 (MQ=255)
                                 aCGGTTAATCATAAGGCACCTTCCTGAAAACAAGGTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCg                                                                                                                    <  2:247055/150‑1 (MQ=255)
                                  cGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGTGCGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCgg                                                                                                                   <  1:137633/150‑1 (MQ=255)
                                   ggTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCggg                                                                                                                  <  1:423432/150‑1 (MQ=255)
                                   ggTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGAAGATCATCGGGTTCATTTCAATTGAGGAAATCggg                                                                                                                  >  2:706107/1‑150 (MQ=255)
                                       aaTCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGaa                                                                                                              >  2:1635300/1‑150 (MQ=255)
                                       aaTCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCAGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGaa                                                                                                              >  1:914258/1‑150 (MQ=255)
                                         tCATTAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGTGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAAtt                                                                                                            <  2:355826/150‑1 (MQ=255)
                                             aaGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGc                                                                                                        <  2:1337667/150‑1 (MQ=255)
                                              aGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGcc                                                                                                       >  2:937685/1‑150 (MQ=255)
                                               ggCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCa                                                                                                      >  1:1483377/1‑150 (MQ=255)
                                                     ttCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTg                                                                                                <  2:759535/150‑1 (MQ=255)
                                                         tGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGAcgc                                                                                            >  2:399501/1‑150 (MQ=255)
                                                          gAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGCTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCAGAcgcc                                                                                           >  2:346696/1‑150 (MQ=255)
                                                                  gTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATc                                                                                   <  1:837108/150‑1 (MQ=255)
                                                                   ttGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCt                                                                                  >  1:1072718/1‑150 (MQ=255)
                                                                       tctcGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACa                                                                              <  1:1696786/150‑1 (MQ=255)
                                                                       tctcGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCAACCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATAGCGTTCATTTCAATAGAGGAAAACGGGAGAAATACGCCACTTCTTACGACGAGCTAACTAACa                                                                              >  1:1447062/1‑150 (MQ=255)
                                                                         tcGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAgc                                                                            <  1:922403/150‑1 (MQ=255)
                                                                         tcGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAgc                                                                            >  2:280723/1‑150 (MQ=255)
                                                                            ttATCGGCAATGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGcc                                                                         <  2:1246734/150‑1 (MQ=255)
                                                                            ttATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGcc                                                                         <  2:1353953/150‑1 (MQ=255)
                                                                             tATCGGCAAGGAGGGGGGAAACATTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCa                                                                        <  1:972509/150‑1 (MQ=255)
                                                                                   cAAGTAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCtt                                                                  <  1:706107/150‑1 (MQ=255)
                                                                                   cAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCAGATCCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTTACGCCGCGCTATCTAACAGCGCCATCTCtt                                                                  >  2:327194/1‑150 (MQ=255)
                                                                                    aaGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTc                                                                 <  1:1635300/150‑1 (MQ=255)
                                                                                          gggggAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGtt                                                           >  2:641336/1‑150 (MQ=255)
                                                                                               aaaCTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTcagca                                                      >  2:29635/1‑150 (MQ=255)
                                                                                               aaaCTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTAACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTcagca                                                      >  1:625029/1‑150 (MQ=255)
                                                                                                 aCTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGt                                                    >  1:19348/1‑150 (MQ=255)
                                                                                                 aCTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGAGAAGATCATCGGGTTCATGTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTGAGCAGt                                                    >  1:13361/1‑150 (MQ=255)
                                                                                                  cTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGtt                                                   >  2:217921/1‑150 (MQ=255)
                                                                                                   tttATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGttt                                                  >  1:1541981/1‑150 (MQ=255)
                                                                                                      aTTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGTAATCGGGAGAATTACGCCTCTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGtttttg                                               <  1:346696/150‑2 (MQ=255)
                                                                                                       ttGCTGATGCCACCCGCCGCGAAATTGAAATAATAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGATATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCg                                              <  2:1602783/150‑1 (MQ=255)
                                                                                                          cTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATg                                           >  1:951716/1‑150 (MQ=255)
                                                                                                          cTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCAGTTCAGCAGTTATTCGATg                                           >  1:1198129/1‑150 (MQ=255)
                                                                                                             aTGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTc                                        <  1:280723/150‑1 (MQ=255)
                                                                                                             aTGCCACCCGCCGCGAAAATGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTc                                        <  1:778008/150‑1 (MQ=255)
                                                                                                               ggcagccgcagCGAAATTGAAATAAAAAACCCGATGCTCAGATCATCGGGTTCATTTCAATTGAGGAAATCTTGAGAATTTCGCCACTTCTGACGACGCGCTATATAACAGTGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCAc                                      <  1:1148679/145‑1 (MQ=255)
                                                                                                                 cACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGGAGTTTTTCGATGTTCACCa                                    >  2:818271/1‑150 (MQ=255)
                                                                                                                 cACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCATCGCTGTTCAGCAGTTTTTCGATGTTCACCa                                    <  2:1509702/150‑1 (MQ=255)
                                                                                                                  aCCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAg                                   <  2:1072718/150‑1 (MQ=255)
                                                                                                                  aCCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAg                                   <  1:1209478/150‑1 (MQ=255)
                                                                                                                     cgccgcGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAAt                                <  2:1013917/150‑1 (MQ=255)
                                                                                                                     cgccgcGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAAt                                <  2:512583/150‑1 (MQ=255)
                                                                                                                       ccgcGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCa                              <  1:923061/150‑1 (MQ=255)
                                                                                                                       ccgcGAAATTGAAATAAAAAACCCGATGCGCAGATCAGCGGGTTCAGTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCGTCTCTTCGCTGTTCAGCAGTTGTTCGATGTTCCCCAGAGTCa                              >  1:1049672/1‑150 (MQ=255)
                                                                                                                        cgcgAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCaa                             <  1:1424484/150‑1 (MQ=255)
                                                                                                                         gcgAATTTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAAc                            <  1:865041/150‑1 (MQ=255)
                                                                                                                         gcgAAATTGAAATAAAAAACCCTATGCGCATCTCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTAGGCCTCTTCTGACGCCGCGCTTTCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCTAc                            <  2:1447062/150‑1 (MQ=255)
                                                                                                                          cgAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACa                           <  1:857855/150‑1 (MQ=255)
                                                                                                                             aaTTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATcc                        <  2:1678035/150‑1 (MQ=255)
                                                                                                                                tGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCTCTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGAGGTTCACCAGAATCAACATCCGGt                     <  1:1003520/150‑1 (MQ=255)
                                                                                                                                   aaTAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGc                  >  1:152083/1‑150 (MQ=255)
                                                                                                                                    gaaatgaaCCCGATGATCTGCGCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTc                    >  2:929642/7‑147 (MQ=255)
                                                                                                                                    gaaatgaaCCCGATGATCTGCGCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTc                    <  1:929642/141‑1 (MQ=255)
                                                                                                                                     tAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGccc                >  1:220546/1‑150 (MQ=255)
                                                                                                                                     tAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGccc                <  1:1645190/150‑1 (MQ=255)
                                                                                                                                      aaaaaaCCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGcccc               >  2:1494718/1‑149 (MQ=255)
                                                                                                                                      aaaaaaCCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCa               >  2:905159/1‑150 (MQ=255)
                                                                                                                                      aaaaaaCCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGATGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCAGGTCGCCCa               >  2:1084369/1‑150 (MQ=255)
                                                                                                                                       aaaaaCCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCat              <  2:13361/150‑2 (MQ=255)
                                                                                                                                       aaaaaCCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAg              <  2:19348/150‑1 (MQ=255)
                                                                                                                                        tgaaCCCGATGATCTGCGCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGCCGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGt             >  1:836797/3‑150 (MQ=255)
                                                                                                                                        aaaaCCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGt             >  1:696133/1‑150 (MQ=255)
                                                                                                                                           aCCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTgcg          >  2:1004070/1‑150 (MQ=255)
                                                                                                                                           aCCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCAGCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTgcg          <  1:1022450/150‑1 (MQ=255)
                                                                                                                                            cccGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACAACCGGTCGCCCAGTgcgc         >  2:1179163/1‑150 (MQ=255)
                                                                                                                                             ccGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGcc        >  1:127564/1‑150 (MQ=255)
                                                                                                                                             ccGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGcc        >  2:90010/1‑150 (MQ=255)
                                                                                                                                             ccGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGcc        >  2:933559/1‑150 (MQ=255)
                                                                                                                                             ccGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGcc        <  2:152083/150‑1 (MQ=255)
                                                                                                                                               gATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCcca      >  1:1367027/1‑150 (MQ=255)
                                                                                                                                               gATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCcca      <  1:217921/150‑1 (MQ=255)
                                                                                                                                                 tGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCccagt    >  2:347048/1‑150 (MQ=255)
                                                                                                                                                   cgcAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCccagtcc  <  1:399501/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCCCAGTCC  >  NC_000913/1972647‑1972943

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