Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,972,799 A→C 8.4% intergenic (‑108/+37) 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,7990AC 8.4% 186.5 / 8.0 72intergenic (‑108/+37)tar/cheWmethyl‑accepting chemotaxis protein II/purine‑binding chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base A (36/30);  new base C (4/2);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.20e-01

TACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCCCAGTCCAGTG  >  NC_000913/1972656‑1972947
                                                                                                                                               |                                                                                                                                                    
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)
                                                                                                                                              gATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACCGCGCCATCTCTTCGCTGTTCAGCAGTTATTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCCCACTCCAGTg  >  2:704184/1‑150 (MQ=255)
                                                                                                                                              gATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCCCAGTCCAGTg  <  1:761057/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TACCAACAGCGTGACTACGCGGATACGGTTAATCATAAGGCACCTTCCTGAAAACAAGTTGATCTCGTTATCGGCAAGGAGGGGGGAAACTTTATTGCTGATGCCACCCGCCGCGAAATTGAAATAAAAAACCCGATGCGCAGATCATCGGGTTCATTTCAATTGAGGAAATCGGGAGAATTACGCCACTTCTGACGCCGCGCTATCTAACAGCGCCATCTCTTCGCTGTTCAGCAGTTTTTCGATGTTCACCAGAATCAACATCCGGTCGCCCAGTGCGCCCAGTCCAGTG  >  NC_000913/1972656‑1972947

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: