Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,972,794 C→A 11.1% intergenic (‑103/+42) 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,7940CA11.1% 145.5 / 8.1 66intergenic (‑103/+42)tar/cheWmethyl‑accepting chemotaxis protein II/purine‑binding chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base A (6/2);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

CCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCC  >  NC_000913/1972645‑1972937
                                                                                                                                                     |                                                                                                                                               
cccAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATgc                                                                                                                                                 >  2:902946/1‑150 (MQ=255)
 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:946414/150‑1 (MQ=255)
                                                                                   ggCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCAtctc                                                              >  2:500515/1‑150 (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:694774/1‑150 (MQ=255)
                                                                                                                           gcgAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAAc                      >  2:1311958/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)
                                                                                                                                                     |                                                                                                                                               
CCCAGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCC  >  NC_000913/1972645‑1972937

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: