Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,564,601 A→G 72.4% F378S (TTT→TCT)  dosP ← oxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,564,6010AG72.4% 177.8 / 69.0 108F378S (TTT→TCT) dosPoxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator
Reads supporting (aligned to +/- strand):  ref base A (13/17);  new base G (44/34);  total (57/51)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 108 total aligned reads displayed.

GATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGG  >  NC_000913/1564507‑1564697
                                                                                              |                                                                                                 
gATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTgg                                                                                              >  2:2492067/1‑100 (MQ=255)
 aTGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGa                                                                                             <  1:2774183/100‑1 (MQ=255)
 aTGAAATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCAAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGa                                                                                             <  1:3690781/100‑1 (MQ=255)
  tGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGAt                                                                                            <  2:3319832/100‑1 (MQ=255)
     gaTACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAg                                                                                         >  2:2378111/1‑100 (MQ=255)
     gaTACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAg                                                                                         >  1:3083075/1‑100 (MQ=255)
       tACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAgtt                                                                                       >  1:2330544/1‑100 (MQ=255)
       tACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAgtt                                                                                       <  2:567830/100‑1 (MQ=255)
         caccacGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAgttgt                                                                                     >  1:2385527/1‑100 (MQ=255)
         caccacGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAgttgt                                                                                     >  1:184174/1‑100 (MQ=255)
           ccacGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTc                                                                                   <  2:2198682/100‑1 (MQ=255)
            cacGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCa                                                                                  >  1:1071589/1‑100 (MQ=255)
            cacGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCa                                                                                  >  1:1051169/1‑100 (MQ=255)
              cGGGTGAGACTGCTTTTTCGACAAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAat                                                                                >  2:404493/1‑100 (MQ=255)
              cGGGAGAGACGGGTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAat                                                                                >  1:2744345/1‑100 (MQ=255)
              cGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAat                                                                                <  1:3042332/100‑1 (MQ=255)
               gggAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAata                                                                               <  1:2917001/100‑1 (MQ=255)
                ggAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAatat                                                                              >  2:2375121/1‑100 (MQ=255)
                     gaCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGa                                                                         >  1:2903718/1‑100 (MQ=255)
                          cTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑t                                                                   >  2:306413/1‑100 (MQ=255)
                          cTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑t                                                                   <  1:2749692/100‑1 (MQ=255)
                            ttGTCGACCAGTGCAGCGAGGTAATTTTCCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑ttt                                                                 <  2:3963754/100‑1 (MQ=255)
                            ttGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑ttt                                                                 <  1:3181739/100‑1 (MQ=255)
                              gTCGACCAGGTCATCTAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑ttttt                                                               <  1:4102376/100‑1 (MQ=255)
                              gTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑ttttt                                                               >  2:2339351/1‑100 (MQ=255)
                              gTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATCTGGCAGACCGGTCATCGGAACAGATTGGATGAGTTGTTCAATATGCTGACGGC‑ttttt                                                               >  2:2262962/1‑100 (MQ=255)
                               tCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑tttttt                                                              <  2:728285/100‑1 (MQ=255)
                                cGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTc                                                             >  2:34916/1‑100 (MQ=255)
                                 gatctggTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTcc                                                            <  2:331427/95‑1 (MQ=255)
                                 gACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTcc                                                            <  2:1472471/100‑1 (MQ=255)
                                 gACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTcc                                                            >  1:3599690/1‑100 (MQ=255)
                                 gACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTcc                                                            <  2:1100189/100‑1 (MQ=255)
                                  aCCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCt                                                           >  2:3720737/1‑100 (MQ=255)
                                   ccAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTg                                                          <  2:3348121/100‑1 (MQ=255)
                                   ccAGGTCATCGAGGTAATTGTGCAGGTTATTGCGACTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTg                                                          >  2:880138/1‑100 (MQ=255)
                                    cAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGt                                                         <  1:2139175/100‑1 (MQ=255)
                                     aGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGtt                                                        <  2:1336366/100‑1 (MQ=255)
                                     aGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGtt                                                        <  1:2674897/100‑1 (MQ=255)
                                     aGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGtt                                                        >  2:667131/1‑100 (MQ=255)
                                     aGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGtt                                                        <  1:3754273/100‑1 (MQ=255)
                                     aGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGtt                                                        >  2:2777352/1‑100 (MQ=255)
                                     aGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGtt                                                        >  1:2451268/1‑100 (MQ=255)
                                     aGGTCATCGAGGTAATCGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGtt                                                        <  2:1255117/100‑1 (MQ=255)
                                      ggTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTccagatc                                                       >  1:331430/1‑95 (MQ=255)
                                       gTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTcc                                                      <  2:2170599/100‑1 (MQ=255)
                                       gTCATCGAGGTAATTGTGCAGGATATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTTCCTGTTc                                                       >  2:1480562/1‑100 (MQ=255)
                                        tCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTcca                                                     <  1:1353926/100‑1 (MQ=255)
                                         cATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTccag                                                    >  2:177073/1‑100 (MQ=255)
                                         cATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTccag                                                    <  1:1613596/100‑1 (MQ=255)
                                          aTCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTccagc                                                   <  1:1452896/100‑1 (MQ=255)
                                          aTCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTccagc                                                   >  1:12919/1‑100 (MQ=255)
                                           tCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTccagcg                                                  <  2:1740580/100‑1 (MQ=255)
                                           tCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTccagcg                                                  >  1:4086174/1‑100 (MQ=255)
                                            cGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTccagcgc                                                 >  2:1097436/1‑100 (MQ=255)
                                            cGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTccagcgc                                                 >  1:3145971/1‑100 (MQ=255)
                                             gAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTccagcgcc                                                <  1:306416/100‑1 (MQ=255)
                                              aGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTccagcgcca                                               >  2:2317875/1‑100 (MQ=255)
                                               ggTAATTGTGCAGGTTATTGCGATTTGGCAGACCTGGCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTccagcgcc                                                <  2:676538/98‑1 (MQ=255)
                                               ggTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTccagcgccag                                              >  2:1690731/1‑100 (MQ=255)
                                               ggTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTccagcgccag                                              >  1:22875/1‑100 (MQ=255)
                                                 tAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTccagcgccagcg                                            >  2:645251/1‑100 (MQ=255)
                                                 tAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTccagcgccagcg                                            >  1:3129114/1‑100 (MQ=255)
                                                    ttGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCgg                                         <  1:2140744/100‑1 (MQ=255)
                                                     tgtgCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGc                                        <  2:3223810/100‑1 (MQ=255)
                                                       tgCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCa                                      <  1:3720764/100‑1 (MQ=255)
                                                        gCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTAACTATGCTGACGGC‑TTTTTTACTGTTCAAGCCCCAACGCGGCCat                                     >  1:1558828/1‑100 (MQ=255)
                                                          aGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCatat                                   >  2:790470/1‑100 (MQ=255)
                                                           ggTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCatat                                   >  2:792691/1‑99 (MQ=255)
                                                           ggTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATg                                  >  1:3336034/1‑100 (MQ=255)
                                                            gTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTGTTCCTGTTTCAGCGCCAGCGAGGCCATATgc                                 >  1:1606379/1‑100 (MQ=255)
                                                              tATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctg                               <  2:1241250/100‑1 (MQ=255)
                                                                ttGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctggc                             <  1:1612111/100‑1 (MQ=255)
                                                                  gtgATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctggctg                           <  1:2078496/98‑1 (MQ=255)
                                                                   cGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGa                          >  1:1642936/1‑100 (MQ=255)
                                                                   cGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGa                          >  1:1652850/1‑100 (MQ=255)
                                                                   cGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGa                          >  1:793465/1‑100 (MQ=255)
                                                                    gATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGat                         >  1:3518121/1‑100 (MQ=255)
                                                                        tGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCt                     <  1:1947249/100‑1 (MQ=255)
                                                                          gCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGc                   <  1:1757251/100‑1 (MQ=255)
                                                                           cAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGcc                  >  1:243468/1‑100 (MQ=255)
                                                                           cAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGcc                  <  1:896025/100‑1 (MQ=255)
                                                                             gACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAc                >  1:2848082/1‑100 (MQ=255)
                                                                              aCCGGTCATCGGATCAGATTGGATGAGCTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAcg               <  1:95953/100‑1 (MQ=255)
                                                                               ccGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAcgc              >  2:3797183/1‑100 (MQ=255)
                                                                               ccGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAcgc              >  1:3338700/1‑100 (MQ=255)
                                                                               ccGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAcgc              >  1:4144188/1‑100 (MQ=255)
                                                                                cGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAcgcg             <  1:2187866/100‑1 (MQ=255)
                                                                                 ggTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCTAGCGCCAGCGAGGCAATATGCTGGCTGATATCTGCAAcgcg             >  2:944804/1‑99 (MQ=255)
                                                                                 ggTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCCGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAcgcg             >  2:3363934/1‑99 (MQ=255)
                                                                                 ggTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCAGCCACGCGt            >  1:1891363/1‑100 (MQ=255)
                                                                                  gTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGtt           <  2:2451253/100‑1 (MQ=255)
                                                                                  gTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCAGATATCAGCCagcatat           <  2:2820687/100‑7 (MQ=255)
                                                                                    cATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCg         <  1:1659678/100‑1 (MQ=255)
                                                                                    cATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCg         <  1:2265067/100‑1 (MQ=255)
                                                                                        ggATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCACAGCGGCCATATGCTGGCTGATATCGCCCACGCGTTCGAt       <  2:2573467/98‑1 (MQ=255)
                                                                                        ggATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATaa     <  1:1311707/100‑1 (MQ=255)
                                                                                         gATCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTTCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATaaa    >  2:1153612/1‑100 (MQ=255)
                                                                                           tCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAgg  <  2:2626213/100‑1 (MQ=255)
                                                                                           tCAGATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAgg  <  2:1040062/100‑1 (MQ=255)
                                                                                            cAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTTTTTTTTTTTTTCCAACGAGGCCATATACTGGCTGATATCTGCCACGGGGTCGATaaa    >  2:3577532/1‑97 (MQ=255)
                                                                                              |                                                                                                 
GATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGC‑TTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGG  >  NC_000913/1564507‑1564697

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: ‑