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.95.8% 290.5 / 6.2 71intergenic (‑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 (2/1);  new base . (32/36);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-01

AGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCC  >  NC_000913/3648040‑3648233
                                                                                                  |                                                                                               
aGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtt                                                                                                <  1:675351/100‑1 (MQ=255)
aGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCCCACTCCAGGCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtt                                                                                                <  1:550088/100‑1 (MQ=255)
  tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttt                                                                                              <  1:32063/100‑1 (MQ=255)
    tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            >  1:266683/1‑100 (MQ=255)
    tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            <  1:295787/100‑1 (MQ=255)
    tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            <  2:1124799/100‑1 (MQ=255)
    tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            <  2:1602205/100‑1 (MQ=255)
    tGGTATTAGCATGGATGCGTTAGATGATGGCTATCTCACTCCAGACAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            <  2:304862/100‑1 (MQ=255)
     ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTTCAGTCAGAGCCACCAACTTAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGtttttta                                                                                           <  2:1333094/100‑2 (MQ=21)
     ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGtttttta                                                                                           <  2:1312091/100‑2 (MQ=21)
     ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttta                                                                                           <  1:1528251/100‑2 (MQ=255)
     ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttta                                                                                           <  1:38762/100‑2 (MQ=255)
      gTACTAGTATTGATGCGTTAGATGATTGGTCTCTCACTCCAGTCAGAGCCACCATCTCAGGGCTGGAAAGTAAAAAACCGACGCGAAGTCGGtttttt                                                                                            <  2:1450289/98‑1 (MQ=17)
        aTTAGCATTGATGCGTTAGATGATGGCTATCTCATTCCAGTCAGCGCCACCCACCCAGAGCTGGAAATTAAAAAACCGACTTTGCGTCGGttttttactt                                                                                        <  2:2218069/100‑5 (MQ=21)
          tAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATcc                                                                                     >  1:878542/1‑100 (MQ=255)
          tAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATcc                                                                                     >  1:333921/1‑100 (MQ=255)
          tAGCATTAATGAGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttaca                                                                                         >  2:1249345/1‑94 (MQ=255)
           agtaTTGAGGCGTTAGATGATAGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTcc                                                                                     <  2:1992401/97‑1 (MQ=14)
           aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCg                                                                                    <  1:2004664/100‑1 (MQ=255)
            gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCgg                                                                                   >  1:989137/1‑100 (MQ=255)
             cATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGa                                                                                  >  1:1674206/1‑100 (MQ=255)
               ttGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGAtt                                                                                >  2:1224632/1‑100 (MQ=255)
                 gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCg                                                                              <  1:2033017/100‑1 (MQ=255)
                 gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCg                                                                              >  1:2302126/1‑100 (MQ=255)
                 gATGCGTTAGATGATGGCTACCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCg                                                                              <  1:2069787/100‑1 (MQ=255)
                  aTGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCgg                                                                             >  1:1788097/1‑100 (MQ=255)
                     cGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACa                                                                          >  1:927611/1‑100 (MQ=255)
                     cGTTAGATGATGGCTATCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTCCTGATTCgg                                                                             >  2:440237/1‑98 (MQ=255)
                          gctgatgGCCATCTCCCTCCAGTGAGAGCCGCCACCTCAGGGCGGGAGAGGCATACACCGCCGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCt                                                                     <  2:125856/98‑1 (MQ=16)
                           atgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGc                                                                      <  2:708857/98‑1 (MQ=255)
                              atgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAat                                                                 <  2:2140523/100‑1 (MQ=255)
                               tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAata                                                                >  1:2097113/1‑100 (MQ=255)
                               tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAat                                                                 <  1:1031375/99‑1 (MQ=255)
                                ggCTATCTCACTCCAGTCAGAGCCACTAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAatat                                                               <  2:17130/100‑1 (MQ=255)
                                ggCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAatat                                                               <  1:1847006/100‑1 (MQ=255)
                                ggCTACCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAatat                                                               <  1:976294/100‑1 (MQ=255)
                                      ctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGAt                                                         >  1:1792961/1‑100 (MQ=255)
                                       tcACTCCAGTCAGAGCCACCGACTCAGGGCTGGAAAGTAAACAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGCCTTAATATGACg                                                           <  2:1066173/97‑1 (MQ=18)
                                        cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGa                                                       >  1:2001667/1‑100 (MQ=255)
                                          cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                        >  2:2009010/1‑97 (MQ=255)
                                          cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAc                                                      >  1:327625/1‑99 (MQ=255)
                                           tCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAccc                                                    <  1:4213/100‑1 (MQ=255)
                                            ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCa                                                   <  2:122330/100‑1 (MQ=255)
                                            ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCa                                                   <  1:1220107/100‑1 (MQ=255)
                                            ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCa                                                   >  1:1841826/1‑100 (MQ=255)
                                               gTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTg                                                >  2:1753328/1‑100 (MQ=255)
                                                 cAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaa                                              >  1:1878508/1‑100 (MQ=255)
                                                   gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                            <  1:1627682/100‑1 (MQ=255)
                                                     gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                          >  1:950591/1‑100 (MQ=255)
                                                     gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                          >  1:122437/1‑100 (MQ=255)
                                                     gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                          >  1:979900/1‑100 (MQ=255)
                                                     gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                          <  1:636965/100‑1 (MQ=255)
                                                     gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                          <  2:285057/100‑1 (MQ=255)
                                                     gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                          <  1:611653/100‑1 (MQ=255)
                                                       caccaACTCAGGGCTGGAAAGTAAAAACCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                        <  1:1224639/100‑1 (MQ=255)
                                                       caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                        <  1:1753342/100‑1 (MQ=255)
                                                       caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCAGTTTTTTACATCCGGATTCGTACACGGCTTAATATGACGATGACCCAGTGAAAg                                            <  2:1970478/96‑1 (MQ=21)
                                                         ccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatata                                      <  1:602835/100‑1 (MQ=255)
                                                          caACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATaa                                     >  1:137340/1‑100 (MQ=255)
                                                             cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                                  >  1:337102/1‑100 (MQ=255)
                                                              tCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                                 <  2:245580/100‑1 (MQ=255)
                                                                 gggCTGGAAAGTAAAAAAACGACGCAAAGTCGGTTTTTTAAATCCGGATTAGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGt                                >  2:2020379/1‑98 (MQ=37)
                                                                  ggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCAc                             <  1:1825039/100‑1 (MQ=255)
                                                                    cTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACt                            >  1:777684/1‑100 (MQ=255)
                                                                    cTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACt                            >  2:1720719/1‑99 (MQ=255)
                                                                        aaaGTAAAAAACCAACGCAAAGTCGGTTTTTTACATCCCGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                       >  2:1468050/1‑100 (MQ=39)
                                                                         aaGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGata                      <  2:1947113/100‑1 (MQ=255)
                                                                          aGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatat                     <  2:2314232/100‑1 (MQ=255)
                                                                           gTAAAAAACCGACTTTGCGTCGGTTTTTTACATCCGGATTCGGATAAGGCTTAATATGACGATGACCCAGTGACAGTATATAAATCGTCACTGCGctata                    >  2:1134195/1‑100 (MQ=25)
                                                                           gTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatata                    <  1:1220865/100‑1 (MQ=255)
                                                                            tAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatatat                   >  1:769123/1‑100 (MQ=255)
                                                                             naaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                  >  1:2345766/2‑100 (MQ=255)
                                                                             aaaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCCTAATATGCCGATGACCCACTGAAAGTATATAAATCGCCACTGAGAtatata                  >  2:225672/1‑100 (MQ=37)
                                                                               aaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc                <  2:211464/100‑1 (MQ=255)
                                                                                     gACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGt           <  1:1582851/99‑1 (MQ=255)
                                                                                     gACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTg          >  1:1282628/1‑100 (MQ=255)
                                                                                     gACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTg          <  1:1349163/100‑1 (MQ=255)
                                                                                     gACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTg          >  2:399350/1‑100 (MQ=255)
                                                                                        gCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc       <  2:1123915/100‑1 (MQ=255)
                                                                                        gCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc       <  2:666104/100‑1 (MQ=255)
                                                                                          tgcgTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc     <  2:753178/97‑1 (MQ=255)
                                                                                          aaaGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc     <  1:646306/100‑1 (MQ=255)
                                                                                          aaaGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGccccc     >  2:2292938/1‑100 (MQ=255)
                                                                                           aaGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccct    <  2:927609/100‑1 (MQ=255)
                                                                                           aaGTCGGTGTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc       <  2:576688/97‑1 (MQ=255)
                                                                                             gTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccctcc  >  1:413544/1‑100 (MQ=255)
                                                                                                  |                                                                                               
AGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCC  >  NC_000913/3648040‑3648233

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: