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‑‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCA  >  NC_000913/1665091‑1665286
                                                                                                  |                                                                                                     
ccTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑t                                                                                                   >  2:1709806/1‑100 (MQ=255)
   gTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tt‑tt                                                                                               >  1:1102726/1‑100 (MQ=255)
   gTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tt‑tt                                                                                               >  2:598502/1‑100 (MQ=255)
     cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tt‑tttt                                                                                             <  2:989928/100‑1 (MQ=255)
     cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tt‑tttt                                                                                             >  1:1376641/1‑100 (MQ=255)
        ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCCCTAATAGGTTACAGGGTCACG‑‑t‑tt‑ttt                                                                                              >  2:96087/1‑96 (MQ=255)
        ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATCGGTTACAGTGTCACG‑‑t‑tt‑tttttat                                                                                          >  1:808459/1‑97 (MQ=255)
         gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tt‑tttttatc                                                                                         <  2:784280/100‑5 (MQ=255)
         gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tt‑tttttatc                                                                                         <  2:1809423/100‑5 (MQ=255)
         gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tt‑tttttatc                                                                                         <  1:1556703/100‑5 (MQ=255)
          ttAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑t‑tt‑tttttttct                                                                                        >  2:332253/1‑95 (MQ=255)
             aaCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTATCACG‑‑TCTT‑TTTTATCTCtt                                                                                      >  2:2110710/1‑100 (MQ=255)
               ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTaa                                                                                    <  2:142383/100‑1 (MQ=255)
                cgatctttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTaaa                                                                                   <  1:1059391/95‑1 (MQ=255)
                cGCTGTTTCGATAGGATCACCGTAATATTGCCGGATTGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTaaa                                                                                   <  2:597935/100‑1 (MQ=255)
                 gCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGAGGATTACGCGCAATCCGGCAATAAGAGGTTAGAGTGTCACGT‑T‑TT‑TTTTATCTGTTaa                                                                                    >  2:1281254/1‑98 (MQ=255)
                 gCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAg                                                                                  >  1:1440994/1‑100 (MQ=255)
                 gCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAg                                                                                  <  1:875496/100‑1 (MQ=255)
                 gCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAg                                                                                  <  1:623991/100‑1 (MQ=255)
                 gCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑T‑TT‑TTTTATCTCTTAAAgc                                                                                 <  1:130143/100‑1 (MQ=255)
                   tGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAgca                                                                                <  1:2011907/100‑1 (MQ=255)
                    gTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAgcac                                                                               >  1:1387699/1‑100 (MQ=255)
                    gTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAgcac                                                                               >  1:393087/1‑100 (MQ=255)
                    gTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAgca                                                                                >  1:2123622/1‑99 (MQ=255)
                      ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTaaaagatcg                                                                             >  2:1059383/1‑94 (MQ=255)
                      ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAgcacgc                                                                             <  1:2273842/100‑1 (MQ=255)
                      ccgatctggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAgcacgc                                                                             <  1:91648/93‑1 (MQ=255)
                            aGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCt                                                                       >  1:1002945/1‑100 (MQ=255)
                             ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAgcacgcagatcgg                                                                      >  2:91647/1‑94 (MQ=255)
                                  aCCGTCCTATTGGCGGATGGTGAGGATTACGCTCTATCTGGCCATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGc                                                                  <  2:1506207/99‑1 (MQ=255)
                                    cGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGc                                                               <  2:1734506/100‑1 (MQ=255)
                                         tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTAGcc                                                          >  2:1893744/1‑100 (MQ=255)
                                          atTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTT‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc                                                          >  1:2339995/1‑100 (MQ=255)
                                           ttGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTAGCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctc                                                        <  1:981737/100‑1 (MQ=255)
                                            tGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCctct                                                       >  1:621287/1‑100 (MQ=255)
                                             gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtt                                                      >  2:1687249/1‑100 (MQ=255)
                                              ccGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCttt                                                     <  2:258287/100‑1 (MQ=255)
                                               cGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGTCCTCTTTTg                                                   >  2:1017591/1‑100 (MQ=255)
                                                ggATGGTGACGATTCCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                                   <  1:859483/100‑1 (MQ=255)
                                                ggATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                                   <  2:107188/100‑1 (MQ=255)
                                                ggATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                                   >  1:792261/1‑100 (MQ=255)
                                                ggATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                                   <  2:1163540/100‑1 (MQ=255)
                                                ggATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                                   >  1:2209931/1‑100 (MQ=255)
                                                   tGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCg                                                <  1:147591/100‑1 (MQ=255)
                                                    ggTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGc                                               <  1:1545517/100‑1 (MQ=255)
                                                    ggTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGc                                               >  1:1389632/1‑100 (MQ=255)
                                                      tGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaa                                             >  1:131405/1‑100 (MQ=255)
                                                       gACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaaa                                            >  1:1072835/1‑100 (MQ=255)
                                                         cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                          >  1:1142291/1‑100 (MQ=255)
                                                           aTTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAg                                        <  1:1775922/100‑1 (MQ=255)
                                                           aTTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAg                                        >  1:369787/1‑100 (MQ=255)
                                                            ttACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGt                                       >  1:414631/1‑100 (MQ=255)
                                                              aCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCg                                     >  1:1202259/1‑100 (MQ=255)
                                                                gcgcAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGcc                                   <  2:1892769/100‑1 (MQ=255)
                                                                gcgcAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGcc                                   <  2:1326733/100‑1 (MQ=255)
                                                                   cAATCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgt                                >  2:1428876/1‑100 (MQ=255)
                                                                    aaTCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgtg                               <  1:1428887/100‑1 (MQ=255)
                                                                     aTCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgtgt                              >  2:2342702/1‑100 (MQ=255)
                                                                      tCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGTACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGtt                             >  2:1944979/1‑100 (MQ=255)
                                                                      tCCGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGtt                             <  1:1368722/100‑1 (MQ=255)
                                                                       ccGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGttt                            <  2:2479330/100‑1 (MQ=255)
                                                                       ccGGCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGttt                            <  1:1324709/100‑1 (MQ=255)
                                                                         ggCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCa                          <  1:138348/100‑1 (MQ=255)
                                                                         ggCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCa                          >  2:1776791/1‑100 (MQ=255)
                                                                         ggCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCa                          <  2:1516212/100‑1 (MQ=255)
                                                                          gCAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCCCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAt                         >  1:1833984/1‑100 (MQ=255)
                                                                           cAATAATAGGTTACAGTGTCACGT‑T‑TT‑TTTTATTTCTTATTTCATTCTCTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAtt                        >  1:1750718/1‑100 (MQ=255)
                                                                            aataatAGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTg                       <  1:598509/100‑1 (MQ=255)
                                                                              taataGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                     <  2:583978/100‑1 (MQ=255)
                                                                               aataGGTTACAGTATCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc                    <  1:621789/100‑1 (MQ=255)
                                                                                ataGGTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCa                   <  1:1607481/100‑1 (MQ=255)
                                                                                    gTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc               >  1:1853385/1‑100 (MQ=255)
                                                                                    gTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc               >  1:1552234/1‑100 (MQ=255)
                                                                                    gTTACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc               >  1:908924/1‑100 (MQ=255)
                                                                                      tACAGTGTCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTg             >  1:2071952/1‑100 (MQ=255)
                                                                                       actgtgtCCCGT‑T‑TT‑TTTTATCTCTCAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc            <  1:2163334/97‑1 (MQ=255)
                                                                                          gtgtCGCGT‑T‑TT‑TTTTATCTCTTAGAGCCCGCACCGCTTTTGCGCCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTG‑CCATTTCTGCTc          <  2:1925571/98‑1 (MQ=255)
                                                                                          gtgtCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCa         <  1:1776804/100‑1 (MQ=255)
                                                                                          gtgtCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCa         <  1:181884/100‑1 (MQ=255)
                                                                                           tgtCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCGCTTTTGCCGCAAAATAGTCGCACGTGTTTCATTGCCCATTTCTGCTCAt        <  2:156725/100‑1 (MQ=255)
                                                                                           tgtCACG‑‑T‑TTATGTTATCTCTTCAAGCACGCACTGCTTTTGCGGCTGGCCTCTTGTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc            <  1:857838/96‑1 (MQ=255)
                                                                                            gtCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg       <  1:948407/100‑1 (MQ=255)
                                                                                            gtCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg       >  2:2074512/1‑100 (MQ=255)
                                                                                            gtCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg       <  1:1747252/100‑1 (MQ=255)
                                                                                             tCACGT‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGc      >  1:2425681/1‑100 (MQ=255)
                                                                                                cG‑‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGcatca  <  2:2017138/100‑1 (MQ=255)
                                                                                                  |                                                                                                     
CCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACG‑‑T‑TT‑TTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATGCATCA  >  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: ‑