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.T94.6% 229.5 / 6.3 74intergenic (+68/‑127)dmsD/clcBtwin‑argninine leader‑binding protein for DmsA and TorA/H(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base T (35/34);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATG  >  NC_000913/1665090‑1665281
                                                                                                   |                                                                                               
tCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGt                                                                                                >  2:1632982/1‑100 (MQ=255)
  cTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttt                                                                                              <  2:1441941/100‑1 (MQ=255)
  cTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttt                                                                                              >  1:2469860/1‑100 (MQ=255)
  cTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttt                                                                                              >  1:1363862/1‑100 (MQ=255)
     tCGCGGTTAAACCGCTGTTTCGATAGTATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGtttttt                                                                                         <  2:652803/100‑1 (MQ=255)
      cgcgGTTAAACCGTTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAAGCCGGCAATAATAGGTTACAGTGTCACGttttttt                                                                                        <  1:1328510/100‑1 (MQ=255)
      cgcgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttt                                                                                        <  2:97896/100‑1 (MQ=255)
       gcgGTTAAACCGCTGTTTCGATAGGAGCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGtttttttt                                                                                       <  2:2159509/100‑2 (MQ=255)
        cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGtttttttta                                                                                      <  1:709616/100‑3 (MQ=255)
        cgGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGtttttttta                                                                                      <  2:2552992/100‑3 (MQ=255)
         ggTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttttat                                                                                     <  1:2126284/100‑4 (MQ=255)
          gTTAAACCGCTGTTTCGATAGTATCACTGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGtttttttt                                                                                       <  1:1453526/97‑2 (MQ=255)
          gTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttttatc                                                                                    <  2:2604592/100‑5 (MQ=255)
           ttAAACCGCTTTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttttatct                                                                                   <  2:883856/100‑6 (MQ=255)
           ttAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGttttttttatc                                                                                    <  2:844637/99‑5 (MQ=255)
            tAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAATGTAACGTTTTTCTTAtctc                                                                                   >  2:2433236/1‑100 (MQ=255)
               aCCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTa                                                                                <  2:1279586/100‑1 (MQ=255)
                ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTaa                                                                               >  1:2534434/1‑100 (MQ=255)
                ccGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTaa                                                                               >  1:2090744/1‑100 (MQ=255)
                  gCTGTTTCTATAGGATCACCGTAATATTGTCGGATGGTGACGATTACTCGCAATCGGGCAATAATAGGTTACAGGGGCACGTTTTTTTTATCTCTTaa                                                                               <  2:127021/98‑1 (MQ=255)
                  gCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTaagg                                                                             >  1:1630527/1‑98 (MQ=255)
                  gCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCAGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAg                                                                             >  1:1989984/1‑100 (MQ=255)
                      tttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacg                                                                         <  2:934119/100‑1 (MQ=255)
                       ttCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACATTGTCACGTTTTTTTTATCTCTTAAAgcactc                                                                        >  1:1424289/1‑98 (MQ=255)
                        tCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacgca                                                                       <  2:2817761/100‑1 (MQ=255)
                         cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacgcac                                                                      >  1:2098577/1‑100 (MQ=255)
                         cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacgcac                                                                      <  2:2315953/100‑1 (MQ=255)
                         cGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAgcacgcac                                                                      <  2:1713246/100‑1 (MQ=255)
                              ggATCACCGTAATATTGCCGGATGGTGTCGATTAGGCGCAATTCGGCAATAATAGGTTACAGTGTCACGTATTTTATCTCTTAAAGCACGCACTGCtttt                                                               >  1:225240/1‑100 (MQ=255)
                              ggATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCtt                                                                 >  1:18685/1‑100 (MQ=255)
                               gATCACCGTAATCTTGCCGGATGGTGACGATTACGCGCAATCCGACAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCttt                                                                >  1:2479373/1‑100 (MQ=255)
                                aTCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCttt                                                                <  2:2010527/99‑1 (MQ=255)
                                  cACCGTAATATTGCCGGATTGTGACGATTACGCTCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGc                                                             <  1:2796696/100‑1 (MQ=255)
                                   actgTAATATTGCCGGATGGTGATGATTACGCGCTTTCCGGCAATAATAGGTTACAGTGTCGCGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                            <  1:1195780/97‑1 (MQ=255)
                                   aCCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTAAAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGc                                                             >  2:252704/1‑99 (MQ=255)
                                    ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCgg                                                           >  1:860148/1‑100 (MQ=255)
                                    ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                            >  2:1630098/1‑99 (MQ=255)
                                    ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                            >  1:2030465/1‑99 (MQ=255)
                                    ccGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTATCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCg                                                            >  1:438587/1‑99 (MQ=255)
                                      atAATATTGCCGGATGGTGAAAAATACGAGAAATAAGGCAAAAATAGGTTAAAGTGTAACGTTTCTTTTATATATTAAAGAACGCACTGCTTTTGCggtt                                                         >  2:1166335/2‑98 (MQ=255)
                                       tAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTg                                                        >  2:1236397/1‑100 (MQ=255)
                                        ttccgatctggATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTgg                                                       <  1:52601/91‑1 (MQ=255)
                                        aaTATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTgg                                                       <  1:1185614/100‑1 (MQ=255)
                                          tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc                                                     <  1:2094089/100‑1 (MQ=255)
                                          tatTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGcc                                                     <  1:2031616/100‑1 (MQ=255)
                                           atTGCCGGATGGTGAGGATTATGCTCAATCCGGCAATAATAGGTTACAATGTCACGTTTGTTTTATCTCTTCAAGCACGCACTGCTTTTGCg                                                            <  2:903553/92‑1 (MQ=255)
                                           atTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCct                                                    <  2:523816/100‑1 (MQ=255)
                                              gCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCtt                                                 <  2:1107726/100‑1 (MQ=255)
                                               ccGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCttt                                                <  2:2461391/100‑1 (MQ=255)
                                                 ggATGGTGACGATTACGCGCGATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTg                                              <  2:463475/100‑1 (MQ=255)
                                                 ggATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTggagatcggaa                                              >  2:52601/1‑91 (MQ=255)
                                                  gATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc                                             <  2:1420859/100‑1 (MQ=255)
                                                  gATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc                                             <  1:1415616/100‑1 (MQ=255)
                                                  gATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGc                                             <  1:1236409/100‑1 (MQ=255)
                                                    tGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCTCGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCg                                           <  1:2557023/100‑1 (MQ=255)
                                                    tGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCATGCACTGCTTTTGCGGCTTGCCTCTTTTGcag                                           <  1:2043599/100‑3 (MQ=255)
                                                     ggTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGc                                          >  2:869964/1‑100 (MQ=255)
                                                      gTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCa                                         >  1:2716531/1‑100 (MQ=255)
                                                       tGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCaa                                        >  1:437049/1‑100 (MQ=255)
                                                          cGATTACGCGCAATCCGGCAATACTAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAACACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                     >  1:842407/1‑100 (MQ=255)
                                                          cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTTCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                     >  2:845904/1‑100 (MQ=255)
                                                          cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                     <  2:2548716/100‑1 (MQ=255)
                                                          cGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAAt                                     <  1:1051697/100‑1 (MQ=255)
                                                           gATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATa                                    >  1:2335215/1‑100 (MQ=255)
                                                               aCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTc                                 >  1:1874308/1‑99 (MQ=255)
                                                               aCGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATATCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCTCAAACTTGTc                                 >  2:852388/1‑99 (MQ=255)
                                                                cgcgcgCTCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGc                               <  1:1195947/100‑1 (MQ=255)
                                                                cgcgCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGc                               >  1:1341497/1‑100 (MQ=255)
                                                                  cgcAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGccc                             >  1:2267414/1‑100 (MQ=255)
                                                                     aaTCCGGCAATAATAGGTAACAGAGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCAGCTGCACTCATTTGACGCAAAACAGTCcccc                             >  2:2551510/1‑97 (MQ=255)
                                                                      aTCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCgtgt                         <  1:2487036/100‑1 (MQ=255)
                                                                           gCAATAATAGGTTCCAGCGTCACGATTTTTTTATCTCTTAAAGCAAGCACTGCTTTTACGGCTTGCTTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCAt                    >  2:2112817/1‑100 (MQ=255)
                                                                             aataatAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTg                  <  2:519994/100‑1 (MQ=255)
                                                                             aataatAGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTACTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTg                  <  1:526598/100‑1 (MQ=255)
                                                                              ataataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGc                 >  1:828322/1‑100 (MQ=255)
                                                                              ataataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGc                 <  2:2879113/100‑1 (MQ=255)
                                                                               taataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                >  1:2967425/1‑100 (MQ=255)
                                                                                aataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGccc               <  1:1264695/100‑1 (MQ=255)
                                                                                 ataGGTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTAGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGcc                >  2:1338483/1‑98 (MQ=255)
                                                                                    ggTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTAGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCAttt           >  1:2742212/1‑100 (MQ=255)
                                                                                     gTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTc          <  2:1363850/100‑1 (MQ=255)
                                                                                     gTTACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGACTGAACACTATTGCCACAAAAAAATCACACGTGTATCATTGCCCATTTc          >  1:2534506/1‑100 (MQ=255)
                                                                                      ttACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCt         <  2:548804/100‑1 (MQ=255)
                                                                                       tACAGTGTCACGTTTTTTTTTTTTTTTTAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCt         >  2:884370/1‑99 (MQ=255)
                                                                                       tACAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTg        <  2:424619/100‑1 (MQ=255)
                                                                                        aCAGTGTCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGc       >  1:2921186/1‑100 (MQ=255)
                                                                                           gtgtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCa    >  1:763759/1‑100 (MQ=255)
                                                                                             gtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTc     <  2:2469831/97‑1 (MQ=255)
                                                                                             gtCACGTTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATg  >  1:1692834/1‑100 (MQ=255)
                                                                                                   |                                                                                               
TCCTGTCGCGGTTAAACCGCTGTTTCGATAGGATCACCGTAATATTGCCGGATGGTGACGATTACGCGCAATCCGGCAATAATAGGTTACAGTGTCACGTTTTTTTATCTCTTAAAGCACGCACTGCTTTTGCGGCTGGCCTCTTTTGCCGCAAAATAGTCGCCCGTGTTTCATTGCCCATTTCTGCTCATG  >  NC_000913/1665090‑1665281

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