Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,260,901 C→T 15.8% intergenic (‑98/+27) dauA ← / ← prs C4‑dicarboxylic acid transporter/phosphoribosylpyrophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,260,9010CT15.8% 331.7 / 45.6 135intergenic (‑98/+27)dauA/prsC4‑dicarboxylic acid transporter/phosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (54/59);  new base T (15/6);  total (69/66)
Fisher's exact test for biased strand distribution p-value = 5.81e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

Only 100 of 144 total aligned reads displayed.

ACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTG  >  NC_000913/1260802‑1260998
                                                                                                   |                                                                                                 
acaATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACccgcc                                                                                                   <  2:1027455/100‑1 (MQ=255)
 caATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACccgccg                                                                                                  >  1:862220/1‑100 (MQ=255)
      atGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCg                                                                                             >  2:1136157/1‑100 (MQ=255)
        ggttCATCCGTGAGCGTATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGAAAAAAAAAAAAAAAAAAACCCGCCGCAGCggg                                                                                           <  1:2296391/98‑1 (MQ=255)
        gTTTCATCCGTGAGCGCATCATCTGTCAACTACATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCggg                                                                                           >  2:533867/1‑100 (MQ=255)
        gTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCggg                                                                                           >  1:243091/1‑100 (MQ=255)
        gTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAAAGGATATTACAGACAAAAAAACCCGCCGCAGCggg                                                                                           >  1:2234694/1‑100 (MQ=255)
        gTTTCAGCCGTGAACGCATCATCTGTCAACTAAATGGCAGGTCAAGGAGGCATAGGTCATACAAATGTAAAAAAAAAAAAAAAAAACCCGCCGCAGCggg                                                                                           <  1:3409426/100‑1 (MQ=255)
          ttCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTc                                                                                         <  1:259248/100‑1 (MQ=255)
          ttCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGAAATTACAAAAAAAAAAACCCGCCGCAGCGGGTc                                                                                         <  2:2514174/100‑1 (MQ=255)
           tCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCt                                                                                        >  1:3331152/1‑100 (MQ=255)
             aTCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCttt                                                                                      <  2:3151536/100‑1 (MQ=255)
              tCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTg                                                                                     <  2:801860/100‑1 (MQ=255)
              tCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTg                                                                                     <  2:203012/100‑1 (MQ=255)
              tCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTg                                                                                     <  1:2576468/100‑1 (MQ=255)
               ccGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCTGCGGCGGGTtttttt                                                                                    <  1:2682312/100‑3 (MQ=255)
               ccGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGa                                                                                    <  1:1324881/100‑1 (MQ=255)
                cGTGAGCGGATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCTGCGGCGGGTttttt                                                                                     <  1:1556225/98‑2 (MQ=255)
                cGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAg                                                                                   <  1:605926/100‑1 (MQ=255)
                cGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAg                                                                                   <  2:1235728/100‑1 (MQ=255)
                cGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAg                                                                                   >  2:2081252/1‑100 (MQ=255)
                 gTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCACACAAATGGATATTACAGACAAAAAAACCCGCTGCGGCGGGTttttttgt                                                                                  <  1:2564159/100‑5 (MQ=255)
                    aGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCTGCGGCGGGTttttttgtc                                                                                 <  1:1082001/98‑6 (MQ=38)
                    aGCGCATCATCTGACAAATAAATGGCAGGTGAAGGAGGCATAGGTCATAAAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAg                                                                                   <  1:1572640/96‑1 (MQ=255)
                     gcgcATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCggg                                                                              <  1:2739872/100‑1 (MQ=255)
                      cgcATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGt                                                                             <  2:2205282/100‑1 (MQ=255)
                       gcATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAGAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGtt                                                                            >  1:2767384/1‑100 (MQ=255)
                        catcatCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTc                                                                           <  2:1021909/100‑1 (MQ=255)
                         atcatcTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCg                                                                          >  2:1219489/1‑100 (MQ=255)
                            atcTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGAtt                                                                       <  2:3134306/100‑1 (MQ=255)
                              cTGTGAACTAAATGGCAGGGGAAGGAGGCATAGTTCATACAAATGGATATTACAGACAAAAAACCCCGCCGCAGCGGGGCTTTGAGCCGGGTTCGATTAg                                                                     <  2:1542565/100‑1 (MQ=255)
                               tGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAgt                                                                    <  1:957181/100‑1 (MQ=255)
                                gTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAgtg                                                                   <  1:2916639/100‑1 (MQ=255)
                                 tcaccaaaaTGGAATGTAAAGGAGGCGGAGGTCAGACAAATGGATATTACAGACAAAAAAACCCGTGGCAGCGGGTCTTTCAGCCGGGTTCGATTAgtg                                                                   <  1:1517149/93‑1 (MQ=255)
                                 tCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAgtg                                                                   <  1:276111/99‑1 (MQ=255)
                                  cAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGtt                                                                 >  2:2511984/1‑100 (MQ=255)
                                    aCTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCg                                                               >  2:2370121/1‑100 (MQ=255)
                                     cTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGa                                                              <  1:517624/100‑1 (MQ=255)
                                       aaaTGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAAc                                                            <  1:2238775/100‑1 (MQ=255)
                                          tGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATg                                                         >  1:2885734/1‑100 (MQ=255)
                                          tGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATg                                                         >  2:2910670/1‑100 (MQ=255)
                                           ggaaGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATg                                                         <  1:523890/96‑1 (MQ=255)
                                           ggCAGGTGAAGGAGGCATAGGTCATACAAATGGATATGACAGACAAAAAAACCTGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATgg                                                        >  1:1904417/1‑100 (MQ=255)
                                            gCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGc                                                       <  2:3331176/100‑1 (MQ=255)
                                              aGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCag                                                     <  1:2909289/100‑1 (MQ=255)
                                               ggTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTATTCGAACATGGCaga                                                    >  2:3529625/1‑100 (MQ=255)
                                                 tGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAgaga                                                  <  2:1891042/100‑1 (MQ=255)
                                                   aaGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTAAGCCGGGTTCGATTCGTGTTCGAACATGGCAGAGATc                                                >  1:2395805/1‑100 (MQ=255)
                                                        ggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttc                                           >  2:1659579/1‑100 (MQ=255)
                                                        ggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttc                                           >  1:3126502/1‑100 (MQ=255)
                                                          cATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttctt                                         >  1:2686006/1‑100 (MQ=255)
                                                             aGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGt                                      >  1:1665860/1‑100 (MQ=255)
                                                              ggTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGtt                                     >  2:55153/1‑100 (MQ=255)
                                                              ggTCATACAAATGGAAAAAAAAAACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGtt                                     <  1:1022767/100‑1 (MQ=255)
                                                                tCATACAAATGGAGATTAAAAAAAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGc                                   <  1:1667075/100‑1 (MQ=255)
                                                                 cATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCt                                  >  2:1650180/1‑100 (MQ=255)
                                                                 cATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCt                                  >  2:3204351/1‑100 (MQ=255)
                                                                  aTACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTg                                 <  1:2748345/100‑1 (MQ=255)
                                                                     cAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGAt                               >  1:3497000/1‑99 (MQ=255)
                                                                     cAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATa                              >  1:178100/1‑100 (MQ=255)
                                                                        aTGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacga                           >  2:3451252/1‑100 (MQ=255)
                                                                        aTGGATACTACGGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacga                           >  1:1865935/1‑100 (MQ=255)
                                                                         tGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacgac                          <  2:406650/100‑1 (MQ=255)
                                                                               ttACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATc                     <  2:3220536/99‑1 (MQ=255)
                                                                               ttACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCg                    >  2:486740/1‑100 (MQ=255)
                                                                                   aGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCtt                 <  1:2053791/99‑1 (MQ=255)
                                                                                   aGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTc                <  1:1318881/100‑1 (MQ=255)
                                                                                    ggctcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGGTTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGAGACGACGAATCGCTTCg               >  2:1882716/10‑100 (MQ=37)
                                                                                    gACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCg               >  1:1639297/1‑100 (MQ=255)
                                                                                    gACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGCGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCg               >  1:365783/1‑100 (MQ=255)
                                                                                    gACAAAAAAACACGCCGCAGCGGGTCTTTGAGCCGGGATCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCg               >  1:1856833/1‑100 (MQ=255)
                                                                                    atctcaaagaCCCGCTGCGGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCg               <  1:3249096/91‑1 (MQ=255)
                                                                                     gctcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCgg              >  2:2700506/9‑100 (MQ=38)
                                                                                     gctcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCg               >  1:716318/9‑99 (MQ=38)
                                                                                     gctcaaagaCCCGCTGCGGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCgg              >  2:111716/9‑100 (MQ=255)
                                                                                     aCAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCgg              <  1:320616/100‑1 (MQ=255)
                                                                                     aCAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCgg              <  2:2308077/100‑1 (MQ=255)
                                                                                     aCAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCgg              <  1:1136150/100‑1 (MQ=255)
                                                                                      ctcaaagaCCCGCTGCGGCGGGTCTTTGCGCCGGGTTCGATTAGTGTTCGCACCTGGCCGCGATCGATTCCTCCTTGCTGATACGACGAATCGCTTCg               >  1:2621466/4‑98 (MQ=255)
                                                                                      ctcaaagaCCCGCTGCGGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGc             >  1:2753594/8‑100 (MQ=255)
                                                                                      cAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGc             >  1:2735286/1‑100 (MQ=255)
                                                                                       tcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGcc            >  1:910348/7‑100 (MQ=38)
                                                                                       aaaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATTGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCg               >  1:1813538/1‑97 (MQ=255)
                                                                                       aaaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGcc            <  2:924903/100‑1 (MQ=255)
                                                                                       aaaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGcc            >  1:3464355/1‑100 (MQ=255)
                                                                                        caaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCa           >  2:485742/2‑100 (MQ=255)
                                                                                        caaagaCCCGCTGCGGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCgagat           >  2:3249118/2‑96 (MQ=255)
                                                                                        aaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGGTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGcc            >  1:2407309/1‑99 (MQ=255)
                                                                                          aagaCCCGCTGCGGCGGGTTATTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGc         >  2:1383133/4‑100 (MQ=255)
                                                                                          aagaCCCGCTGCGGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGc         >  2:926029/4‑100 (MQ=255)
                                                                                           aaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCa        <  2:2373623/100‑1 (MQ=255)
                                                                                           aaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCa        >  1:1645980/1‑100 (MQ=255)
                                                                                             aCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATa      <  2:808370/100‑1 (MQ=255)
                                                                                              cccGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATAc     >  2:442022/1‑100 (MQ=255)
                                                                                              cccGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATAc     <  1:515254/100‑1 (MQ=255)
                                                                                              cccGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATAc     >  1:280620/1‑100 (MQ=255)
                                                                                               ccgccgCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATAcc    >  2:2210105/1‑100 (MQ=255)
                                                                                               ccgccgCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATAcc    >  1:331665/1‑100 (MQ=255)
                                                                                                cgccgcAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCt   >  1:1988940/1‑100 (MQ=255)
                                                                                                 gccgcAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTg  <  2:1686544/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
ACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTG  >  NC_000913/1260802‑1260998

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: