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

GTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCC  >  NC_000913/3648041‑3648229
                                                                                                  |                                                                                           
gTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttt                                                                                           <  1:825092/100‑1 (MQ=255)
gTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttt                                                                                           >  1:1225481/1‑100 (MQ=255)
gTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttt                                                                                           <  1:1931207/100‑1 (MQ=255)
   tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                        <  1:748859/100‑1 (MQ=255)
   tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                        <  2:1225477/100‑1 (MQ=255)
   tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                        <  1:1572409/100‑1 (MQ=255)
      taatAGCATTGATGCGTTAGATGAGCGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAACCGACTTTGCGTCGGttttttac                                                                                      <  2:33631/97‑3 (MQ=17)
      tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGttttttact                                                                                     <  1:108998/100‑4 (MQ=21)
         tAGCATTGATGCGTTAGATGATGGCTAACTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCCAAGTCGGTTTTTTACATcc                                                                                 >  2:1901490/1‑100 (MQ=18)
          aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCg                                                                                <  1:941162/100‑1 (MQ=255)
           ncATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCgg                                                                               >  1:1241696/2‑100 (MQ=255)
           gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCg                                                                                <  1:1936745/99‑1 (MQ=255)
            cATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACAGACGCACAGTCGGTTGTTTACATCCGGa                                                                              >  1:396472/1‑100 (MQ=21)
                gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCg                                                                          >  1:1283538/1‑100 (MQ=255)
                 aTGCGTTAGATGATGGCTATCTCACTCCATTCATAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCgg                                                                         <  2:1850132/100‑1 (MQ=21)
                  tGCGTTAGATTATGGCTATCTCACTCCAGTCAGAGCCACAAACTCAGTGCTGGAAAGTAAAAAACCGACGCATAGCTGGTTTTTTACATCCGGATTCg                                                                          <  2:1299188/98‑1 (MQ=18)
                  tGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGa                                                                        <  2:636854/100‑1 (MQ=255)
                     gTTATATCATGGCTATCTCACTCCAGTCAGACCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACaa                                                                     <  2:445565/100‑1 (MQ=18)
                     gTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACaa                                                                     >  1:1502911/1‑100 (MQ=255)
                     gTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACaa                                                                     <  2:2008342/100‑1 (MQ=255)
                       tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCTGATTCGGACAAg                                                                    >  1:1232118/1‑99 (MQ=255)
                       tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAgg                                                                   >  1:1624561/1‑100 (MQ=255)
                       tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAgg                                                                   >  1:1201735/1‑100 (MQ=255)
                       tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAgg                                                                   <  2:1690420/100‑1 (MQ=255)
                       tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATACGGATTCGGACAAgg                                                                   >  1:1000539/1‑100 (MQ=255)
                        aGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAg                                                                    >  1:1887575/1‑98 (MQ=255)
                        aGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGc                                                                  <  2:1858757/100‑1 (MQ=255)
                        aGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGc                                                                  <  1:1184479/100‑1 (MQ=255)
                        aGATGATAGCTATCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCAGAcctcc                                                                   <  2:1232113/100‑5 (MQ=14)
                         gatgatGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCt                                                                 >  1:478997/1‑100 (MQ=255)
                              tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAata                                                            >  2:548374/1‑100 (MQ=255)
                                 cTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGa                                                         >  1:620830/1‑100 (MQ=255)
                                  tATTTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                        <  1:1750791/100‑1 (MQ=255)
                                  tATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                        <  1:2019266/100‑1 (MQ=255)
                                    tctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACg                                                       >  2:701601/1‑99 (MQ=255)
                                      tcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                    <  2:678451/100‑1 (MQ=255)
                                       cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGa                                                   >  1:1633287/1‑100 (MQ=255)
                                        aCTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAc                                                  >  2:279573/1‑100 (MQ=255)
                                         cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                 >  2:1831336/1‑100 (MQ=255)
                                         cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAc                                                  <  2:1600934/99‑1 (MQ=255)
                                           ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAccc                                                >  1:518379/1‑99 (MQ=255)
                                             aGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAg                                              <  2:1245004/99‑1 (MQ=255)
                                             aGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGt                                             <  2:1326618/100‑1 (MQ=255)
                                              gTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTg                                            >  1:909283/1‑100 (MQ=255)
                                                  gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAACGTCGGTTTTTTACATCCGGATTCGGACAAGGCTGAATATGACGATGCCCCGTTGaaa                                         >  2:669835/1‑99 (MQ=18)
                                                   agCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATAACCCAGTGAAAGt                                       >  1:360351/1‑100 (MQ=255)
                                                    gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                      <  2:1633276/100‑1 (MQ=255)
                                                    gCCACCAACGCAGGGCTGGAAAGTAAAAATCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaaat                                        <  2:1962587/98‑2 (MQ=25)
                                                     ccaccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtat                                     >  1:2246376/1‑100 (MQ=255)
                                                     ccaccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtat                                     >  2:1978780/1‑100 (MQ=255)
                                                      caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                    <  2:2211153/100‑1 (MQ=255)
                                                         caACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATaa                                 >  2:303452/1‑100 (MQ=255)
                                                           aCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAAt                               <  1:2097957/100‑1 (MQ=255)
                                                            cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATaaa                                >  2:2040295/1‑98 (MQ=255)
                                                            cTCAGAGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATaaa                                <  1:68668/98‑1 (MQ=255)
                                                             tCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                             <  1:583403/100‑1 (MQ=255)
                                                              cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGt                            >  1:190069/1‑100 (MQ=255)
                                                               aGGGCTGGAAAGTAAAAAACAGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAAGATGACACAGTGAAAGTATATAAATTGt                            >  2:391605/1‑99 (MQ=25)
                                                                gggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGtga                          >  1:2036737/1‑98 (MQ=255)
                                                                  gCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACTATGACCCAGAAAAAGTATAAAAATCTTCAc                         >  2:1063509/1‑99 (MQ=25)
                                                                  gCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACt                        >  2:1820261/1‑100 (MQ=255)
                                                                   cTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTg                       >  1:1648125/1‑100 (MQ=255)
                                                                    tGGAAAGTAAAAAACCGACGCAAAGTCTGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCTTGTTAAAGTATATAAATCGTCACTga                      <  2:1139333/100‑2 (MQ=37)
                                                                    tGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGc                      <  1:532920/100‑1 (MQ=255)
                                                                      gAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGa                    >  2:805883/1‑100 (MQ=255)
                                                                       aaaGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                   <  2:1614059/100‑1 (MQ=255)
                                                                       aaaGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                   <  1:1414672/100‑1 (MQ=255)
                                                                       aaaGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                   <  1:697143/100‑1 (MQ=255)
                                                                       aaaGTAAAAAACCGACGAAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                   <  1:353815/100‑1 (MQ=255)
                                                                        aaGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGata                  <  1:805883/100‑1 (MQ=255)
                                                                             aaaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatatat               >  1:2178558/1‑99 (MQ=255)
                                                                              aaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAc             >  1:90288/1‑100 (MQ=255)
                                                                               aaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc            >  1:2246/1‑100 (MQ=255)
                                                                                  aCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGaa         >  2:1680344/1‑100 (MQ=255)
                                                                                  aCCGACGCAAAGTCGGATTTTTAAATCCGGATTCGGACAAGGATTAATATGATGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGa          >  2:315409/1‑99 (MQ=255)
                                                                                   ccGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg        >  2:1692814/1‑100 (MQ=255)
                                                                                   ccGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg        <  1:2098377/100‑1 (MQ=255)
                                                                                   ccGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg        >  1:246603/1‑100 (MQ=255)
                                                                                    cGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg        <  2:1749729/99‑1 (MQ=255)
                                                                                    cGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGt       <  2:2102869/100‑1 (MQ=255)
                                                                                    cGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGt       >  2:1559840/1‑100 (MQ=255)
                                                                                     gACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGCGAAAGTATATAAATCGTCACTTCGATATATACCGaa         >  1:995148/1‑97 (MQ=255)
                                                                                      aCGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATACATACCGAAGTGc     >  2:1848695/1‑100 (MQ=255)
                                                                                       cGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGct    >  2:2081726/1‑100 (MQ=255)
                                                                                       cGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGct    >  1:2187705/1‑100 (MQ=255)
                                                                                        gCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc   <  1:629714/100‑1 (MQ=255)
                                                                                        gCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc   <  2:48920/100‑1 (MQ=255)
                                                                                        gCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc   <  1:270185/100‑1 (MQ=255)
                                                                                         cAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcc  >  2:416633/1‑100 (MQ=255)
                                                                                         cAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcc  >  1:391378/1‑100 (MQ=255)
                                                                                         cAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcc  <  1:138023/100‑1 (MQ=255)
                                                                                         cAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcc  >  1:103396/1‑100 (MQ=255)
                                                                                         cAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGct    >  2:1585561/1‑98 (MQ=255)
                                                                                         aaaaGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcc  >  2:455339/2‑100 (MQ=255)
                                                                                           aaGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc   <  1:711770/97‑1 (MQ=255)
                                                                                                  |                                                                                           
GTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCC  >  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: