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.93.5% 300.8 / 14.7 77intergenic (‑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 (3/2);  new base . (37/35);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTC  >  NC_000913/3648040‑3648232
                                                                                                  |                                                                                              
aGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGATAAAAAAAAAAAAACCGACGCAAAGTCGGt                                                                                                <  2:38669/99‑1 (MQ=25)
  tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttt                                                                                             >  2:1001573/1‑100 (MQ=255)
  tATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACCCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttt                                                                                             <  2:254814/100‑1 (MQ=255)
   aTGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttt                                                                                            <  1:503926/100‑1 (MQ=255)
    tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                           <  2:2057236/100‑1 (MQ=255)
    tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                           <  1:11650/100‑1 (MQ=255)
    tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                           <  1:506908/100‑1 (MQ=255)
    tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                           <  2:837981/100‑1 (MQ=255)
     ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttta                                                                                          >  2:270443/1‑99 (MQ=255)
     ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttta                                                                                          <  1:358547/100‑2 (MQ=255)
     ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttta                                                                                          >  1:1642118/1‑99 (MQ=255)
     ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttta                                                                                          <  1:1510963/100‑2 (MQ=255)
     ggTATTAGCATTGATGCGTTAGATAAGAGCTATCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                           <  2:837830/99‑1 (MQ=14)
     gataatAGCATTGAGGCGTTAGATGATGGCTATCTCACTCCAGCCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttt                                                                                          <  2:1579038/95‑1 (MQ=17)
      ggagtAGCATTGATGCGTTGGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACTTTGCGTCGGttttttac                                                                                         <  2:996023/96‑3 (MQ=18)
      gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTa                                                                                         >  2:936897/1‑100 (MQ=255)
      gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACAGAAAGTACGttttttat                                                                                         >  1:1740926/1‑98 (MQ=21)
        aTTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttacat                                                                                       <  1:1972947/100‑5 (MQ=255)
         ttAGCATTGATGCGTTAGATGATGGCTATCTCACTCCGGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttt                                                                                          >  2:1395222/1‑96 (MQ=255)
         ttAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATc                                                                                     <  1:200599/100‑1 (MQ=255)
         ttAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATc                                                                                     >  1:759078/1‑100 (MQ=255)
          tAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATcc                                                                                    >  2:1349537/1‑100 (MQ=255)
           aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCg                                                                                   <  1:1058582/100‑1 (MQ=255)
            gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCTGAGTCGG‑TTTTTTACATCCg                                                                                   >  2:1948534/1‑99 (MQ=21)
             cATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGa                                                                                 >  1:1702605/1‑100 (MQ=255)
                 gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACGTc                                                                                     >  2:1923300/1‑93 (MQ=255)
                 gATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCg                                                                             >  2:2449501/1‑100 (MQ=255)
                  aTGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCgg                                                                            <  2:1990549/100‑1 (MQ=255)
                   tGCGTTAGATGATGGCTATCTCACTCCAGTCAGATCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGCATTCGGa                                                                           <  2:1370496/100‑1 (MQ=18)
                   tGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGa                                                                           >  2:2388285/1‑100 (MQ=255)
                   tGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGa                                                                           <  1:786918/100‑1 (MQ=255)
                        tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAgg                                                                      <  1:1773678/100‑1 (MQ=255)
                         aGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGc                                                                     <  2:1854735/100‑1 (MQ=255)
                         aGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGc                                                                     <  2:840506/100‑1 (MQ=255)
                             gatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTaa                                                                 <  1:432776/100‑1 (MQ=255)
                               tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAata                                                               <  1:1094931/100‑1 (MQ=255)
                               tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAata                                                               <  1:803851/100‑1 (MQ=255)
                               tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAata                                                               >  2:79220/1‑100 (MQ=255)
                               tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAata                                                               <  2:622690/100‑1 (MQ=255)
                               tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAata                                                               >  2:1510463/1‑100 (MQ=255)
                                ggCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAACGCTTAatat                                                              >  1:1787765/1‑100 (MQ=255)
                                 gCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑ATTTTTACATCCGGATTCGGACAAGGCTTAatat                                                              >  2:1572753/1‑99 (MQ=255)
                                  cTATCTCACTTCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATg                                                             <  2:1704112/100‑1 (MQ=255)
                                   tATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                           >  2:2557747/1‑100 (MQ=255)
                                   tATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                           <  2:1273268/100‑1 (MQ=255)
                                   tATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGAc                                                           >  1:501750/1‑100 (MQ=255)
                                    aTCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACg                                                          <  1:1194111/100‑1 (MQ=255)
                                    aTCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACg                                                          >  2:286325/1‑100 (MQ=255)
                                     tctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGa                                                         >  2:1879576/1‑100 (MQ=255)
                                       tcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCACAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATg                                                       <  1:484326/100‑1 (MQ=255)
                                         aCTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAGACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAc                                                     <  1:2495556/100‑1 (MQ=255)
                                          cTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAcc                                                    >  1:2334106/1‑100 (MQ=255)
                                            ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCa                                                  <  2:2495380/100‑1 (MQ=255)
                                             cAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAg                                                 <  1:1120936/100‑1 (MQ=255)
                                               gTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTg                                               <  1:1998539/100‑1 (MQ=255)
                                               gTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTg                                               >  2:445433/1‑100 (MQ=255)
                                                 cAGCGCCACCAACTCAGGGCTGGAAAGTAAAAACCCGACGCCGAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTg                                               <  2:652254/98‑1 (MQ=25)
                                                 cAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaa                                             <  1:2449514/100‑1 (MQ=255)
                                                  agagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaaa                                            >  1:958011/1‑100 (MQ=255)
                                                   gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                           <  2:1491675/100‑1 (MQ=255)
                                                   gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                           <  2:439261/100‑1 (MQ=255)
                                                   gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                           <  2:1152565/100‑1 (MQ=255)
                                                   gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATGCGGACAAGGCTTAATATGACGATGACCCAGTGaaan                                           <  2:1008609/100‑2 (MQ=255)
                                                    agCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGt                                          <  1:577003/100‑1 (MQ=255)
                                                     gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                         <  1:1879585/100‑1 (MQ=255)
                                                       caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtata                                       <  1:499178/100‑1 (MQ=255)
                                                            aCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAAt                                  >  1:971068/1‑100 (MQ=255)
                                                               cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGt                               >  1:1578532/1‑100 (MQ=255)
                                                               cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGt                               >  1:817751/1‑100 (MQ=255)
                                                                    cTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTg                          >  1:1900139/1‑100 (MQ=255)
                                                                     tGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGc                         <  1:79220/100‑1 (MQ=255)
                                                                     tGGAAAGGAAAATACCGTTGCAATTTC‑GTTTTTTTACATCCGGATTCGGATAGGGCTAAATCTGACGATGCCCTAGTGAAAGTCTATAAATCGTCACTGc                         <  2:788599/100‑1 (MQ=255)
                                                                      ggAAAGTAAAAGGCCGACTGAAAGGCGG‑GTTTTGACATCCGGATTCGGACAAGGCTGAATATGCCGATCACCCAGGGAAAGTATATAAATCGTCACTGc                         <  2:1523349/99‑1 (MQ=255)
                                                                      ggAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCg                        >  1:613931/1‑100 (MQ=255)
                                                                        naaGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                      >  2:910423/2‑100 (MQ=255)
                                                                           gTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatata                   >  1:2270779/1‑100 (MQ=255)
                                                                           gTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatata                   >  2:1812336/1‑100 (MQ=255)
                                                                           gTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatata                   >  1:2398715/1‑100 (MQ=255)
                                                                             aaaaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                 >  1:232081/1‑100 (MQ=255)
                                                                             aaaaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                 >  2:1051110/1‑100 (MQ=255)
                                                                              aaaaaCCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                 >  1:2321976/1‑100 (MQ=255)
                                                                              aaaaaCCGACGCAAAGTCGG‑TTTTTTACCTCCGGAATCGAACCCGGCTTAATATGACGATGACGCCGTGAAAATATATAAATAGTCACTGCGatatatgc                >  2:289933/1‑98 (MQ=25)
                                                                              aaaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAc                >  1:2362126/1‑100 (MQ=255)
                                                                              aaaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAc                >  1:1242208/1‑100 (MQ=255)
                                                                               aaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAcc               >  1:58894/1‑100 (MQ=255)
                                                                                 aaCCGACGCAAATTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGGCACTGCGATATATACCGa             <  2:530449/100‑1 (MQ=255)
                                                                                 aaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGa             <  2:527620/100‑1 (MQ=255)
                                                                                   ccGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg           <  1:1740118/100‑1 (MQ=255)
                                                                                   ccGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg           <  1:2013052/100‑1 (MQ=255)
                                                                                   ccGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg           >  2:2048963/1‑100 (MQ=255)
                                                                                    ngACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGt          >  1:747656/2‑100 (MQ=255)
                                                                                    cGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATCTACCGAAg           <  2:716613/99‑1 (MQ=255)
                                                                                       cGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGct       >  2:1101035/1‑100 (MQ=255)
                                                                                          aaaGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccc    >  2:677610/1‑100 (MQ=255)
                                                                                            aGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccctc  >  2:340256/1‑100 (MQ=255)
                                                                                                  |                                                                                              
AGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTC  >  NC_000913/3648040‑3648232

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