Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,915,105 C→T P17S (CCG→TCG)  barA → hybrid sensory histidine kinase, in two‑component regulatory system with UvrY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,915,1050CT98.5% 206.2 / ‑5.3 65P17S (CCG→TCG) barAhybrid sensory histidine kinase, in two‑component regulatory system with UvrY
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/34);  minor base G (0/1);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCACCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGCCGGTGCCAGCATTATTGAGCCGCTTGCAGTTTCTACTGAATATGGCA  >  NC_000913/2914970‑2915234
                                                                                                                                       |                                                                                                                                 
tGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCatcga                                                                                                                                <  2:1530824/139‑1 (MQ=255)
  cgtcgtTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCatcgacc                                                                                                                              <  1:457439/139‑1 (MQ=255)
  cgtcgtTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAAGTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCatcgacc                                                                                                                              <  1:451567/139‑1 (MQ=255)
       ttttGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCt                                                                                                                         >  2:1286461/1‑139 (MQ=255)
        tttGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCtt                                                                                                                        <  2:78606/139‑1 (MQ=255)
                  aGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTAt                                                                                                              <  1:761544/139‑1 (MQ=255)
                      ttGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCatcgac                                                                                                                               >  2:589025/1‑118 (MQ=255)
                      ttGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCatcgac                                                                                                                               <  1:589025/118‑1 (MQ=255)
                      ttGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCt                                                                                                          <  1:1610008/139‑1 (MQ=255)
                       tGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTg                                                                                                         <  1:1286461/139‑1 (MQ=255)
                           ccATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTa                                                                                                     >  2:244742/1‑139 (MQ=255)
                                   gagaCTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCatcgac                                                                                                                               <  2:1809560/105‑1 (MQ=255)
                                   gagaCTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCatcgac                                                                                                                               >  1:1809560/1‑105 (MQ=255)
                                         ttCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTg                                                                                       <  2:798241/139‑1 (MQ=255)
                                         ttCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTg                                                                                       <  2:1893365/139‑1 (MQ=255)
                                           ctcAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTg                                                                                       >  1:556992/1‑137 (MQ=255)
                                           ctcAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTg                                                                                       <  2:556992/137‑1 (MQ=255)
                                           ctcAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCa                                                                                     <  1:244742/139‑1 (MQ=255)
                                            tcAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCAt                                                                                    <  2:1906109/139‑1 (MQ=255)
                                            tcAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCAt                                                                                    >  2:1410831/1‑139 (MQ=255)
                                                tttAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCatcga                                                                                                                                >  1:824328/1‑91 (MQ=255)
                                                tttAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCatcga                                                                                                                                <  2:824328/91‑1 (MQ=255)
                                                        tgtgACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTg                                                                                       >  2:1455476/1‑124 (MQ=255)
                                                        tgtgACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGAc                                                                        >  1:1552761/1‑139 (MQ=255)
                                                        tgtgACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCAGCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTg                                                                                       <  1:1455476/124‑1 (MQ=255)
                                                          tgACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACtt                                                                      >  2:1898679/1‑139 (MQ=255)
                                                           gACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTg                                                                     <  1:1898679/139‑1 (MQ=255)
                                                             ccTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCatcga                                                                                                                                <  2:912012/78‑1 (MQ=255)
                                                             ccTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCatcga                                                                                                                                >  1:912012/1‑78 (MQ=255)
                                                                tAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCg                                                                <  1:1252059/139‑1 (MQ=255)
                                                                tAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCg                                                                <  2:1552761/139‑1 (MQ=255)
                                                                      tCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACt                                                          >  1:881325/1‑139 (MQ=255)
                                                                       cccATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCatcgacc                                                                                                                              >  1:1320668/1‑70 (MQ=255)
                                                                       cccATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCatcgacc                                                                                                                              <  2:1320668/70‑1 (MQ=255)
                                                                        ccATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCatcga                                                                                                                                <  2:1887329/67‑1 (MQ=255)
                                                                        ccATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCatcga                                                                                                                                >  1:1887329/1‑67 (MQ=255)
                                                                           tAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAg                                                     >  1:1286761/1‑139 (MQ=255)
                                                                               ggAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGc                                                 <  1:569981/139‑1 (MQ=255)
                                                                                gAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGGATCGACCGTCCTTAt                                                                                                                      <  1:1374924/69‑1 (MQ=255)
                                                                                gAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTAt                                                                                                                      >  2:1374924/1‑69 (MQ=255)
                                                                                gAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGcc                                                >  2:772961/1‑139 (MQ=255)
                                                                                gAACTCCATGACCAACTACAGCCTGCGCGCACGCATGAAGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGcc                                                <  2:881325/139‑1 (MQ=255)
                                                                                   cTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCg                                                                                  <  2:1588779/102‑1 (MQ=255)
                                                                                   cTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCg                                                                                  >  1:1588779/1‑102 (MQ=255)
                                                                                     ccATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCACCTGGAAGATGCCGGTGc                                           <  2:1286761/139‑1 (MQ=255)
                                                                                        tGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTg                                                                                       <  2:1981154/92‑1 (MQ=255)
                                                                                        tGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTg                                                                                       >  1:1981154/1‑92 (MQ=255)
                                                                                        tGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGc                                                                                      <  2:1792693/93‑1 (MQ=255)
                                                                                        tGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGc                                                                                      >  1:1792693/1‑93 (MQ=255)
                                                                                                     ccTGCGCGCACGCATGATGATTCTGATCCTGGCatcga                                                                                                                                <  2:1035936/38‑1 (MQ=255)
                                                                                                     ccTGCGCGCACGCATGATGATTCTGATCCTGGCatcga                                                                                                                                >  1:1035936/1‑38 (MQ=255)
                                                                                                       tGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCAt                                                                                    <  1:948404/80‑1 (MQ=255)
                                                                                                       tGCGCGCACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCAt                                                                                    >  2:948404/1‑80 (MQ=255)
                                                                                                         cgcgcACTCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCt                                                                                                          <  1:1775787/56‑1 (MQ=255)
                                                                                                         cgcgcACGCATGATGATTCTGATCCTGGCatcgacc                                                                                                                              <  1:580104/36‑1 (MQ=255)
                                                                                                         cgcgcACGCATGATGATTCTGATCCTGGCatcgacc                                                                                                                              <  2:1678874/36‑1 (MQ=255)
                                                                                                         cgcgcACGCATGATGATTCTGATCCTGGCatcgacc                                                                                                                              >  2:580104/1‑36 (MQ=255)
                                                                                                         cgcgcACGCATGATGATTCTGATCCTGGCatcgacc                                                                                                                              >  1:1678874/1‑36 (MQ=255)
                                                                                                         cgcgcACGCATGATGATTCTGATCCTGGCATCGACCGTc                                                                                                                           >  2:761532/1‑39 (MQ=255)
                                                                                                         cgcgcACGCATGATGATTCTGATCCTGGCATCGACCGTc                                                                                                                           <  1:761532/39‑1 (MQ=255)
                                                                                                         cgcgcACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCt                                                                                                          >  2:1775787/1‑56 (MQ=255)
                                                                                                         cgcgcACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCtttt                                                                                               >  1:1062016/1‑67 (MQ=255)
                                                                                                         cgcgcACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCtttt                                                                                               <  2:1062016/67‑1 (MQ=255)
                                                                                                         cgcgcACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCa                                                                   >  2:1497238/1‑95 (MQ=255)
                                                                                                         cgcgcACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCa                                                                   <  1:1497238/95‑1 (MQ=255)
                                                                                                         cgcgcACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGCCGGTGCCAGCATTATTGAg                              >  1:229854/1‑132 (MQ=255)
                                                                                                         cgcgcACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGCCGGTGCCAGCATTATTGAg                              <  2:229854/132‑1 (MQ=255)
                                                                                                         cgcgcACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGCCGGTGCCAGCATTATTGAGCCGCTTg                       >  2:642283/1‑139 (MQ=255)
                                                                                                         cgcgcACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGCCGGTGCCAGCATTATTGAGCCGCTTg                       >  2:1053796/1‑139 (MQ=255)
                                                                                                          gcgcACGCATGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGCCGGTGCCAGCATTATTGAGCCGCTTGc                      <  1:1410831/139‑1 (MQ=255)
                                                                                                               cgcaTGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCt                                                                                                          <  2:1130390/50‑1 (MQ=255)
                                                                                                               cgcaTGATGATTCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCt                                                                                                          >  1:1130390/1‑50 (MQ=255)
                                                                                                                  atgatgATTCTGATCCTGGTATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGCCGGTGCCAGCATTATTGAGCCGCTTGCAGTTTCTa              <  1:1053796/139‑1 (MQ=255)
                                                                                                                          tCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCg                                                                <  2:1711752/81‑1 (MQ=255)
                                                                                                                          tCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCg                                                                >  1:1711752/1‑81 (MQ=255)
                                                                                                                          tCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGCCGGTGCCAGCATTATTGAGCCGCTTGCa                     <  2:1137208/124‑1 (MQ=255)
                                                                                                                          tCTGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGCCGGTGCCAGCATTATTGAGCCGCTTGCa                     >  1:1137208/1‑124 (MQ=255)
                                                                                                                            tGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCAt                                                                                    <  2:796648/59‑1 (MQ=255)
                                                                                                                            tGATCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCAt                                                                                    >  1:796648/1‑59 (MQ=255)
                                                                                                                              aTCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGCCGGTGCCCGCATTATTGAGCCGCTTGCAGTTTCTa              >  1:1805629/1‑127 (MQ=255)
                                                                                                                              aTCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGCCGGTGCCAGCATTATTGAGCCGCTTGCAGTTTCTa              <  2:1805629/127‑1 (MQ=255)
                                                                                                                              aTCCTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGCCGGTGCCAGCATTATTGAGCCGCTTGCAGTTTCTACTGAATATGGCa  >  2:1311304/1‑139 (MQ=255)
                                                                                                                                ccTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGCCGGTGc                                           >  1:982244/1‑96 (MQ=255)
                                                                                                                                ccTGGCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGCCGGTGc                                           <  2:982244/96‑1 (MQ=255)
                                                                                                                                   ggCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGc                                                                                      <  2:1983657/50‑1 (MQ=255)
                                                                                                                                   ggCATCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGc                                                                                      >  1:1983657/1‑50 (MQ=255)
                                                                                                                                       |                                                                                                                                 
TGCGTCGTTTTGCAGAGTAGAATTGCGCCATTGGCGAGACTTTCTCAATTTAACAGTGTGACCTTAATTGTCCCATAACGGAACTCCATGACCAACTACAGCCTGCGCGCACGCATGATGATTCTGATCCTGGCACCGACCGTCCTTATTGGTTTATTGCTGAGTATCTTTTTCGTCGTGCATCGCTATAACGACTTGCAGCGTCAACTGGAAGATGCCGGTGCCAGCATTATTGAGCCGCTTGCAGTTTCTACTGAATATGGCA  >  NC_000913/2914970‑2915234

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