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.97.7% 353.9 / 1.5 86intergenic (‑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 (1/1);  new base . (41/43);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

CAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCC  >  NC_000913/3648039‑3648233
                                                                                                   |                                                                                               
cAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGt                                                                                                 >  1:2804999/1‑100 (MQ=255)
  gTATGGTATTAGCATTGATGCGTTAGATGATGGTTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGCCGGtt                                                                                                >  1:1925317/1‑99 (MQ=255)
   tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCAATCCAGTCAGAGCCACCAACTCATGGATGAAAAAAAAAAAACCGACGCAAAGTCGGtttt                                                                                              <  1:316683/100‑1 (MQ=25)
    aTGGCATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttt                                                                                             <  2:2484520/100‑1 (MQ=255)
     tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTGACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtatttt                                                                                            >  2:692759/1‑100 (MQ=255)
     tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            <  2:1832536/100‑1 (MQ=255)
     tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            <  2:1194606/100‑1 (MQ=255)
     tGGTAGTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAACACCGACGCAAAGTCGGtttttt                                                                                            <  1:1100806/100‑1 (MQ=255)
      ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttta                                                                                           <  1:2834614/100‑2 (MQ=255)
      ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttta                                                                                           <  1:1166518/100‑2 (MQ=255)
       gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttac                                                                                          <  2:1285757/100‑3 (MQ=255)
        tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttaca                                                                                         >  2:1249719/1‑97 (MQ=255)
        tATTAGCATTGATGCGTTAGATGATCGCTATCTCACTCCAGTCAGCGCCACCACCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttact                                                                                         <  1:2322466/100‑4 (MQ=21)
            aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCg                                                                                    <  1:2559994/100‑1 (MQ=255)
            aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTAGAAAATAAAAAACCGACGCAAAGTCGGTTTTTTACATCCg                                                                                    >  1:1856269/1‑100 (MQ=21)
             gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCgg                                                                                   >  2:1334608/1‑100 (MQ=255)
             gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCgg                                                                                   >  2:23024/1‑100 (MQ=255)
              cATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGa                                                                                  >  1:487862/1‑100 (MQ=255)
               aTTGATTCGTTAGATGATGGCTATCTCACTCCAGTCATAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGAt                                                                                 <  1:1426527/100‑1 (MQ=18)
                ttGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGAtt                                                                                <  2:1327881/100‑1 (MQ=255)
                 tGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTc                                                                               >  2:2898310/1‑100 (MQ=255)
                   aTGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCgg                                                                             <  1:1537568/100‑1 (MQ=255)
                   aTGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCgg                                                                             <  2:1566398/100‑1 (MQ=255)
                     gCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTAACTCCGGGTTCg                                                                              >  2:1000744/1‑97 (MQ=18)
                     gCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGAc                                                                           <  1:1360957/100‑1 (MQ=255)
                        ttAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCAGAtctgaagagc                                                                        >  2:2914623/1‑91 (MQ=17)
                         tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAgg                                                                       >  2:2764390/1‑100 (MQ=255)
                          aGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGc                                                                      <  2:1317454/100‑1 (MQ=255)
                           gatgatGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCt                                                                     <  1:535407/100‑1 (MQ=255)
                            atgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCtt                                                                    >  2:1396746/1‑100 (MQ=255)
                             tgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTa                                                                   >  1:2553383/1‑100 (MQ=255)
                             tgatgGCTACCTAACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTa                                                                   <  1:545047/100‑1 (MQ=18)
                             ctcttccgaTCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTa                                                                   <  2:1996141/92‑1 (MQ=18)
                              gatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTaa                                                                  <  2:2387124/100‑1 (MQ=255)
                               atgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAat                                                                 <  2:2466798/100‑1 (MQ=255)
                               atgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAat                                                                 <  2:1354899/100‑1 (MQ=255)
                                  gCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAATGCTTAata                                                                >  2:1890065/1‑98 (MQ=255)
                                  gCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATg                                                              <  2:2005907/100‑1 (MQ=255)
                                   cTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAACTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGa                                                             >  1:571630/1‑100 (MQ=255)
                                    tATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                            <  2:625816/100‑1 (MQ=255)
                                     aTCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACg                                                           >  2:1612639/1‑100 (MQ=255)
                                       ctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGAt                                                         <  1:2143970/100‑1 (MQ=255)
                                       ctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGAt                                                         <  1:202257/100‑1 (MQ=255)
                                       ctcACTCCAGTCACAGCCACCAACTCAGGGCGGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGAt                                                         >  2:1960382/1‑100 (MQ=21)
                                        tcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                        >  1:2839283/1‑100 (MQ=255)
                                        tcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                        >  2:2557098/1‑100 (MQ=255)
                                          aCTCCAGACAGAGCCACCAACTCAGGGCTGGAAAGTCAACAACGGAAGCAAAGACGGTTTTATACATGCGGACTCGGAGAAGACCTAATATgcccaggac                                                      >  1:2375102/1‑92 (MQ=9)
                                           cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                     >  1:1162525/1‑100 (MQ=255)
                                             ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAccc                                                    <  1:1005014/100‑1 (MQ=255)
                                              cAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGTCATGGCTTACTATGACGATgc                                                       >  2:1261499/1‑94 (MQ=18)
                                                gTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTg                                                <  2:1388887/100‑1 (MQ=255)
                                                gTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTg                                                <  2:356409/100‑1 (MQ=255)
                                                gTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTg                                                <  2:253551/100‑1 (MQ=255)
                                                     agCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGt                                           <  2:261214/100‑1 (MQ=255)
                                                      gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCCGTAAAAGta                                          >  2:332011/1‑100 (MQ=21)
                                                      gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                          >  2:344426/1‑100 (MQ=255)
                                                      gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                          <  1:1231314/100‑1 (MQ=255)
                                                        caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                        <  2:892899/100‑1 (MQ=255)
                                                        caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                        >  2:1186722/1‑100 (MQ=255)
                                                           caACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATACGACGATGACCCAGTGAAAGTATATaa                                     <  2:1772004/100‑1 (MQ=255)
                                                             aCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAAt                                   >  2:1139969/1‑100 (MQ=255)
                                                              cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGAATCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAAt                                   >  1:1088791/1‑99 (MQ=255)
                                                               tCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                                 <  2:1644738/100‑1 (MQ=255)
                                                               tCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                                 <  2:344272/100‑1 (MQ=255)
                                                                cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGt                                <  2:2490029/100‑1 (MQ=255)
                                                                cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACACCCGGATTCGGCCACCAATTAATATAACGATGACCCCGTGAAAGTCCATAAAta                                  >  2:2433246/1‑97 (MQ=17)
                                                                cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTAGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGAGGACCCAGTGAAAGTATATAAATCGt                                >  2:2771099/1‑100 (MQ=37)
                                                                  gggCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCa                              <  2:287983/100‑1 (MQ=255)
                                                                     cTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTg                           >  1:105826/1‑100 (MQ=255)
                                                                     cTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTg                           <  2:1269613/100‑1 (MQ=255)
                                                                      tGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGc                          <  2:2830248/100‑1 (MQ=255)
                                                                          aaGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGata                      <  2:463662/100‑1 (MQ=255)
                                                                          aaGTAAAAAACCGACGCAAAGTCGGTTCTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGata                      >  2:2689103/1‑100 (MQ=255)
                                                                           tgTAAAAAACCGACTTTGCGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAGAGTATATAAATCGTCACTGCGata                      >  1:1277763/2‑99 (MQ=25)
                                                                           aGTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatat                     <  1:2121647/100‑1 (MQ=255)
                                                                            gTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatata                    >  2:694748/1‑100 (MQ=255)
                                                                            gTAAAAAACCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatata                    <  2:2865032/100‑1 (MQ=255)
                                                                             tAAAAAACCGACTTTGCGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACGCCGTGAAAGGATATAAATCGTCACTGCGatata                    >  1:841855/1‑99 (MQ=25)
                                                                              aaaaaaCCGACGCAAAGTCGGTTTTTTACATTCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                  <  1:1786635/100‑1 (MQ=255)
                                                                              aaaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                  >  2:1442371/1‑100 (MQ=255)
                                                                              aaaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                  <  1:1249714/100‑1 (MQ=255)
                                                                              aaaaaaCCGACGCAAAGTCGGTTTTTTACAGCCAGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACAGCGAtatata                  >  1:1917926/1‑100 (MQ=38)
                                                                               aaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCATTGCGAtatata                  >  1:1932280/1‑99 (MQ=255)
                                                                               aaaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAc                 >  2:833995/1‑100 (MQ=255)
                                                                                aaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc                <  2:2426988/100‑1 (MQ=255)
                                                                                aaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc                <  2:367134/100‑1 (MQ=255)
                                                                                aaaaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAc                 >  2:1770520/1‑99 (MQ=255)
                                                                                  aaCCGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGa              <  2:2655157/100‑1 (MQ=255)
                                                                                     cGACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCTAGTGAAAGTATATAAATCGTCACTGCGGTATATACCGAAGt           >  2:260556/1‑100 (MQ=255)
                                                                                      gACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGt           >  1:13457/1‑99 (MQ=255)
                                                                                      gACGCAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTg          >  1:1042447/1‑100 (MQ=255)
                                                                                         ncAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc       >  1:2785388/2‑100 (MQ=255)
                                                                                          cAAAGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg            <  1:635160/94‑1 (MQ=255)
                                                                                           aaaGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc     <  2:2839302/100‑1 (MQ=255)
                                                                                           aaaGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc     <  2:27881/100‑1 (MQ=255)
                                                                                            aaGTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccct    <  1:2565625/100‑1 (MQ=255)
                                                                                              gTCGGTTTTTTACATCCGGATTCTGACAAGTCTTAATATGATGATGACCCAGTGAAAGTATATAAAGCGTCACTGCGATATATACCGAAGTGct        <  1:1865440/94‑1 (MQ=255)
                                                                                              gTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccctcc  >  2:2515938/1‑100 (MQ=255)
                                                                                              gTCGGTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATAGTAACTGCGATATATACCGAAGTGctccc     >  1:1028047/1‑97 (MQ=255)
                                                                                                   |                                                                                               
CAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCC  >  NC_000913/3648039‑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: