Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,665,195:1 (T)7→8 100% 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.T93.3% 238.0 / 10.1 75intergenic (+68/‑127)dmsD/clcBtwin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base T (33/37);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCA  >  NC_000913/1665091‑1665286
                                                                                                  |                                                                                                     
ccTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑tt                                                                                                   >  2:1709806/1‑100 (MQ=255)
   gTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttt                                                                                               >  1:1102726/1‑100 (MQ=255)
   gTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttt                                                                                               >  2:598502/1‑100 (MQ=255)
     cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttt                                                                                             <  2:989928/100‑1 (MQ=255)
     cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttt                                                                                             >  1:1376641/1‑100 (MQ=255)
        ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCCCTAATAGGTTACAGGGTCACG‑‑tttttt                                                                                              >  2:96087/1‑96 (MQ=255)
        ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATCGGTTACAGTGTCACG‑‑ttttttttat                                                                                          >  1:808459/1‑97 (MQ=255)
         gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttttatc                                                                                         <  2:784280/100‑5 (MQ=255)
         gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttttatc                                                                                         <  2:1809423/100‑5 (MQ=255)
         gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttttatc                                                                                         <  1:1556703/100‑5 (MQ=255)
          ttAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑ttttttttttct                                                                                        >  2:332253/1‑95 (MQ=255)
             aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTATCACG‑‑TCTTTTTTATCTCtt                                                                                      >  2:2110710/1‑100 (MQ=255)
               ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTaa                                                                                    <  2:142383/100‑1 (MQ=255)
                cgatctttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTaaa                                                                                   <  1:1059391/95‑1 (MQ=255)
                cGCTGTTTCGATAGGATCACCGTAATATTGCCGGATTGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTaaa                                                                                   <  2:597935/100‑1 (MQ=255)
                 gCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGAGGATTACGCGCAATCCGGCAATAAGAGGTTAGAGTGTCACGTTTTTTTTATCTGTTaa                                                                                    >  2:1281254/1‑98 (MQ=255)
                 gCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAg                                                                                  >  1:1440994/1‑100 (MQ=255)
                 gCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAg                                                                                  <  1:875496/100‑1 (MQ=255)
                 gCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAg                                                                                  <  1:623991/100‑1 (MQ=255)
                 gCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑TTTTTTTATCTCTTAAAgc                                                                                 <  1:130143/100‑1 (MQ=255)
                   tGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgca                                                                                <  1:2011907/100‑1 (MQ=255)
                    gTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcac                                                                               >  1:1387699/1‑100 (MQ=255)
                    gTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcac                                                                               >  1:393087/1‑100 (MQ=255)
                    gTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgca                                                                                >  1:2123622/1‑99 (MQ=255)
                      ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTaaaagatcg                                                                             >  2:1059383/1‑94 (MQ=255)
                      ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacgc                                                                             <  1:2273842/100‑1 (MQ=255)
                      ccgatctggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacgc                                                                             <  1:91648/93‑1 (MQ=255)
                            aGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCt                                                                       >  1:1002945/1‑100 (MQ=255)
                             ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacgcagatcgg                                                                      >  2:91647/1‑94 (MQ=255)
                                  aCCGTCCTATTGGCGGATGGTGAGGATTACGCTCTATCTGGCCATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGc                                                                  <  2:1506207/99‑1 (MQ=255)
                                    cGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGc                                                               <  2:1734506/100‑1 (MQ=255)
                                         tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTAGcc                                                          >  2:1893744/1‑100 (MQ=255)
                                          atTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc                                                          >  1:2339995/1‑100 (MQ=255)
                                           ttGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTAGCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctc                                                        <  1:981737/100‑1 (MQ=255)
                                            tGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctct                                                       >  1:621287/1‑100 (MQ=255)
                                             gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtt                                                      >  2:1687249/1‑100 (MQ=255)
                                              ccGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCttt                                                     <  2:258287/100‑1 (MQ=255)
                                               cGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGTCCTCTTTTg                                                   >  2:1017591/1‑100 (MQ=255)
                                                ggATGGTGACGATTCCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                                   <  1:859483/100‑1 (MQ=255)
                                                ggATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                                   <  2:107188/100‑1 (MQ=255)
                                                ggATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                                   >  1:792261/1‑100 (MQ=255)
                                                ggATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                                   <  2:1163540/100‑1 (MQ=255)
                                                ggATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                                   >  1:2209931/1‑100 (MQ=255)
                                                   tGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCg                                                <  1:147591/100‑1 (MQ=255)
                                                    ggTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGc                                               <  1:1545517/100‑1 (MQ=255)
                                                    ggTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGc                                               >  1:1389632/1‑100 (MQ=255)
                                                      tGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaa                                             >  1:131405/1‑100 (MQ=255)
                                                       gACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaaa                                            >  1:1072835/1‑100 (MQ=255)
                                                         cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                          >  1:1142291/1‑100 (MQ=255)
                                                           aTTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAg                                        <  1:1775922/100‑1 (MQ=255)
                                                           aTTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAg                                        >  1:369787/1‑100 (MQ=255)
                                                            ttACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGt                                       >  1:414631/1‑100 (MQ=255)
                                                              aCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCg                                     >  1:1202259/1‑100 (MQ=255)
                                                                gcgcAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGcc                                   <  2:1892769/100‑1 (MQ=255)
                                                                gcgcAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGcc                                   <  2:1326733/100‑1 (MQ=255)
                                                                   cAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgt                                >  2:1428876/1‑100 (MQ=255)
                                                                    aaTCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgtg                               <  1:1428887/100‑1 (MQ=255)
                                                                     aTCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgtgt                              >  2:2342702/1‑100 (MQ=255)
                                                                      tCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGTACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGtt                             >  2:1944979/1‑100 (MQ=255)
                                                                      tCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGtt                             <  1:1368722/100‑1 (MQ=255)
                                                                       ccGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGttt                            <  2:2479330/100‑1 (MQ=255)
                                                                       ccGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGttt                            <  1:1324709/100‑1 (MQ=255)
                                                                         ggCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCa                          <  1:138348/100‑1 (MQ=255)
                                                                         ggCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCa                          >  2:1776791/1‑100 (MQ=255)
                                                                         ggCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCa                          <  2:1516212/100‑1 (MQ=255)
                                                                          gCAATAATAGGTTACAGTGTCACGTTTTTTTTATCCCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAt                         >  1:1833984/1‑100 (MQ=255)
                                                                           cAATAATAGGTTACAGTGTCACGTTTTTTTTATTTCTTATTTCATTCTCTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAtt                        >  1:1750718/1‑100 (MQ=255)
                                                                            aataatAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTg                       <  1:598509/100‑1 (MQ=255)
                                                                              taataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                     <  2:583978/100‑1 (MQ=255)
                                                                               aataGGTTACAGTATCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc                    <  1:621789/100‑1 (MQ=255)
                                                                                ataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCa                   <  1:1607481/100‑1 (MQ=255)
                                                                                    gTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc               >  1:1853385/1‑100 (MQ=255)
                                                                                    gTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc               >  1:1552234/1‑100 (MQ=255)
                                                                                    gTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc               >  1:908924/1‑100 (MQ=255)
                                                                                      tACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTg             >  1:2071952/1‑100 (MQ=255)
                                                                                       actgtgtCCCGTTTTTTTTATCTCTCAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc            <  1:2163334/97‑1 (MQ=255)
                                                                                          gtgtCGCGTTTTTTTTATCTCTTAGAGCCCGCACCGCTTTTGCGCCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCATTTCTGCTc          <  2:1925571/98‑1 (MQ=255)
                                                                                          gtgtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCa         <  1:1776804/100‑1 (MQ=255)
                                                                                          gtgtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCa         <  1:181884/100‑1 (MQ=255)
                                                                                           tgtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCGCTTTTGCCGCAAAATAGTCGCACGTGTTTCATTGCCCATTTCTGCTCAt        <  2:156725/100‑1 (MQ=255)
                                                                                           tgtCACG‑‑TTTATGTTATCTCTTCAAGCACGCACTGCTTTTGCGGCTGGCCTCTTGTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc            <  1:857838/96‑1 (MQ=255)
                                                                                            gtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg       <  1:948407/100‑1 (MQ=255)
                                                                                            gtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg       >  2:2074512/1‑100 (MQ=255)
                                                                                            gtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg       <  1:1747252/100‑1 (MQ=255)
                                                                                             tCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGc      >  1:2425681/1‑100 (MQ=255)
                                                                                                cG‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGcatca  <  2:2017138/100‑1 (MQ=255)
                                                                                                  |                                                                                                     
CCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑TTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCA  >  NC_000913/1665091‑1665286

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