Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,056,038 T→G 8.1% D101A (GAC→GCC)  glnL ← sensory histidine kinase in two‑component regulatory system with GlnG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,056,0380TG 8.1% 224.2 / 5.7 90D101A (GAC→GCC) glnLsensory histidine kinase in two‑component regulatory system with GlnG
Reads supporting (aligned to +/- strand):  ref base T (38/44);  new base G (3/5);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAA  >  NC_000913/4055891‑4056186
                                                                                                                                                   |                                                                                                                                                    
ccACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTc                                                                                                                                                     >  2:888773/1‑149 (MQ=255)
ccACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTc                                                                                                                                                     <  1:546331/149‑1 (MQ=255)
 cACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTACGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTcc                                                                                                                                                    >  2:1188267/1‑149 (MQ=255)
   cGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCgg                                                                                                                                                  >  2:205518/1‑149 (MQ=255)
        aCCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCCGCAGGATCATGCCGTCCGGCATAc                                                                                                                                             >  1:13905/1‑149 (MQ=255)
          cGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTAAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGGACGGCAtgatc                                                                                                                                           <  2:460180/149‑5 (MQ=255)
            ccAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGGACGGCAtgatcct                                                                                                                                         <  1:1292800/149‑7 (MQ=255)
                  ggATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACTGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGt                                                                                                                                   <  1:1561876/149‑1 (MQ=255)
                   gATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTc                                                                                                                                  <  1:1334259/149‑1 (MQ=255)
                     tttttAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTcac                                                                                                                                <  1:1525136/149‑1 (MQ=255)
                      ttttAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTcaca                                                                                                                               <  1:1300626/149‑1 (MQ=255)
                          aaTCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGaga                                                                                                                           <  1:483874/149‑3 (MQ=255)
                          aaTCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGaaa                                                                                                                           <  1:1571030/149‑1 (MQ=255)
                           aTCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCCACCTCCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGTAGCCATCTCCAGCAGGAGCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAg                                                                                                                          <  1:1273523/149‑1 (MQ=255)
                           aTCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAg                                                                                                                          <  1:391625/149‑1 (MQ=255)
                            tctcATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAgg                                                                                                                         >  1:763378/1‑149 (MQ=255)
                              tcATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGat                                                                                                                       <  2:84435/149‑1 (MQ=255)
                               cATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGATCCTGACTTATGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCACGCCGTCCGGCATCCGCTGGGCCGTCACAGCAAGGata                                                                                                                      >  2:277100/1‑149 (MQ=255)
                                    gCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAATGATATGCGa                                                                                                                 <  1:904732/149‑1 (MQ=255)
                                          ccGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGccc                                                                                                           >  1:1384944/1‑149 (MQ=255)
                                            gcgcTCTAGATCCCGGGCAGCAACCTGCTGGGCGGGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCAAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGt                                                                                                         <  2:1075997/149‑1 (MQ=255)
                                              gcACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAACGAGATGCGCGCGCCCGtgg                                                                                                       >  2:378229/1‑147 (MQ=255)
                                                aCTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAt                                                                                                     >  2:609961/1‑149 (MQ=255)
                                                aCTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAt                                                                                                     >  1:476213/1‑149 (MQ=255)
                                                aCTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGAt                                                                                                     >  1:4621/1‑149 (MQ=255)
                                                   aaaTCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGACGATCAc                                                                                                  >  1:392622/1‑149 (MQ=255)
                                                    aaTCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGAcc                                                                                                 <  1:555153/149‑1 (MQ=255)
                                                      tCAGGGGGAGCAACCTGCTGGGCTTGCTGTTGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCTTCGGAGCCATCTCCAGGAGGATCCTGCCGGCCGGCATACGCTGGGCCGTCAGAGAAAGGATATGCGAGCGCCCGTCGATGACCAg                                                                                               <  2:1359730/149‑1 (MQ=255)
                                                      tCACGGGCAGCAACCTGCTGGGCGTGCTGTAGGTGTTCCTGCCTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCATGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGATCGCCCGTCGATGACCAg                                                                                               >  1:707092/1‑149 (MQ=255)
                                                       cACGGGCAGCACCCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGc                                                                                              <  1:875/149‑1 (MQ=255)
                                                       cACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGc                                                                                              >  1:41388/1‑149 (MQ=255)
                                                              cgcaACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCCCAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACtt                                                                                       >  2:312085/2‑149 (MQ=255)
                                                              agcaACCTGCTGGGCGTGCAGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATATCCAGCAGGATCATGCCGTCAGGCATACCCTGGGGCGTCACAGCAAGGATATGCGAGCGCCAGTCGCTGACCCGCGTCACtg                                                                                       >  1:1121745/1‑148 (MQ=255)
                                                               gcaACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTc                                                                                      <  1:888773/149‑1 (MQ=255)
                                                                caACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCTTAACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCg                                                                                     >  1:928348/1‑149 (MQ=255)
                                                                caACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCg                                                                                     >  2:1244627/1‑149 (MQ=255)
                                                                caACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCg                                                                                     >  1:126090/1‑149 (MQ=255)
                                                                  cccTGCTGGGCGTGCTGTCGCTGTTCCGTACTTAAGCGGCGCTGGTTATCCATCCGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAACGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGtt                                                                                   <  2:1085163/148‑1 (MQ=255)
                                                                    gtgctgGGCGTGCTGTAGCTGTTCCTGACTTACGCGGCGCGGGTTATCCATCGTAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTAt                                                                                 <  2:1515064/148‑1 (MQ=255)
                                                                      gctgGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGGGCTGGTTATCCATCGGAGCCATCTCCAGCATGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGTATATGCGAGAGCCCGTCGATGACCAGCGTCACTTCGTTATCg                                                                               <  2:1228407/149‑1 (MQ=255)
                                                                      gctgGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCg                                                                               <  1:1143102/149‑1 (MQ=255)
                                                                      gctgGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGACCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCg                                                                               <  1:459229/149‑1 (MQ=255)
                                                                        tgGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGt                                                                             <  2:715410/149‑1 (MQ=255)
                                                                        tgGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGt                                                                             <  1:258008/149‑1 (MQ=255)
                                                                        tgGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGt                                                                             >  1:808265/1‑149 (MQ=255)
                                                                          ggCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTaa                                                                           >  1:1095544/1‑149 (MQ=255)
                                                                                  gTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGCCCAGCGTCACTTCGTTATCGGTAAAACCTTGc                                                                   >  1:755714/1‑149 (MQ=255)
                                                                                  gTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGc                                                                   >  1:547826/1‑149 (MQ=255)
                                                                                    aGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGccc                                                                 <  2:547826/149‑1 (MQ=255)
                                                                                     gCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCTTCTCCAGCAGGATCATGCCGTCCTGCATACGCTGGGCCGTCACAGAATGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGcccc                                                                <  1:395714/149‑1 (MQ=255)
                                                                                        gTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGcc                                                             >  2:1301479/1‑149 (MQ=255)
                                                                                            cTGACTTAAGCTGCGCTGGTTATCCATCGGAGCCATCTCCTGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACGTCGTGAACGGGAAAACCTTGCCCGGCGTCCa                                                         >  2:1075004/1‑149 (MQ=255)
                                                                                              gACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGa                                                       <  2:1384944/149‑1 (MQ=255)
                                                                                                   aaGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGTCCGTCACAGAAAGGATCTGCGAGCGCCCGTCGATGACCCGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTc                                                  >  2:1433283/1‑149 (MQ=255)
                                                                                                   aaGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTc                                                  <  2:860985/149‑1 (MQ=255)
                                                                                                   aaGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGACACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTc                                                  <  1:704369/149‑1 (MQ=255)
                                                                                                    aGCGGCGCTGGGTATCCATCGGAGCCAGCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCt                                                 <  1:727377/149‑1 (MQ=255)
                                                                                                     gcggcgCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCtt                                                <  1:828481/149‑1 (MQ=255)
                                                                                                       ggcgCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGc                                              <  1:1454803/149‑1 (MQ=255)
                                                                                                       ggcgCGGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGc                                              <  1:884058/149‑1 (MQ=255)
                                                                                                            tGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAg                                         <  2:425031/149‑1 (MQ=255)
                                                                                                              gTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCt                                       >  2:893886/1‑149 (MQ=255)
                                                                                                                  tCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGAGTCACTTCGTTATCGGTAACACCTTGCCCCGCCTCCAGCCTTTCTTGCATCAGCTCGat                                   >  1:446110/1‑149 (MQ=255)
                                                                                                                  tCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTAACAGAAAGGATATGCGAGCGCCCGTCGATGACCCGCGTCACTTCGATATCGGTAAAACCTTGCCCCGCCACCAGACTTTATTGCATCAGCTCGaa                                   >  2:1170672/1‑148 (MQ=255)
                                                                                                                  tCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCAGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTCCATCAGCTCGat                                   <  1:394535/149‑1 (MQ=255)
                                                                                                                   ccATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGata                                  >  1:1173657/1‑149 (MQ=255)
                                                                                                                   ccATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGata                                  <  1:1549833/149‑1 (MQ=255)
                                                                                                                    cATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTGGTGATCGGTAAAACCTTGCCCCGCCTCCAGCCTTTCTTGCATCTGCTCGatat                                 >  2:597277/1‑149 (MQ=255)
                                                                                                                        ggAGCCATCTCCAGCAGGATCATGCCGGCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTaa                             <  1:1151186/149‑1 (MQ=255)
                                                                                                                        ggAGCCATCTCAAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCCTCAGCTCGATATTTaa                             <  1:1027997/149‑1 (MQ=255)
                                                                                                                         gAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAAt                            <  1:721127/149‑1 (MQ=255)
                                                                                                                           gCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATga                          <  1:589120/149‑1 (MQ=255)
                                                                                                                            ccATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATgaa                         >  1:1371300/1‑148 (MQ=255)
                                                                                                                            ccATCTCCAGCAGGAACATGCCGGCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTGCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATgag                         <  2:872663/149‑1 (MQ=255)
                                                                                                                                 tcCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTa                    <  2:476213/149‑1 (MQ=255)
                                                                                                                                  ccAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACTGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTGGCCCCGCCTCCAGACTTTATTGGATCAGGTCTAGATTTAATGAGaagaag                   >  2:1128586/1‑149 (MQ=255)
                                                                                                                                   cagcagGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGc                  >  2:1538496/1‑149 (MQ=255)
                                                                                                                                    agcagGATCATGCCGTCCGGCATACACTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAACGCCTTGCCCCGCCTCCAGACTTTCTTGCATCCGCTCGATATTTAATGAGAAGTAGCt                 >  1:466516/1‑149 (MQ=255)
                                                                                                                                     gcagGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTc                <  1:1269420/149‑1 (MQ=255)
                                                                                                                                      ccagcgtaTGCCGGACGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGCCTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCa               >  2:676360/8‑149 (MQ=255)
                                                                                                                                        agcgtaTGCCGGACGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCCGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGCGTTTCTTGCTTCAGCGCGATATTTAATGAGAAGTAGCTCAAc             >  1:322057/6‑149 (MQ=255)
                                                                                                                                        agcgtaTGCCGGACGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAAc             >  1:1162669/6‑149 (MQ=255)
                                                                                                                                         gATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACa            >  1:8296/1‑149 (MQ=255)
                                                                                                                                             aTGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGCGCGCCCGGCGATGCCCAGCGTCACTTCGTGATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTc        >  2:628603/1‑149 (MQ=255)
                                                                                                                                             aTGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTc        <  2:825817/149‑1 (MQ=255)
                                                                                                                                             aTGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTc        <  1:1188267/149‑1 (MQ=255)
                                                                                                                                              tGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTcc       <  2:763378/149‑1 (MQ=255)
                                                                                                                                               gCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCg      >  1:802028/1‑149 (MQ=255)
                                                                                                                                                ccGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCgg     <  2:4621/149‑1 (MQ=255)
                                                                                                                                                 cGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTTATGAGAAGTAGCTCAACAGTTCCGGt    <  2:1519829/149‑1 (MQ=255)
                                                                                                                                                   tCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTaa  >  1:1350452/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCACGTAAACCGCCAAGCGGATTTTTAATCTCATGTGCCAGGCCGCGCACTAAATCACGGGCAGCAACCTGCTGGGCGTGCTGTAGCTGTTCCTGACTTAAGCGGCGCTGGTTATCCATCGGAGCCATCTCCAGCAGGATCATGCCGTCCGGCATACGCTGGGCCGTCACAGAAAGGATATGCGAGCGCCCGTCGATGACCAGCGTCACTTCGTTATCGGTAAAACCTTGCCCCGCCTCCAGACTTTCTTGCATCAGCTCGATATTTAATGAGAAGTAGCTCAACAGTTCCGGTAA  >  NC_000913/4055891‑4056186

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: ‑