Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,620,224 C→G 15.5% intergenic (+17/+14) marB → / ← eamA periplasmic mar operon regulator/cysteine and O‑acetyl‑L‑serine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,620,2240CG15.5% 359.6 / 57.4 136intergenic (+17/+14)marB/eamAperiplasmic mar operon regulator/cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base C (61/54);  new base G (13/8);  total (74/62)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01

Only 100 of 140 total aligned reads displayed.

CGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAA  >  NC_000913/1620125‑1620322
                                                                                                   |                                                                                                  
cGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAc                                                                                                    <  1:2583956/100‑1 (MQ=255)
cGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAc                                                                                                    <  1:2583789/100‑1 (MQ=255)
cGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAc                                                                                                    <  1:2227716/100‑1 (MQ=255)
 gTTTGATTTACATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCACCACGCTATGTAGTTTGTTCTGGCCCCGACa                                                                                                   <  1:2887318/100‑1 (MQ=255)
 gTTTGATGTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCCACACGCTATGTAGTTTGTTCTGGCCCCGAc                                                                                                    <  1:2542987/99‑1 (MQ=255)
     gATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCGGACAttt                                                                                                >  1:3741/1‑97 (MQ=255)
     gATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGAAGCGCTCGGAGTGCCCTAGCATAATCAACACGCTATGTAGTTTGTTGTGGCGCCGACa                                                                                                   >  1:2146007/1‑96 (MQ=255)
      aTTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCACTCGGCGTGCCCTATTATAATCAACACGCTATAGAGTTTGGTATGGCCCCGACATCTCg                                                                                              >  1:466263/1‑100 (MQ=255)
       tttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCgg                                                                                             >  2:278492/1‑100 (MQ=255)
       tttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCgg                                                                                             >  1:1003539/1‑100 (MQ=255)
          aaaTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGc                                                                                          <  2:4378432/100‑1 (MQ=255)
           aaTCACATGGGTACTGGCAGTGATAAGTCGGATTCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGcc                                                                                         <  1:1932580/100‑2 (MQ=255)
            aTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGcca                                                                                        <  2:304143/100‑3 (MQ=255)
            aTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGAGTGCCGTATTATAATCAACACGCTATGTAGTTTGGTCTGGCCCCGACATCTCGGGGc                                                                                          >  1:1885378/1‑98 (MQ=255)
             tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccag                                                                                       <  2:1574579/100‑4 (MQ=255)
             tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccag                                                                                       <  2:2601909/100‑4 (MQ=255)
             tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGc                                                                                          <  2:2851256/97‑1 (MQ=255)
              cacaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttat                                                                                      >  1:1203757/1‑100 (MQ=255)
              cacaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttat                                                                                      <  1:3964326/100‑1 (MQ=255)
              cacaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttat                                                                                      >  2:769539/1‑100 (MQ=255)
               acaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccagaa                                                                                     <  1:502973/100‑6 (MQ=255)
               acaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                     >  2:1895388/1‑100 (MQ=255)
                caTGGGTACTGGCAGTGATAAGTCGGATGCGCTAGGCGTGCCCTATTATCATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttatta                                                                                    >  1:97487/1‑100 (MQ=255)
                 aTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTaa                                                                                   <  2:3851345/100‑1 (MQ=255)
                 aTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTaa                                                                                   >  2:964954/1‑100 (MQ=255)
                  tGGGTACTGGCAGTGATAGGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCCATGTAGGTGGCTCTGGCCCCGACATCTCGGGGCttat                                                                                      >  1:3323216/1‑96 (MQ=255)
                  tGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAAc                                                                                  <  2:4305037/100‑1 (MQ=255)
                   gggTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGcca                                                                                        <  1:172419/93‑3 (MQ=255)
                   gggTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACt                                                                                 >  2:4252126/1‑100 (MQ=255)
                    ggTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccagaacaaac                                                                                <  2:4088902/100‑11 (MQ=38)
                      tACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTcc                                                                              >  1:878548/1‑100 (MQ=255)
                        cTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACtagatc                                                                            >  1:3924382/1‑95 (MQ=255)
                        cTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCa                                                                            <  2:2745905/100‑1 (MQ=255)
                         tGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCAc                                                                           >  1:4096441/1‑100 (MQ=255)
                         tGGCAGGGATAAGTCGGATGCGCTCGGCGGGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGTGGCTTATTAACTTCCCAc                                                                           <  1:1592004/100‑1 (MQ=255)
                          ggCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCAcc                                                                          <  1:3275998/100‑1 (MQ=255)
                            cAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCtt                                                                        >  2:1514148/1‑100 (MQ=255)
                              gTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTa                                                                      >  1:2272804/1‑100 (MQ=255)
                                gATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAcc                                                                    >  1:3891485/1‑100 (MQ=255)
                                gATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAcc                                                                    >  2:1434493/1‑100 (MQ=255)
                                 aTAAGTCGGATGCGCTCGGCGTGCCCTATTATAATAAACACGCTATGTTGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCg                                                                   >  1:2515853/1‑100 (MQ=255)
                                  tAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGc                                                                  <  1:310604/100‑1 (MQ=255)
                                   aaGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCt                                                                 >  2:3296640/1‑100 (MQ=255)
                                    aGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCtt                                                                >  1:893426/1‑100 (MQ=255)
                                      tCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTa                                                              <  2:2352297/100‑1 (MQ=255)
                                       cGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTAGGTCGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTAc                                                             <  1:3090047/100‑1 (MQ=255)
                                        ggATGTGATCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACg                                                            <  1:2235650/100‑1 (MQ=255)
                                        ggATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACg                                                            <  2:2226756/100‑1 (MQ=255)
                                               cTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                     >  2:1679635/1‑100 (MQ=255)
                                               cTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                     <  2:1889680/100‑1 (MQ=255)
                                                  ggCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAg                                                  >  1:512580/1‑100 (MQ=255)
                                                    cGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                                <  1:369663/100‑1 (MQ=255)
                                                      tGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGTCACGACATCTCGGGGCTTATTAACTTTCCACCTTTACCGCTTTACGCCACCGCAAGCCa                                               >  1:2698751/1‑99 (MQ=255)
                                                       gCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCaaa                                             <  2:3499801/100‑1 (MQ=255)
                                                        cccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAAt                                            <  2:1043490/100‑1 (MQ=255)
                                                         ccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                           <  2:3781177/100‑1 (MQ=255)
                                                         ccTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATa                                           <  1:1895382/100‑1 (MQ=255)
                                                         ccTATTATAATCAACACGCTATGTAGCTCGTTCTGGCCACGACATCCCGGGGATTATTACCTTCCCCCCTTTACCCCTCTACGCCACCGCAATCCAAATa                                           >  1:1983281/1‑100 (MQ=255)
                                                            attatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAt                                        >  1:4017347/1‑100 (MQ=255)
                                                            attatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAt                                        >  1:3340696/1‑100 (MQ=255)
                                                             ttatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTCATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAtt                                       <  2:1717087/100‑1 (MQ=255)
                                                              tataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAtt                                       <  2:2484334/99‑1 (MQ=255)
                                                               ataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGa                                     <  1:3522503/100‑1 (MQ=255)
                                                               ataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGa                                     >  1:3109788/1‑100 (MQ=255)
                                                               ataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGa                                     <  2:319385/100‑1 (MQ=255)
                                                               ataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGAAAGCCAAATACATTGa                                     >  2:367681/1‑100 (MQ=255)
                                                                 aaTCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGata                                   <  2:934079/100‑1 (MQ=255)
                                                                  aTCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGatat                                  <  2:2060297/100‑1 (MQ=255)
                                                                    cAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAc                                <  2:2916955/100‑1 (MQ=255)
                                                                     aaCACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACa                               >  2:3484563/1‑100 (MQ=255)
                                                                     aaCACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACa                               >  1:3644327/1‑100 (MQ=255)
                                                                      acacGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAg                              <  1:2835065/100‑1 (MQ=255)
                                                                       cacGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGc                             >  2:1766949/1‑100 (MQ=255)
                                                                        acGCTATGTAGTTTGTTCTGGCTCCGACATCTGGGGGCTTATTAACTTCTTACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACa                               >  1:3327165/1‑97 (MQ=255)
                                                                          gCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCg                          >  2:3879755/1‑100 (MQ=255)
                                                                             aTGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                       <  1:986777/100‑1 (MQ=255)
                                                                               gTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCAt                     <  2:78302/100‑1 (MQ=255)
                                                                                 aGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATaa                   >  2:4282998/1‑100 (MQ=255)
                                                                                   tttGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATg                 >  2:2643131/1‑100 (MQ=255)
                                                                                    gttaataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGa                >  2:1059438/9‑100 (MQ=38)
                                                                                     ttaataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAg               >  2:1176042/8‑100 (MQ=255)
                                                                                     tGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAg               >  1:3243445/1‑100 (MQ=255)
                                                                                      taataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgc              >  2:3404041/7‑100 (MQ=255)
                                                                                      taataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgc              >  2:3499324/7‑100 (MQ=255)
                                                                                      gTTCTGGCCCCGTCATCTCGGCGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCTCGGTCATAATGAgc              >  1:412830/1‑100 (MQ=255)
                                                                                       ttCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgca             >  1:1033750/1‑100 (MQ=255)
                                                                                       aataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgca             >  2:2594894/6‑100 (MQ=255)
                                                                                       aataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgca             >  2:613287/6‑100 (MQ=255)
                                                                                       aataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgca             >  2:3262539/6‑100 (MQ=255)
                                                                                        tCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcac            >  2:1652984/1‑100 (MQ=255)
                                                                                        tCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcac            <  2:3810228/100‑1 (MQ=255)
                                                                                        ataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgc              >  2:2581672/5‑98 (MQ=255)
                                                                                           ggCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgc         >  2:2418928/1‑100 (MQ=255)
                                                                                           ggCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgc         >  2:4128063/1‑100 (MQ=255)
                                                                                           ggCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgc         <  1:3502381/100‑1 (MQ=255)
                                                                                           ggCCCCGACATCTCGGGGCTTATTAACTTACCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgaacagc         >  1:381799/1‑100 (MQ=255)
                                                                                           agCCCCGAGATGTAGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgc         >  2:2984415/2‑100 (MQ=255)
                                                                                            gCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgca        >  2:206876/1‑100 (MQ=255)
                                                                                                cGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTa    >  2:4192698/1‑100 (MQ=255)
                                                                                                  aCATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTaaa  <  1:823308/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
CGTTTGATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAA  >  NC_000913/1620125‑1620322

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