Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,648,042 T→C 100% intergenic (‑254/‑486) 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,0420TC96.4% 443.9 / 2.3 135intergenic (‑254/‑486)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/3);  new base C (76/54);  total (78/57)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 135 total aligned reads displayed.

TATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGT  >  NC_000913/3647945‑3648138
                                                                                                 |                                                                                                
tATGCTCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAAGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCAt                                                                                                >  2:1492095/1‑100 (MQ=255)
 aTGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATg                                                                                               >  2:1321634/1‑100 (MQ=255)
 aTGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATg                                                                                               >  1:1141491/1‑100 (MQ=255)
  tGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATgg                                                                                              <  1:743613/100‑1 (MQ=255)
    cgcATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTa                                                                                            >  1:1652960/1‑100 (MQ=255)
    cgcATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTa                                                                                            >  1:576442/1‑100 (MQ=255)
         tGCTAGAGTGAAAATGCCTCAGATCTAACTTGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTTCAATAGTTTCAGCATGGTATTAg                                                                                        <  2:2054136/99‑1 (MQ=255)
          gCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCa                                                                                      >  2:1036563/1‑100 (MQ=255)
          gCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCa                                                                                      >  2:2179297/1‑100 (MQ=255)
           ccAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCAt                                                                                     >  1:892308/1‑100 (MQ=255)
           ccAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCAt                                                                                     >  1:229892/1‑100 (MQ=255)
               agTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGAt                                                                                 <  1:414509/100‑1 (MQ=255)
                ggtaaaaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATg                                                                                >  1:1501464/4‑100 (MQ=255)
                gTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACAAGCACGTCCTTAAAATAGTTTCAGCATGGTATTAGCATTGAt                                                                                 >  1:807654/1‑99 (MQ=255)
                    acaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGtt                                                                            <  2:1455819/98‑1 (MQ=255)
                    aaaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTGTCAGCATGGTATTAGCATTGATGCGt                                                                             <  2:2169272/99‑1 (MQ=255)
                     aaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTa                                                                           <  1:1074952/100‑1 (MQ=255)
                      aTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAg                                                                          <  1:1647986/100‑1 (MQ=255)
                       tGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCATCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAg                                                                          <  2:2117122/99‑1 (MQ=255)
                        gCCTCAGATCTAACTGGATAATCATACCGTACATGCAGTTTATAAAACCAGCACGTCCTGTCGATAGTTTCAGCATGGTATTAGCATTGATACGTTAga                                                                         <  2:1577994/99‑1 (MQ=255)
                          cTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAgatga                                                                      >  1:1987588/1‑100 (MQ=255)
                          cTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAgatga                                                                      <  1:2354418/100‑1 (MQ=255)
                             agagCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATgg                                                                   <  1:394473/100‑1 (MQ=255)
                             agagCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATgg                                                                   <  2:1628665/100‑1 (MQ=255)
                             agagCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATgg                                                                   <  2:276718/100‑1 (MQ=255)
                                   aaCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtct                                                             >  2:1181845/1‑100 (MQ=255)
                                   aaCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtct                                                             <  1:1437554/100‑1 (MQ=255)
                                    aCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtctc                                                            <  1:111375/100‑1 (MQ=255)
                                    aCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtctc                                                            >  1:104338/1‑100 (MQ=255)
                                      tGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCAc                                                          <  2:927518/100‑1 (MQ=255)
                                       ggATAATCATACAGTACATGCAGGTTATAAAACCAGTGCGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACt                                                         >  2:1276068/1‑100 (MQ=255)
                                       ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACt                                                         >  2:894545/1‑100 (MQ=255)
                                         aTAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCCTGGTATTAGCCTTGATGCGTTAGATGATTGCTATCTCACTcc                                                       <  1:1616292/100‑1 (MQ=255)
                                          tAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTca                                                       >  2:1713301/1‑98 (MQ=255)
                                            aTCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAg                                                     >  1:919599/1‑99 (MQ=255)
                                              cATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCa                                                  <  1:1701682/100‑1 (MQ=255)
                                              cATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCa                                                  >  1:625578/1‑100 (MQ=255)
                                               aTACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCag                                                 >  2:195003/1‑100 (MQ=255)
                                                  cAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGc                                              <  2:1263686/100‑1 (MQ=255)
                                                   aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGcc                                             <  1:595411/100‑1 (MQ=255)
                                                   aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGc                                              >  1:972796/1‑99 (MQ=255)
                                                      aCATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccacc                                          <  1:316470/100‑1 (MQ=255)
                                                      aCATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccacc                                          <  1:2066849/100‑1 (MQ=255)
                                                      aCATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTACAGTCAGAGccacc                                          >  2:2094186/1‑100 (MQ=255)
                                                       cgtGCTGGTTTTATAACCTGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccacca                                         >  1:132311/3‑100 (MQ=255)
                                                       cgtGCTGGTTTTATAACCTGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGTCTATCTCACTCCAGTCAGAGccacca                                         >  2:1068833/3‑100 (MQ=255)
                                                       cATTCAGGTTATAAAACCAGCACGTCCTGGCAAGATTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccacca                                         <  2:1843215/100‑1 (MQ=255)
                                                       cATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccacca                                         >  2:19897/1‑100 (MQ=255)
                                                       cATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCAGGGAATTAGCATTGATGCGTGAGATGATGGCTATCTCACTCCAGTCAGAGccacc                                          >  2:1639180/1‑99 (MQ=255)
                                                         tGCTGGTTTTATAACCTGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAAc                                       >  2:1355210/1‑100 (MQ=255)
                                                         tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAAc                                       >  2:1046980/1‑100 (MQ=255)
                                                         tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAAc                                       >  2:1668080/1‑100 (MQ=255)
                                                         tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAAc                                       >  1:468633/1‑100 (MQ=255)
                                                          gCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACt                                      <  1:1683267/100‑1 (MQ=255)
                                                            aGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCa                                    >  1:1835862/1‑100 (MQ=255)
                                                             ggTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAg                                   <  1:655523/100‑1 (MQ=255)
                                                             ggTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAg                                   >  2:889255/1‑100 (MQ=255)
                                                             ggTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAg                                   >  1:2377550/1‑100 (MQ=255)
                                                             ggTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAg                                   >  1:2172975/1‑100 (MQ=255)
                                                              gTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAgg                                  <  2:1237479/100‑1 (MQ=255)
                                                               ttATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAggg                                 <  1:1004478/100‑1 (MQ=255)
                                                                tataAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGc                                >  2:1364580/1‑100 (MQ=255)
                                                                tataAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGc                                <  1:1305973/100‑1 (MQ=255)
                                                                  taAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTg                              >  2:885214/1‑100 (MQ=255)
                                                                  taAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGACt                               >  1:106259/1‑99 (MQ=255)
                                                                  taAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGAAGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTg                              >  2:1688193/1‑100 (MQ=255)
                                                                    aaaCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGa                            <  2:1342697/100‑1 (MQ=255)
                                                                     aaCTAGCACGTTCTTGCAATAGTTTCAGCATGGTATTAGCATTGAGGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTgg                             <  2:576439/98‑1 (MQ=38)
                                                                     aaCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGaa                           >  2:914923/1‑100 (MQ=255)
                                                                     aaCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGaa                           >  2:1003379/1‑100 (MQ=255)
                                                                      aCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGaaa                          <  2:1321736/100‑1 (MQ=255)
                                                                       ccAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAg                         <  2:62492/100‑1 (MQ=255)
                                                                       ccAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAg                         <  2:1941376/100‑1 (MQ=255)
                                                                        cAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGt                        <  2:748408/100‑1 (MQ=255)
                                                                        cAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGt                        >  2:2112466/1‑100 (MQ=255)
                                                                         aGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTa                       >  2:6275/1‑100 (MQ=255)
                                                                         aGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTa                       <  2:2063428/100‑1 (MQ=255)
                                                                         aGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGAATTGATGCGTTAGATGAAGGCTATCTCACTCCAGTCAGAGCCACCAAATCAGGGCTGGaaa                          >  2:1653761/1‑97 (MQ=255)
                                                                           cACGTCCTTGCAATAGTTTCAGCATTGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAg                         <  2:445273/96‑1 (MQ=255)
                                                                             cGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTaaaaa                   >  1:955061/1‑100 (MQ=255)
                                                                             cGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTaaaaa                   >  1:890364/1‑100 (MQ=255)
                                                                             cGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTaaaaa                   >  2:1155908/1‑100 (MQ=255)
                                                                             cGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTaaaaa                   >  1:2020203/1‑100 (MQ=255)
                                                                                ccTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAcc                >  1:729400/1‑100 (MQ=255)
                                                                                 cTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCg               >  1:1313943/1‑100 (MQ=255)
                                                                                  ttGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGACTGGAAAGTAAAAAACCGa              >  1:1857998/1‑100 (MQ=255)
                                                                                   tGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGAc             >  2:1769655/1‑100 (MQ=255)
                                                                                   tGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAGCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAATAAACCGAc             >  1:350938/1‑100 (MQ=255)
                                                                                    gCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACg            <  2:668774/100‑1 (MQ=255)
                                                                                    gCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACg            <  2:10830/100‑1 (MQ=255)
                                                                                     cAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGAc             >  2:708431/1‑98 (MQ=255)
                                                                                      aaTAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGcg          <  1:307420/100‑2 (MQ=255)
                                                                                       aTAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCaa         <  1:2022407/100‑1 (MQ=255)
                                                                                          gTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGt      >  1:2268055/1‑100 (MQ=255)
                                                                                          gTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTAACTCCAGTCAGAGCCACCAACTCAGGACTGGAAAGTAAAAAACCAACgaacag       >  2:1463289/1‑94 (MQ=25)
                                                                                           tttCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTc     <  1:1401591/100‑1 (MQ=255)
                                                                                           tttCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTc     >  1:435210/1‑100 (MQ=255)
                                                                                            ttCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAACGACGCAAAGt      >  1:2160319/1‑98 (MQ=255)
                                                                                              cAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGt  >  2:2190552/1‑100 (MQ=255)
                                                                                              cAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGt  >  2:336254/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TATGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGGT  >  NC_000913/3647945‑3648138

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: