Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,972,795 G→T 9.1% 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 9.1% 172.7 / 8.0 66intergenic (‑104/+41)tar/cheWmethyl‑accepting chemotaxis protein II/purine‑binding chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base G (32/28);  new base T (4/2);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

CCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCCCAGTCCA  >  NC_000913/1972646‑1972944
                                                                                                                                                     |                                                                                                                                                     
ccAGCACCATTACCAACAGCGTGACTACGCGTATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGTGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATgcg                                                                                                                                                       <  2:1221506/150‑1 (MQ=255)
  aGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCa                                                                                                                                                     <  2:431188/150‑1 (MQ=255)
          tACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTTATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCg                                                                                                                                             <  1:560057/150‑1 (MQ=255)
             cAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGt                                                                                                                                          <  2:486887/150‑1 (MQ=255)
              aaCAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGtt                                                                                                                                         >  2:807489/1‑150 (MQ=255)
                 aGCGTGACTTCGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCAt                                                                                                                                      <  1:161915/150‑1 (MQ=255)
                   cGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGCTGCGCAGATCATCGGGTTCAttt                                                                                                                                    >  1:1073373/1‑150 (MQ=255)
                   cGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCAttt                                                                                                                                    >  1:611945/1‑150 (MQ=255)
                   cGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGATCTGCGCATCGGGtttttta                                                                                                                                    <  2:993368/150‑6 (MQ=255)
                   cGTGACTACGCGGATACGGTAAATCATCAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCGGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGATCTGCGCATCGGGtttttta                                                                                                                                    <  2:639170/150‑6 (MQ=255)
                    gTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCAttta                                                                                                                                   >  2:1316389/1‑149 (MQ=255)
                    gTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTc                                                                                                                                   <  1:720862/150‑1 (MQ=255)
                    gTGACTAAGCGGATAAGGTTAATCATAAGGCACCTTCCTTAAAACAAGTTTATCTCGTTATCGGCAAGGAGGGTGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTc                                                                                                                                   <  1:869860/150‑1 (MQ=255)
                     tGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGAAGATCATCGGGTTCATTTCa                                                                                                                                  >  1:545511/1‑150 (MQ=255)
                      gACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCaa                                                                                                                                 >  2:557323/1‑150 (MQ=255)
                        cTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGACGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAAtt                                                                                                                               >  1:794704/1‑150 (MQ=255)
                           cgcgGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAg                                                                                                                            >  1:519192/1‑150 (MQ=255)
                            gcgGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAgg                                                                                                                           <  2:431207/150‑1 (MQ=255)
                                   cGGTTAATCATAAGGCACCTTCCTGAAACCAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCgg                                                                                                                    <  1:1029149/150‑1 (MQ=255)
                                   cGGGTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCgg                                                                                                                    <  2:589365/150‑1 (MQ=255)
                                    ggTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCggg                                                                                                                   >  1:545531/1‑150 (MQ=255)
                                         gtCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAAt                                                                                                              >  1:762616/2‑150 (MQ=255)
                                         aTCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTAATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGGCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAAt                                                                                                              <  1:1007042/150‑1 (MQ=255)
                                          tCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTACATTTCAATTGAGGAAATCGGGAGAAtt                                                                                                             >  1:893048/1‑150 (MQ=255)
                                             aaaGGCACCTTCCTGTAAACAAGTTGATCTTGATATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAATTTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGATGAAATCGGGAGAATTACg                                                                                                          <  2:663800/149‑1 (MQ=255)
                                                   aCCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACtt                                                                                                    <  2:1314150/150‑1 (MQ=255)
                                                    ccTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACtta                                                                                                   >  2:551208/1‑149 (MQ=255)
                                                    ccTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTc                                                                                                   >  2:1058603/1‑150 (MQ=255)
                                                      ttCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCTCCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTg                                                                                                 <  1:992086/150‑1 (MQ=255)
                                                        ccTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGAc                                                                                               >  2:1154634/1‑150 (MQ=255)
                                                        ccTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATAGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGAc                                                                                               >  2:1144547/1‑150 (MQ=255)
                                                         atGAAATCAAGTTGATCTCGTTATCGGCTAGGAGGGTGGATACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGTTGGGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGAcg                                                                                              <  1:1105165/149‑1 (MQ=255)
                                                           gAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGAGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGAcgcc                                                                                            >  2:966318/1‑150 (MQ=255)
                                                                   gtggATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATAGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATc                                                                                    <  1:1114595/147‑1 (MQ=255)
                                                                    ttGATCTCGTTATCGGCAAGGAGGGGTGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTTAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCt                                                                                   <  2:595076/150‑1 (MQ=255)
                                                                     tGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTa                                                                                  <  2:996066/150‑1 (MQ=255)
                                                                            gTTATCGGCATGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACTCCGCGCTATCTAACAgcgc                                                                           <  2:301016/150‑1 (MQ=255)
                                                                               aTCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCAt                                                                        >  1:1097089/1‑150 (MQ=255)
                                                                               aTCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGAAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGACAt                                                                        >  1:213624/1‑150 (MQ=255)
                                                                                  ggCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCAtctc                                                                     >  2:500515/1‑150 (MQ=255)
                                                                                  ggCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCAtctc                                                                     <  2:946414/150‑1 (MQ=255)
                                                                                  ggCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCAtcta                                                                     >  1:870276/1‑149 (MQ=255)
                                                                                  ggCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGAGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCAtctc                                                                     >  2:252563/1‑150 (MQ=255)
                                                                                       ggaggGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGc                                                                <  2:81308/150‑1 (MQ=255)
                                                                                          gtggggAAAAGTTATTGCTGATGCCACCCGCCGCGAAATTGAACTAAAAAACCCGCTGCGCAGAGCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACATCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGt                                                             <  1:967175/148‑1 (MQ=255)
                                                                                                aatcTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTcagca                                                       <  2:997831/147‑1 (MQ=255)
                                                                                                aaaCTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTcagca                                                       <  2:555869/150‑1 (MQ=255)
                                                                                                   cTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGtt                                                    >  1:47566/1‑150 (MQ=255)
                                                                                                          gCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGTGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGAt                                             >  2:769585/1‑150 (MQ=255)
                                                                                                          gCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGAt                                             >  1:918916/1‑150 (MQ=255)
                                                                                                           cTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATg                                            >  1:516392/1‑150 (MQ=255)
                                                                                                             gATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGtt                                          >  2:278816/1‑150 (MQ=255)
                                                                                                                  cACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCACTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCa                                     <  1:902946/150‑1 (MQ=255)
                                                                                                                     ccgccgCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGaa                                  >  1:1238187/1‑150 (MQ=255)
                                                                                                                          gcgAAATTGAAATAAAAAACCCTATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGATATCTAACAGCGACATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAAc                             <  2:1263168/150‑1 (MQ=255)
                                                                                                                          gcgAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAAc                             >  2:1311958/1‑150 (MQ=255)
                                                                                                                          gcgAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAAc                             >  2:694774/1‑150 (MQ=255)
                                                                                                                           cgAAATTGAAATAAAAAACCCGCTGCGCAGATCATCGGGTTAATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACa                            >  1:864344/1‑150 (MQ=255)
                                                                                                                              aaTTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATcc                         >  2:256101/1‑150 (MQ=255)
                                                                                                                                     gaaatgaaCCCGATGATCTGCGCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGcc                  >  2:505547/7‑150 (MQ=255)
                                                                                                                                     gaaatgaaCCCGATGATCTGCGCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGcc                  >  2:1000411/7‑150 (MQ=255)
                                                                                                                                     gaaatgaaCCCGATGATCTGCGCATCGGGTACATTTCAATTGACGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCACTTTTTCGATGTTCACCAGAATCAACATCCGGTCGcc                  >  2:880329/7‑150 (MQ=255)
                                                                                                                                         aaaaCCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGt              <  1:252563/150‑1 (MQ=255)
                                                                                                                                          gaaCCCGATGATCTGCGCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATGTAACAGCGCCATCTCTTCGCTGATCAGAAGTTTTGCGATGTACACCAGAATCAACGTCCGGTCGCGCAGTg             >  2:709656/2‑150 (MQ=255)
                                                                                                                                           aaCCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTgc            <  1:1154634/150‑1 (MQ=255)
                                                                                                                                            aCCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTgcg           <  1:557323/150‑1 (MQ=255)
                                                                                                                                              ccGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGcc         >  1:886023/1‑150 (MQ=255)
                                                                                                                                                     gcAGATAATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATACGGTCGCCCAGTGCGCccagtcca  >  1:414713/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCCCAGTCCA  >  NC_000913/1972646‑1972944

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: