Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,620,227 C→G 15.8% intergenic (+20/+11) 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,2270CG15.8% 349.5 / 58.6 133intergenic (+20/+11)marB/eamAperiplasmic mar operon regulator/cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base C (59/53);  new base G (13/8);  total (72/61)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

Only 100 of 135 total aligned reads displayed.

ATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATT  >  NC_000913/1620131‑1620326
                                                                                                |                                                                                                   
 tttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCgg                                                                                                 >  1:1003539/1‑100 (MQ=255)
 tttAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCgg                                                                                                 >  2:278492/1‑100 (MQ=255)
    aaaTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGc                                                                                              <  2:4378432/100‑1 (MQ=255)
     aaTCACATGGGTACTGGCAGTGATAAGTCGGATTCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGcc                                                                                             <  1:1932580/100‑2 (MQ=255)
      aTCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGAGTGCCGTATTATAATCAACACGCTATGTAGTTTGGTCTGGCCCCGACATCTCGGGGc                                                                                              >  1:1885378/1‑98 (MQ=255)
       tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccag                                                                                           <  2:2601909/100‑4 (MQ=255)
       tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccag                                                                                           <  2:1574579/100‑4 (MQ=255)
       tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGc                                                                                              <  2:2851256/97‑1 (MQ=255)
       tCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttn                                                                                           <  1:4223657/100‑2 (MQ=255)
        cacaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttat                                                                                          <  1:3964326/100‑1 (MQ=255)
        cacaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttat                                                                                          >  2:769539/1‑100 (MQ=255)
        cacaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttat                                                                                          >  1:1203757/1‑100 (MQ=255)
         acaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccagaa                                                                                         <  1:502973/100‑6 (MQ=255)
         acaTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttatt                                                                                         >  2:1895388/1‑100 (MQ=255)
          caTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttatta                                                                                        >  1:2924725/1‑100 (MQ=255)
          caTGGGTACTGGCAGTGATAAGTCGGATGCGCTAGGCGTGCCCTATTATCATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCttatta                                                                                        >  1:97487/1‑100 (MQ=255)
           aTGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTaa                                                                                       >  2:964954/1‑100 (MQ=255)
            tGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAAc                                                                                      <  2:4305037/100‑1 (MQ=255)
            tGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAAc                                                                                      <  1:884498/100‑1 (MQ=255)
             gggTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGcca                                                                                            <  1:172419/93‑3 (MQ=255)
             gggTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACt                                                                                     >  2:4252126/1‑100 (MQ=255)
             gatctcTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACt                                                                                     <  2:3924407/95‑1 (MQ=255)
              ggTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGAGATGTCGGGGccagaacaaac                                                                                    <  2:4088902/100‑11 (MQ=38)
               gTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTc                                                                                   >  1:3538061/1‑100 (MQ=255)
                tACTGGCAGTGATACGTCGGATTCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTc                                                                                   <  2:3642043/99‑1 (MQ=255)
                tACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTcc                                                                                  <  2:1129688/100‑1 (MQ=255)
                tACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTcc                                                                                  >  1:878548/1‑100 (MQ=255)
                  cTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACtagatc                                                                                >  1:3924382/1‑95 (MQ=255)
                  cTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCa                                                                                <  1:2497027/100‑1 (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)
                          gATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAcc                                                                        >  2:1002085/1‑100 (MQ=255)
                          gATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTAcc                                                                        >  1:3891485/1‑100 (MQ=255)
                           aTAAGTCGGATGCGCTCGGCGTGCCCTATTATAATAAACACGCTATGTTGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCg                                                                       >  1:2515853/1‑100 (MQ=255)
                            tAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGc                                                                      <  1:310604/100‑1 (MQ=255)
                             aaGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCt                                                                     >  2:3797414/1‑100 (MQ=255)
                                tCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTa                                                                  <  2:2352297/100‑1 (MQ=255)
                                  ggATGTGATCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACg                                                                <  1:2235650/100‑1 (MQ=255)
                                  ggATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACg                                                                <  1:3494635/100‑1 (MQ=255)
                                         cTCGTCGTTCCCTATTATAATCAACACGATATTTAGTTTGTTCTGGCCCCGACATCTCGTGGCTCATTAACTTCCCCCCTTTACCGCTTTACGCCAcc                                                           <  1:2712165/98‑1 (MQ=255)
                                         cTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                         <  2:1889680/100‑1 (MQ=255)
                                         cTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGc                                                         >  2:1679635/1‑100 (MQ=255)
                                            ggCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAg                                                      >  1:512580/1‑100 (MQ=255)
                                             gCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGc                                                     <  2:2320234/100‑1 (MQ=255)
                                              cGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGcc                                                    <  1:2148307/100‑1 (MQ=255)
                                                 gCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCaaa                                                 <  2:3499801/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:3340696/1‑100 (MQ=255)
                                                      attatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAt                                            >  1:4017347/1‑100 (MQ=255)
                                                       ttatAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTCATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAtt                                           <  2:1717087/100‑1 (MQ=255)
                                                        tataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACAtt                                           <  2:2484334/99‑1 (MQ=255)
                                                         ataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGa                                         >  1:3109788/1‑100 (MQ=255)
                                                         ataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGa                                         <  1:3522503/100‑1 (MQ=255)
                                                         ataATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGa                                         <  2:319385/100‑1 (MQ=255)
                                                           aaTCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGata                                       <  2:934079/100‑1 (MQ=255)
                                                            aTCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGatat                                      <  2:2060297/100‑1 (MQ=255)
                                                              cAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAc                                    <  2:2916955/100‑1 (MQ=255)
                                                              cAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATAc                                    <  2:1784631/100‑1 (MQ=255)
                                                               aaCACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACa                                   >  1:3644327/1‑100 (MQ=255)
                                                               aaCACGCTATGTAGTTTGTTCTGGCCCCGACAACTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACa                                   >  2:4176962/1‑100 (MQ=255)
                                                                acacGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAg                                  <  1:2835065/100‑1 (MQ=255)
                                                                 cacGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGc                                 <  1:3309180/100‑1 (MQ=255)
                                                                  acGCTATGTAGTTTGTTCTGGCTCCGACATCTGGGGGCTTATTAACTTCTTACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACa                                   >  1:3327165/1‑97 (MQ=255)
                                                                   cGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGccc                               <  2:1679343/100‑1 (MQ=255)
                                                                    gCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCg                              >  2:3879755/1‑100 (MQ=255)
                                                                       aTGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTc                           <  1:986777/100‑1 (MQ=255)
                                                                         gTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCAt                         <  2:78302/100‑1 (MQ=255)
                                                                            gattGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAAt                      >  2:2991620/3‑100 (MQ=255)
                                                                            gTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTATGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAAt                      <  1:3797390/100‑1 (MQ=255)
                                                                             tttGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATg                     >  2:2643131/1‑100 (MQ=255)
                                                                              ttGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGa                    <  2:3432203/100‑1 (MQ=255)
                                                                              gttaataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGa                    >  2:3265320/9‑100 (MQ=38)
                                                                              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:3262539/6‑100 (MQ=255)
                                                                                 aataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgca                 >  2:2594894/6‑100 (MQ=255)
                                                                                 aataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgca                 >  2:613287/6‑100 (MQ=255)
                                                                                  tCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcac                <  2:3810228/100‑1 (MQ=255)
                                                                                  ataagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgc                  >  2:2581672/5‑98 (MQ=255)
                                                                                    aagCCCCGAGATGTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccg              >  2:1566928/3‑100 (MQ=255)
                                                                                     ggCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgc             >  2:2418928/1‑100 (MQ=255)
                                                                                     ggCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgc             >  2:4300603/1‑100 (MQ=255)
                                                                                     ggCCCCGACATCTCGGGGCTTATTAACTTACCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgaacagc             >  1:381799/1‑100 (MQ=255)
                                                                                     agCCCCGAGATGTAGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgc             >  2:2984415/2‑100 (MQ=255)
                                                                                      gCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgca            >  2:206876/1‑100 (MQ=255)
                                                                                      gCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAgcaccgca            >  1:2575760/1‑100 (MQ=255)
                                                                                          cGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTa        >  2:4192698/1‑100 (MQ=255)
                                                                                            aCATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTaaa      <  1:823308/100‑1 (MQ=255)
                                                                                                ctcGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAAtt  >  2:2611844/1‑100 (MQ=255)
                                                                                                ctcGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAAtt  >  1:1725008/1‑100 (MQ=255)
                                                                                                |                                                                                                   
ATTTAAATCACATGGGTACTGGCAGTGATAAGTCGGATGCGCTCGGCGTGCCCTATTATAATCAACACGCTATGTAGTTTGTTCTGGCCCCGACATCTCGGGGCTTATTAACTTCCCACCTTTACCGCTTTACGCCACCGCAAGCCAAATACATTGATATACAGCCCGGTCATAATGAGCACCGCACCTAAAAATT  >  NC_000913/1620131‑1620326

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: