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% 348.7 / 1.3 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 (2/0);  new base . (48/36);  total (50/36)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

AGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGCTCCCTCC  >  NC_000913/3648040‑3648233
                                                                                                  |                                                                                               
gatctGGTATTAGCATTGATGCGTTAGATGATGGCGATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtt                                                                                                <  1:1029220/96‑1 (MQ=255)
aGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtt                                                                                                >  1:1606032/1‑100 (MQ=255)
aGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtt                                                                                                <  2:1200797/100‑1 (MQ=255)
    tGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            <  2:547857/100‑1 (MQ=255)
     ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGtttttt                                                                                            <  1:1519122/99‑1 (MQ=255)
     ggTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttagatc                                                                                           >  2:1029223/1‑95 (MQ=255)
      gTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttac                                                                                          >  1:606371/1‑98 (MQ=255)
       tATTAGTATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTTGAAAGTAAAAAACCGACTTTGCGTCGGttttttac                                                                                          <  2:214377/99‑3 (MQ=18)
       tATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttaca                                                                                         <  2:1303969/100‑4 (MQ=255)
        aTTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGttttttacat                                                                                        >  1:921976/1‑96 (MQ=255)
           aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCg                                                                                    <  2:2567144/100‑1 (MQ=255)
           aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCg                                                                                    >  2:1124880/1‑100 (MQ=255)
           aGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCg                                                                                    >  1:214761/1‑100 (MQ=255)
            gCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCgg                                                                                   >  2:279893/1‑100 (MQ=255)
                tGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTAAGTCCTGAt                                                                                 >  1:846887/1‑99 (MQ=255)
                  aTGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCgg                                                                             <  2:1715747/100‑1 (MQ=255)
                   tGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGa                                                                            <  2:214763/100‑1 (MQ=255)
                    gCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGAc                                                                           <  2:1925154/100‑1 (MQ=255)
                    gCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGAc                                                                           >  2:1819599/1‑100 (MQ=255)
                     cGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACa                                                                          <  2:416057/100‑1 (MQ=255)
                     cGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACa                                                                          <  1:2262349/100‑1 (MQ=255)
                        tAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAgg                                                                       >  1:1131204/1‑100 (MQ=255)
                           atgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCtt                                                                    <  1:2463667/100‑1 (MQ=255)
                           atgatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑ATTTTTACATCCGGATTTGGACAAGGCtt                                                                    >  2:1147066/1‑100 (MQ=18)
                             gatgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTaa                                                                  >  1:1949373/1‑100 (MQ=255)
                               ttccgaTCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAata                                                                <  1:926755/95‑1 (MQ=255)
                               tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAata                                                                >  2:1852018/1‑100 (MQ=255)
                               tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAata                                                                <  1:772306/100‑1 (MQ=255)
                               tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAata                                                                <  1:772180/100‑1 (MQ=255)
                               tgGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAata                                                                >  2:1501737/1‑100 (MQ=255)
                                    aTCTCACTCCAGTCAGAGCCACGAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACg                                                           <  1:651544/100‑1 (MQ=255)
                                     tctcACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCCAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGa                                                          >  2:1029964/1‑100 (MQ=255)
                                        cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAataagatcggaa                                                       >  2:926758/1‑91 (MQ=18)
                                        cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGa                                                       <  2:83014/100‑1 (MQ=255)
                                        cACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGa                                                       >  2:2219966/1‑100 (MQ=255)
                                         aCTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAc                                                      <  2:1100585/100‑1 (MQ=255)
                                           tCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGAccc                                                    <  1:30343/100‑1 (MQ=255)
                                            ccAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCa                                                   >  1:1310592/1‑100 (MQ=255)
                                             cAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAg                                                  <  1:898270/100‑1 (MQ=255)
                                              aGTCAGAGCCACCAACTCAGGGCTGTAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAg                                                  <  2:508353/99‑1 (MQ=255)
                                              aGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACGTACAGCATCGGACAAATCTTAATATAACGATGACACAGt                                                 >  2:2030529/1‑100 (MQ=14)
                                               gTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTg                                                >  1:300339/1‑100 (MQ=255)
                                               gTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTg                                                >  2:1420797/1‑100 (MQ=255)
                                                 cAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaa                                              >  2:2215252/1‑100 (MQ=255)
                                                 cAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGaa                                              >  2:2273709/1‑100 (MQ=255)
                                                   gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                            <  1:867334/100‑1 (MQ=255)
                                                   gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                            <  1:341751/100‑1 (MQ=255)
                                                   gagCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTGTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAg                                            <  1:1675741/100‑1 (MQ=255)
                                                     gCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGACTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGta                                          >  2:527699/1‑100 (MQ=255)
                                                      ccaccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtat                                         <  2:921979/100‑1 (MQ=255)
                                                       caccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGTTTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtat                                         >  1:2550773/1‑100 (MQ=255)
                                                         ccaACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGtatata                                      >  2:1805144/1‑100 (MQ=255)
                                                           aaCTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATaaa                                    <  1:2365229/100‑1 (MQ=255)
                                                             cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                                  >  2:383735/1‑100 (MQ=255)
                                                             cTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                                  >  1:2093787/1‑100 (MQ=255)
                                                              tCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATc                                  >  2:803319/1‑99 (MQ=255)
                                                              tCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                                 <  2:2213235/100‑1 (MQ=255)
                                                              tCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCg                                 <  2:287658/100‑1 (MQ=255)
                                                               cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGGTTAATATGACGATGACCCAGTGAAAGTATATAAATCGt                                >  2:835544/1‑100 (MQ=255)
                                                               cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGt                                >  2:2469316/1‑100 (MQ=255)
                                                               cAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTAAATTCGGAATCCGACACGCCTTAAAATTATTATGACACTTTAAAAGtataacattc                                  >  1:546050/1‑92 (MQ=2)
                                                                 gggCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCa                              >  2:1535851/1‑100 (MQ=255)
                                                                 gggCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑ATGTTTACAGCGGAATGCAGACAAGGCGTAATATGTCGATGACCCAGGGAAAGTATATAAATAGTCa                              >  1:1741578/1‑100 (MQ=16)
                                                                  ggCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCAc                             >  2:1039263/1‑100 (MQ=255)
                                                                  ggCTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCAc                             <  1:2273819/100‑1 (MQ=255)
                                                                    cTGGAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGCCACTg                           >  2:1836461/1‑100 (MQ=255)
                                                                     ccggatgTAAAAAACCGACTTTGCGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTACTATGACGATGACCCAGTAAAAGTATATAACTCGTCACTGc                          >  1:2385556/7‑100 (MQ=21)
                                                                      ggAAAGTAAAAAACCGACGCGAAGTCGG‑TTTTTTACATCCAGATTCGGACAAAGCTTAATATTATGATGATCCAGAGAAAGTATATAAATCGTCACTgtg                         >  1:1582379/1‑98 (MQ=25)
                                                                      ggAAAGTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCg                         >  2:99178/1‑100 (MQ=255)
                                                                        atcttAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGat                       <  1:2326562/96‑1 (MQ=255)
                                                                          tgTAAAAAACCGACTTTGCGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatat                     >  1:1625264/2‑100 (MQ=37)
                                                                           gTAAAAAACCGACTTTGCGTCGG‑TTTTTGACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatata                    >  1:1894833/1‑100 (MQ=37)
                                                                           gTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatata                    >  2:1681023/1‑100 (MQ=255)
                                                                           gTAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatata                    >  2:2541925/1‑100 (MQ=255)
                                                                            tAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatatat                   >  1:148596/1‑100 (MQ=255)
                                                                            tAAAAAACCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGatagat                   >  2:2326576/1‑97 (MQ=255)
                                                                             aaaaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGAtatata                  >  1:672193/1‑100 (MQ=255)
                                                                              aaaaaCCGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATAc                 <  1:2478206/100‑1 (MQ=255)
                                                                              aaaaaCCGACGCAAAGTCGG‑TTTTTTAAATCCGGAATCGGACAAGGCTTAAAATGAATTTGCCACATTGAAAGTATATCAATCTTCCCTGCGctatat                   >  1:599932/1‑98 (MQ=16)
                                                                                   ccGACGCAAAGTCGG‑TTTTTTACATCCGGATTTTGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg            <  1:1147062/100‑1 (MQ=255)
                                                                                   ccGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg            <  2:275697/100‑1 (MQ=255)
                                                                                   ccGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg            >  1:1634590/1‑100 (MQ=255)
                                                                                   ccGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAg            >  2:263029/1‑100 (MQ=255)
                                                                                    cGACGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGt           <  1:1139688/100‑1 (MQ=255)
                                                                                      aCGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGc         <  1:964073/100‑1 (MQ=255)
                                                                                       cGCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGct        <  1:1781715/100‑1 (MQ=255)
                                                                                        gCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc       <  2:1949383/100‑1 (MQ=255)
                                                                                        gCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc       <  1:2525211/100‑1 (MQ=255)
                                                                                        gCAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc       <  1:99177/100‑1 (MQ=255)
                                                                                         cAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCTAGGGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcc      >  1:1654799/1‑100 (MQ=255)
                                                                                         cAAAGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGGGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctcc      >  1:1656025/1‑100 (MQ=255)
                                                                                         acacgTCGG‑TTTTTTACATCAGGATTCGGACTAGGCTTAATATGACTATGATCCAGTGAAAGTATATAATTTTTCACTGCGATATATACCAAATTGc         >  1:1772176/5‑97 (MQ=255)
                                                                                          aaaGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctc       <  1:412420/98‑1 (MQ=255)
                                                                                           aaGTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccct    <  1:606637/100‑1 (MQ=255)
                                                                                             gTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccctcc  >  1:822679/1‑100 (MQ=255)
                                                                                             gTCGG‑TTTTTTACATCCGGATTCGGACAAGGCTTAATATGACGATGACCCAGTGAAAGTATATAAATCGTCACTGCGATATATACCGAAGTGctccctcc  >  2:1325573/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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑