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,0420TC99.1% 483.8 / ‑4.6 137intergenic (‑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/0);  new base C (70/64);  total (73/64)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

TGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG  >  NC_000913/3647947‑3648137
                                                                                               |                                                                                               
tGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATgg                                                                                             <  1:3386708/100‑1 (MQ=255)
  cgcATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTa                                                                                           >  2:3319722/1‑100 (MQ=255)
  cgcATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTa                                                                                           >  1:2691548/1‑100 (MQ=255)
    cATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTAtt                                                                                         >  2:2347946/1‑100 (MQ=255)
     aTTGCCAGAGTGAAAATGCCTCAGTGCTAACTGGATAATCATTCAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTAtt                                                                                         <  1:3299313/99‑1 (MQ=255)
        gCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCa                                                                                     >  2:1319072/1‑100 (MQ=255)
        gCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCa                                                                                     >  1:51373/1‑100 (MQ=255)
        gCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCa                                                                                     >  2:546240/1‑100 (MQ=255)
        gCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCa                                                                                     >  1:2525657/1‑100 (MQ=255)
         ccAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGATATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCAt                                                                                    >  2:1396440/1‑100 (MQ=255)
          caagggtaaaaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCa                                                                                     <  2:1748010/91‑1 (MQ=255)
          cAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCAtt                                                                                   <  1:3319972/100‑1 (MQ=255)
           agagTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCGCGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTg                                                                                  <  2:2395293/100‑1 (MQ=255)
           agagTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTg                                                                                  <  2:2032813/100‑1 (MQ=255)
           agagTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTg                                                                                  <  2:436701/100‑1 (MQ=255)
             agTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGAt                                                                                <  2:197682/100‑1 (MQ=255)
              gTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATg                                                                               >  2:1053820/1‑100 (MQ=255)
              gTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATg                                                                               <  2:12590/100‑1 (MQ=255)
                gAAAATGCCTCAGAGCTAACTGGATAATCATACCGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCg                                                                             >  1:2031630/1‑100 (MQ=255)
                gAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCg                                                                             >  1:790579/1‑100 (MQ=255)
                gAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCg                                                                             <  2:1665537/100‑1 (MQ=255)
                  aaaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGtt                                                                           <  2:1762443/100‑1 (MQ=255)
                   aaTGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTa                                                                          <  2:1107833/100‑1 (MQ=255)
                       ccTCAGCGCCAACTGGATAATCATGCAGTACATGCAGGGTATAAAACGAGCTCCTCCTTGCAATAGTTTCAGCATGATATTCACCTTGATGCGTTAgat                                                                       >  1:116263/1‑99 (MQ=255)
                       ccTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAgatg                                                                      >  1:719841/1‑100 (MQ=255)
                       ccTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAgatg                                                                      >  2:602579/1‑100 (MQ=255)
                          cAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGatgatg                                                                   <  2:532978/100‑1 (MQ=255)
                           agagCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGatgatg                                                                   <  1:2504047/99‑1 (MQ=255)
                           agagCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATgg                                                                  <  2:2713828/100‑1 (MQ=255)
                            gagCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGc                                                                 >  1:2378820/1‑100 (MQ=255)
                                 aaCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtct                                                            <  1:1442760/100‑1 (MQ=255)
                                 aaCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtct                                                            >  2:1524052/1‑100 (MQ=255)
                                 aaCTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTAtct                                                            <  2:2416660/100‑1 (MQ=255)
                                   cTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCa                                                          >  2:2118782/1‑100 (MQ=255)
                                    tGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCAc                                                         <  1:1152085/100‑1 (MQ=255)
                                    tGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCAc                                                         <  2:286128/100‑1 (MQ=255)
                                     ggATAATCATACAGTACATGTAGGTTATAAAACCAGCACGTCTTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATTATGGCTAGCTCACt                                                        >  1:2195359/1‑100 (MQ=255)
                                     ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCaat                                                        >  2:1381304/1‑98 (MQ=255)
                                     ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACt                                                        >  2:694713/1‑100 (MQ=255)
                                     ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACt                                                        >  1:3254045/1‑100 (MQ=255)
                                     ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACt                                                        >  2:2420902/1‑100 (MQ=255)
                                     ggATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACt                                                        <  1:1018422/100‑1 (MQ=255)
                                       aTAATCATACAGTACATGCAGGTTATAAAACCAACACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATAGATGCGTTAGATGATGGCTATCTCACTcc                                                      >  1:2031547/1‑100 (MQ=255)
                                        tAATCATACAGTATATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTcc                                                      <  1:265574/99‑1 (MQ=255)
                                        tAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCa                                                     >  2:3004621/1‑100 (MQ=255)
                                           tCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTc                                                  <  2:2795065/100‑1 (MQ=255)
                                           tCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTc                                                  <  2:190261/100‑1 (MQ=255)
                                           tCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTc                                                  <  1:2845127/100‑1 (MQ=255)
                                            cATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCa                                                 >  1:118023/1‑100 (MQ=255)
                                             aTACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCag                                                >  1:236938/1‑100 (MQ=255)
                                              tACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCa                                                 >  2:3059472/1‑98 (MQ=255)
                                                cAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGc                                             >  1:1027250/1‑100 (MQ=255)
                                                cAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGc                                             <  2:3069749/100‑1 (MQ=255)
                                                 aGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGcc                                            >  1:77285/1‑100 (MQ=255)
                                                   tACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccac                                          <  2:1301326/100‑1 (MQ=255)
                                                    aCATGCATGTTATAAATCCATTACGTCCTTGCAATAGTGTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccacc                                         <  2:1125078/100‑1 (MQ=255)
                                                    aCATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccacc                                         <  1:1867253/100‑1 (MQ=255)
                                                     cgtGCTGGTTTTATAACCTGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccacca                                        >  2:1742687/3‑100 (MQ=255)
                                                     cATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccacca                                        <  2:734523/100‑1 (MQ=255)
                                                     cATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGccacca                                        >  2:2710124/1‑100 (MQ=255)
                                                       tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAAc                                      >  2:2677867/1‑100 (MQ=255)
                                                       tGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCATTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAAc                                      >  2:2336997/1‑100 (MQ=255)
                                                        gCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACt                                     >  1:3053204/1‑100 (MQ=255)
                                                         cAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTc                                    >  2:2776385/1‑100 (MQ=255)
                                                          aGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCa                                   >  2:1583521/1‑100 (MQ=255)
                                                            gTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAgg                                 <  2:1188231/100‑1 (MQ=255)
                                                            gTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAgg                                 <  1:994736/100‑1 (MQ=255)
                                                              tataAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGc                               >  1:2426838/1‑100 (MQ=255)
                                                              tataAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGc                               >  2:2246481/1‑100 (MQ=255)
                                                              tataAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGc                               <  2:655030/100‑1 (MQ=255)
                                                               ataAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGTATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCt                              <  1:2168920/100‑1 (MQ=255)
                                                               ataAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCt                              >  1:763679/1‑100 (MQ=255)
                                                                taAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTg                             >  2:1879711/1‑100 (MQ=255)
                                                                taAAACCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTg                             >  1:1720790/1‑100 (MQ=255)
                                                                 aaaaCCAGCACGTCCTTGCACTAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTgg                            >  1:2078511/1‑100 (MQ=255)
                                                                    aCCAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGaaa                         <  2:1179700/100‑1 (MQ=255)
                                                                     ccAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAg                        <  2:938102/100‑1 (MQ=255)
                                                                     ccAGCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAg                        <  2:2635522/100‑1 (MQ=255)
                                                                      ctgCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAg                        <  1:3051574/97‑1 (MQ=255)
                                                                        gCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTaa                     <  1:1564991/100‑1 (MQ=255)
                                                                        gCACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTaa                     >  2:3390813/1‑100 (MQ=255)
                                                                         cACGTCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTaaa                    <  1:1663350/100‑1 (MQ=255)
                                                                             tCCTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGACACCAACTCAGGGCTGGAAAGTAAAAAAc                <  1:2875200/100‑1 (MQ=255)
                                                                              ccTTGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAAcc               >  2:926609/1‑100 (MQ=255)
                                                                                 tGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGAc            >  1:2761158/1‑100 (MQ=255)
                                                                                 tGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGAc            <  2:1334493/100‑1 (MQ=255)
                                                                                 tGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGAc            >  1:976856/1‑100 (MQ=255)
                                                                                 tGCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGAc            <  2:51373/100‑1 (MQ=255)
                                                                                  gCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACg           <  2:1303568/100‑1 (MQ=255)
                                                                                  gCAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACg           >  1:3166831/1‑100 (MQ=255)
                                                                                   cAATAGTTTCAGCCTGGTATTAGCATTGATGCGTTAGATGATGGATATCTCACTCAAGTCAGAGCCACAAACTCAGGGCTGGAAAGTAAAAAACCGACGc          >  1:1054667/1‑100 (MQ=25)
                                                                                   cAATAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGc          >  2:714669/1‑100 (MQ=255)
                                                                                    aaTAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCa         >  2:2215066/1‑100 (MQ=255)
                                                                                      tAGTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCaa        >  2:1962317/1‑99 (MQ=255)
                                                                                        gTTTCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCCATCTCGCACCAGTCAGAGCGACTAAAAAAGAGAAGGAAAGTAAAAAACCGACGCAAAGt     >  2:1467559/1‑100 (MQ=255)
                                                                                        gTTTCAGCATGCTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAACACCGACGCAAAGt     >  1:634090/1‑100 (MQ=255)
                                                                                         tttCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTc    <  2:2632153/100‑1 (MQ=255)
                                                                                           tCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCgg  <  1:2990422/100‑1 (MQ=255)
                                                                                           tCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCgg  >  2:1792080/1‑100 (MQ=255)
                                                                                           tCAGCATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCgg  >  1:1418114/1‑100 (MQ=255)
                                                                                               |                                                                                               
TGCGCATTGCCAGAGTGAAAATGCCTCAGAGCTAACTGGATAATCATACAGTACATGCAGGTTATAAAACCAGCACGTCCTTGCAATAGTTTCAGTATGGTATTAGCATTGATGCGTTAGATGATGGCTATCTCACTCCAGTCAGAGCCACCAACTCAGGGCTGGAAAGTAAAAAACCGACGCAAAGTCGG  >  NC_000913/3647947‑3648137

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