Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,665,195 (T)7→8 intergenic (+75/‑120) dmsD → / → clcB twin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,665,1881.T96.7% 369.8 / 1.4 92intergenic (+68/‑127)dmsD/clcBtwin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base T (35/54);  total (35/57)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTT  >  NC_000913/1665039‑1665332
                                                                                                                                                      |                                                                                                                                                 
gggTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑t                                                                                                                                                  >  1:253855/1‑151 (MQ=255)
              gCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCtt                                                                                                                                    >  1:461618/1‑151 (MQ=255)
                 tAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑T‑TTTTTTATCTCTTAAAg                                                                                                                                <  1:1470876/151‑1 (MQ=255)
                  aaCGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAg                                                                                                                                <  2:1448696/151‑1 (MQ=255)
                   aCGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAgc                                                                                                                               <  2:365844/151‑1 (MQ=255)
                        ggCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAgcacgca                                                                                                                          >  2:840841/1‑151 (MQ=255)
                          cgcAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAgcacgcac                                                                                                                         <  2:206611/150‑1 (MQ=255)
                          cgcAATGGCAATCACAACGGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAgcacgcac                                                                                                                         <  1:1097170/150‑1 (MQ=255)
                             aaTGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCt                                                                                                                     <  2:837860/151‑1 (MQ=255)
                               tGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCttt                                                                                                                   <  2:1447728/151‑1 (MQ=255)
                                tgCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGGTCCGAACACCGTAATATTGCCGGATGGAGACAATTACGTGTAATCCAGCAATCATAGGTTACAGAGCCACGTT‑TTTTTTATCTCTTAAAgactgaactgctctt                                                                                                                  >  1:159102/2‑137 (MQ=255)
                                ggCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCtttt                                                                                                                  >  1:128221/1‑151 (MQ=255)
                                  cAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTg                                                                                                                 <  1:639862/150‑1 (MQ=255)
                                  cAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTg                                                                                                                 <  1:198068/150‑1 (MQ=255)
                                  cAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTg                                                                                                                 <  1:1495570/150‑1 (MQ=255)
                                  cAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGGTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGc                                                                                                                <  1:766651/151‑1 (MQ=255)
                                    aTCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTGCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                                                                              <  1:448538/151‑1 (MQ=255)
                                     tCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                                                                              >  2:55995/1‑150 (MQ=255)
                                        caACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTgg                                                                                                          >  1:933810/1‑151 (MQ=255)
                                        caACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTgg                                                                                                          >  2:1477320/1‑151 (MQ=255)
                                        caACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTgg                                                                                                          >  2:1437640/1‑151 (MQ=255)
                                           cTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGc                                                                                                         >  2:850865/1‑149 (MQ=255)
                                             gTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCct                                                                                                       >  2:1167664/1‑149 (MQ=255)
                                             gTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGACTGGCct                                                                                                       >  2:1330890/1‑149 (MQ=255)
                                               tAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCttt                                                                                                   <  1:1442108/151‑1 (MQ=255)
                                                             gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAg                                                                                      <  2:1121063/150‑1 (MQ=255)
                                                             gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAg                                                                                      <  2:94564/150‑1 (MQ=255)
                                                              ttAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGt                                                                                     <  2:209593/150‑1 (MQ=255)
                                                                      cTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgtgt                                                                            >  2:409027/1‑151 (MQ=255)
                                                                         tttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGttt                                                                          <  2:1440192/150‑1 (MQ=255)
                                                                         tttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGttt                                                                          <  2:461780/150‑1 (MQ=255)
                                                                           tCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAt                                                                       <  2:1012599/151‑1 (MQ=255)
                                                                           tCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAt                                                                       <  2:800054/151‑1 (MQ=255)
                                                                           tCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAt                                                                       <  2:895728/151‑1 (MQ=255)
                                                                            cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAtt                                                                      >  1:430415/1‑151 (MQ=255)
                                                                                 ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc                                                                  >  1:489610/1‑150 (MQ=255)
                                                                                 ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCa                                                                 >  2:1369185/1‑151 (MQ=255)
                                                                                 ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCa                                                                 >  2:847818/1‑151 (MQ=255)
                                                                                       ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt                                                              >  2:446546/1‑148 (MQ=255)
                                                                                       ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt                                                              >  2:782433/1‑148 (MQ=255)
                                                                                       ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt                                                              >  1:399945/1‑148 (MQ=255)
                                                                                       ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt                                                              >  1:226282/1‑148 (MQ=255)
                                                                                       ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt                                                              >  1:154219/1‑148 (MQ=255)
                                                                                       ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt                                                              >  2:1495395/1‑148 (MQ=255)
                                                                                        cGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc                                                          <  1:674025/151‑1 (MQ=255)
                                                                                        cGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc                                                          <  1:944680/151‑1 (MQ=255)
                                                                                        cGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGACTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc                                                          <  2:1017431/151‑1 (MQ=255)
                                                                                          tAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTc                                                        <  2:312511/151‑1 (MQ=255)
                                                                                               ttGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGca                                                   <  2:451289/151‑1 (MQ=255)
                                                                                                         gTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTaca                                          >  2:559223/1‑150 (MQ=255)
                                                                                                          tGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTacac                                         <  2:30078/150‑1 (MQ=255)
                                                                                                          tGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTacac                                         <  1:872972/150‑1 (MQ=255)
                                                                                                            aCGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATc                                      <  2:1506652/151‑1 (MQ=255)
                                                                                                              gATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTa                                    >  1:1362236/1‑151 (MQ=255)
                                                                                                                 tACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATcc                                 <  2:1493218/151‑1 (MQ=255)
                                                                                                                       cAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGAt                             >  2:577729/1‑149 (MQ=255)
                                                                                                                       cAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATc                            <  2:529960/150‑1 (MQ=255)
                                                                                                                       cAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATc                            <  1:548093/150‑1 (MQ=255)
                                                                                                                       cAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATc                            <  2:910634/150‑1 (MQ=255)
                                                                                                                       cAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATc                            <  1:1142298/150‑1 (MQ=255)
                                                                                                                       cAATCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATc                            <  2:223233/150‑1 (MQ=255)
                                                                                                                         aTCCGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTg                          >  1:122493/1‑150 (MQ=255)
                                                                                                                           ccGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTgcgc                       <  1:198318/151‑1 (MQ=255)
                                                                                                                           ccGGCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTgcgc                       <  2:1476646/151‑1 (MQ=255)
                                                                                                                              gCAATAATAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCAc                     <  1:1117224/150‑1 (MQ=255)
                                                                                                                                aataatAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTAt                  <  2:994875/151‑1 (MQ=255)
                                                                                                                                aataatAGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTAt                  <  1:224645/151‑1 (MQ=255)
                                                                                                                                 ataataGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATg                 >  1:1030050/1‑151 (MQ=255)
                                                                                                                                  taataGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGt                <  1:189166/151‑1 (MQ=255)
                                                                                                                                    ataGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTc                                                        <  2:499433/109‑1 (MQ=255)
                                                                                                                                    ataGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTc                                                        >  1:499241/1‑109 (MQ=255)
                                                                                                                                    ataGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGGCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTATACATCTATCCGGATCTGCGCACTATGtt               >  1:729753/1‑150 (MQ=255)
                                                                                                                                     taGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTcc             <  2:335209/151‑1 (MQ=255)
                                                                                                                                     taGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTcc             <  1:116166/151‑1 (MQ=255)
                                                                                                                                     taGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTcc             <  1:1496369/151‑1 (MQ=255)
                                                                                                                                     taGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTcc             <  2:629533/151‑1 (MQ=255)
                                                                                                                                     taGGTTACAGTGTCACG‑‑‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCAc           <  1:1116084/151‑1 (MQ=255)
                                                                                                                                      aGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTTATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCa            >  1:969269/1‑151 (MQ=255)
                                                                                                                                      aGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCa            >  2:1231791/1‑151 (MQ=255)
                                                                                                                                      aGGTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCa            >  2:414549/1‑151 (MQ=255)
                                                                                                                                       ggTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCa            >  2:243553/1‑150 (MQ=255)
                                                                                                                                        gTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCAc           <  1:623876/150‑1 (MQ=255)
                                                                                                                                        gTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCAc           <  1:640276/150‑1 (MQ=255)
                                                                                                                                        gTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCAc           <  1:1396158/150‑1 (MQ=255)
                                                                                                                                        gTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCAc           <  1:1301572/150‑1 (MQ=255)
                                                                                                                                        gTTACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCAc           <  2:430560/150‑1 (MQ=255)
                                                                                                                                         ttACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCg         <  1:540751/151‑1 (MQ=255)
                                                                                                                                         ttACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCg         <  1:497615/151‑1 (MQ=255)
                                                                                                                                         ttACAGTGTCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCg         <  2:843864/151‑1 (MQ=255)
                                                                                                                                              gtgtCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTg     >  1:917826/1‑150 (MQ=255)
                                                                                                                                              gtgtCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTg     >  2:1164553/1‑150 (MQ=255)
                                                                                                                                               tttCACG‑TGTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAGCAGTCGCCCGTGTTGCATTGCCCATTTCTGATCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCGGc    <  1:780790/148‑1 (MQ=255)
                                                                                                                                                 tCACGTT‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCtt  <  1:1302626/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
GGGTGAACTGGCTCGCTTAACGCTGGCGCAATGGCAATCACAACTGTTAATTCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑T‑TTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCATCTACACATCTATCCGGATCTGCGCACTATGTTCCACCGTCTGCTT  >  NC_000913/1665039‑1665332

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑