Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,972,797 A→T 9.0% intergenic (‑106/+39) 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,7970AT 9.0% 171.8 / 6.6 67intergenic (‑106/+39)tar/cheWmethyl‑accepting chemotaxis protein II/purine‑binding chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base A (33/28);  new base T (4/2);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

AGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCCCAGTCCAGT  >  NC_000913/1972648‑1972946
                                                                                                                                                     |                                                                                                                                                     
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)
                                                                                                                                                    cAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCccagtccag   <  1:1144547/150‑1 (MQ=255)
                                                                                                                                                     aGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCccagtccagt  >  1:264018/1‑150 (MQ=255)
                                                                                                                                                     aGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCAACATCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTATTTCTATGTTCACCATAATCCACATCAGGTCGCCCAGTGAGCccagtccagt  >  1:232730/1‑150 (MQ=255)
                                                                                                                                                     aGATCATCGGGTACATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCccagtccagt  <  2:893048/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCACCATTACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCCCAGTCCAGT  >  NC_000913/1972648‑1972946

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: