Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,648,144 (T)7→6 100% intergenic (‑356/‑384) dinQ ← / → arsR UV‑inducible membrane toxin, DinQ‑AgrB type I toxin‑antitoxin system/arsenical resistance operon transcriptional repressor; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,648,1380T.98.9% 380.1 / ‑2.4 87intergenic (‑350/‑390)dinQ/arsRUV‑inducible membrane toxin, DinQ‑AgrB type I toxin‑antitoxin system/arsenical resistance operon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (49/37);  total (49/38)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

GTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCC  >  NC_000913/3648041‑3648229
                                                                                                  |                                                                                           
gTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGttt                                                                                           <  1:825092/100‑1 (MQ=255)
gTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGttt                                                                                           >  1:1225481/1‑100 (MQ=255)
gTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGttt                                                                                           <  1:1931207/100‑1 (MQ=255)
   tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGtttttt                                                                                        <  1:748859/100‑1 (MQ=255)
   tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGtttttt                                                                                        <  2:1225477/100‑1 (MQ=255)
   tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGtttttt                                                                                        <  1:1572409/100‑1 (MQ=255)
      taatAGCATTGATGCGTTAGATGAGCGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAACCGACTTTGCGTCGGttttttac                                                                                      <  2:33631/97‑3 (MQ=17)
      tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACTTTGCGTCGGttttttact                                                                                     <  1:108998/100‑4 (MQ=21)
         tAGCATTGATGCGTTAGATGATGGCTAACTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCCAAGTCGG‑TTTTTTACATcc                                                                                 >  2:1901490/1‑100 (MQ=18)
          aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCg                                                                                <  1:941162/100‑1 (MQ=255)
           ncATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCgg                                                                               >  1:1241696/2‑100 (MQ=255)
           gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCg                                                                                <  1:1936745/99‑1 (MQ=255)
            cATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACAGACGCACAGTCGG‑TTGTTTACATCCGGa                                                                              >  1:396472/1‑100 (MQ=21)
                gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCg                                                                          >  1:1283538/1‑100 (MQ=255)
                 aTGCGTTAGATGATGGCTATCTCACTCCATTCATAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCgg                                                                         <  2:1850132/100‑1 (MQ=21)
                  tGCGTTAGATTATGGCTATCTCACTCCAGTCAGAGCCACAAACTCAGTGCTGGAAAGT‑AAAAAACCGACGCATAGCTGG‑TTTTTTACATCCGGATTCg                                                                          <  2:1299188/98‑1 (MQ=18)
                  tGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGa                                                                        <  2:636854/100‑1 (MQ=255)
                     gTTATATCATGGCTATCTCACTCCAGTCAGACCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACaa                                                                     <  2:445565/100‑1 (MQ=18)
                     gTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACaa                                                                     >  1:1502911/1‑100 (MQ=255)
                     gTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACaa                                                                     <  2:2008342/100‑1 (MQ=255)
                       tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCTGATTCGGACAAg                                                                    >  1:1232118/1‑99 (MQ=255)
                       tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAgg                                                                   >  1:1624561/1‑100 (MQ=255)
                       tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAgg                                                                   >  1:1201735/1‑100 (MQ=255)
                       tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAgg                                                                   <  2:1690420/100‑1 (MQ=255)
                       tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATACGGATTCGGACAAgg                                                                   >  1:1000539/1‑100 (MQ=255)
                        aGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAg                                                                    >  1:1887575/1‑98 (MQ=255)
                        aGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGc                                                                  <  2:1858757/100‑1 (MQ=255)
                        aGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGc                                                                  <  1:1184479/100‑1 (MQ=255)
                        aGATGATAGCTATCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGAcctcc                                                                   <  2:1232113/100‑5 (MQ=14)
                         gatgatGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCt                                                                 >  1:478997/1‑100 (MQ=255)
                              tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAata                                                            >  2:548374/1‑100 (MQ=255)
                                 cTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGa                                                         >  1:620830/1‑100 (MQ=255)
                                  tATTTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                        <  1:1750791/100‑1 (MQ=255)
                                  tATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                        <  1:2019266/100‑1 (MQ=255)
                                    tctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACg                                                       >  2:701601/1‑99 (MQ=255)
                                      tcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                    <  2:678451/100‑1 (MQ=255)
                                       cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGa                                                   >  1:1633287/1‑100 (MQ=255)
                                        aCTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAc                                                  >  2:279573/1‑100 (MQ=255)
                                         cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                 >  2:1831336/1‑100 (MQ=255)
                                         cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAc                                                  <  2:1600934/99‑1 (MQ=255)
                                           ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAccc                                                >  1:518379/1‑99 (MQ=255)
                                             aGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAg                                              <  2:1245004/99‑1 (MQ=255)
                                             aGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGt                                             <  2:1326618/100‑1 (MQ=255)
                                              gTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTg                                            >  1:909283/1‑100 (MQ=255)
                                                  gagCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAACGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTGAATATGACGATGCCCCGTTGaaa                                         >  2:669835/1‑99 (MQ=18)
                                                   agCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATAACCCAGTGAAAGt                                       >  1:360351/1‑100 (MQ=255)
                                                    gCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                      <  2:1633276/100‑1 (MQ=255)
                                                    gCCACCAACGCAGGGCTGGAAAGT‑AAAAATCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaaat                                        <  2:1962587/98‑2 (MQ=25)
                                                     ccaccaACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtat                                     >  1:2246376/1‑100 (MQ=255)
                                                     ccaccaACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtat                                     >  2:1978780/1‑100 (MQ=255)
                                                      caccaACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                    <  2:2211153/100‑1 (MQ=255)
                                                         caACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATaa                                 >  2:303452/1‑100 (MQ=255)
                                                           aCTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAAt                               <  1:2097957/100‑1 (MQ=255)
                                                            cTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATaaa                                >  2:2040295/1‑98 (MQ=255)
                                                            cTCAGAGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATaaa                                <  1:68668/98‑1 (MQ=255)
                                                             tCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                             <  1:583403/100‑1 (MQ=255)
                                                              cAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGt                            >  1:190069/1‑100 (MQ=255)
                                                               aGGGCTGGAAAGT‑AAAAAACAGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGAAGATGACACAGTGAAAGTATATAAATTGt                            >  2:391605/1‑99 (MQ=25)
                                                                gggCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGtga                          >  1:2036737/1‑98 (MQ=255)
                                                                  gCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACTATGACCCAGAAAAAGTATAAAAATCTTCAc                         >  2:1063509/1‑99 (MQ=25)
                                                                  gCTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACt                        >  2:1820261/1‑100 (MQ=255)
                                                                   cTGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTg                       >  1:1648125/1‑100 (MQ=255)
                                                                    tGGAAAGT‑AAAAAACCGACGCAAAGTCTG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCTTGTTAAAGTATATAAATCGTCACTga                      <  2:1139333/100‑2 (MQ=37)
                                                                    tGGAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGc                      <  1:532920/100‑1 (MQ=255)
                                                                      gAAAGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGa                    >  2:805883/1‑100 (MQ=255)
                                                                       aaaGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                   <  2:1614059/100‑1 (MQ=255)
                                                                       aaaGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                   <  1:1414672/100‑1 (MQ=255)
                                                                       aaaGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                   <  1:697143/100‑1 (MQ=255)
                                                                       aaaGT‑AAAAAACCGACGAAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                   <  1:353815/100‑1 (MQ=255)
                                                                        aaGT‑AAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGata                  <  1:805883/100‑1 (MQ=255)
                                                                             aaaaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatatat               >  1:2178558/1‑99 (MQ=255)
                                                                              aaaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAc             >  1:90288/1‑100 (MQ=255)
                                                                               aaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc            >  1:2246/1‑100 (MQ=255)
                                                                                  aCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGaa         >  2:1680344/1‑100 (MQ=255)
                                                                                  aCCGACGCAAAGTCGG‑ATTTTTAAATCCGGATTCGGACAAGGATTAATATGATGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGa          >  2:315409/1‑99 (MQ=255)
                                                                                   ccGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg        >  2:1692814/1‑100 (MQ=255)
                                                                                   ccGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg        <  1:2098377/100‑1 (MQ=255)
                                                                                   ccGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg        >  1:246603/1‑100 (MQ=255)
                                                                                    cGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg        <  2:1749729/99‑1 (MQ=255)
                                                                                    cGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGt       <  2:2102869/100‑1 (MQ=255)
                                                                                    cGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGt       >  2:1559840/1‑100 (MQ=255)
                                                                                     gACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGCGAAAGTATATAAATCGTCACTTCGATATATACCGaa         >  1:995148/1‑97 (MQ=255)
                                                                                      aCGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATACATACCGAAGTGc     >  2:1848695/1‑100 (MQ=255)
                                                                                       cGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGct    >  2:2081726/1‑100 (MQ=255)
                                                                                       cGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGct    >  1:2187705/1‑100 (MQ=255)
                                                                                        gCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc   <  1:629714/100‑1 (MQ=255)
                                                                                        gCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc   <  2:48920/100‑1 (MQ=255)
                                                                                        gCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc   <  1:270185/100‑1 (MQ=255)
                                                                                         cAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcc  >  2:416633/1‑100 (MQ=255)
                                                                                         cAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcc  >  1:391378/1‑100 (MQ=255)
                                                                                         cAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcc  <  1:138023/100‑1 (MQ=255)
                                                                                         cAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcc  >  1:103396/1‑100 (MQ=255)
                                                                                         cAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGct    >  2:1585561/1‑98 (MQ=255)
                                                                                         aaaaGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcc  >  2:455339/2‑100 (MQ=255)
                                                                                           aaGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc   <  1:711770/97‑1 (MQ=255)
                                                                                                  |                                                                                           
GTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGT‑AAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCC  >  NC_000913/3648041‑3648229

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