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.

GATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGG  >  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)
                          cTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCt                                                                   >  2:306413/1‑100 (MQ=255)
                          cTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCt                                                                   <  1:2749692/100‑1 (MQ=255)
                            ttGTCGACCAGTGCAGCGAGGTAATTTTCCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCttt                                                                 <  2:3963754/100‑1 (MQ=255)
                            ttGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCttt                                                                 <  1:3181739/100‑1 (MQ=255)
                              gTCGACCAGGTCATCTAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCttttt                                                               <  1:4102376/100‑1 (MQ=255)
                              gTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCttttt                                                               >  2:2339351/1‑100 (MQ=255)
                              gTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATCTGGCAGACCGGTCATCGGAACAGATTGGATGAGTTGTTCAATATGCTGACGGCttttt                                                               >  2:2262962/1‑100 (MQ=255)
                               tCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCtttttt                                                              <  2:728285/100‑1 (MQ=255)
                                cGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTc                                                             >  2:34916/1‑100 (MQ=255)
                                 gatctggTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTcc                                                            <  2:331427/95‑1 (MQ=255)
                                 gACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTcc                                                            <  2:1472471/100‑1 (MQ=255)
                                 gACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTcc                                                            >  1:3599690/1‑100 (MQ=255)
                                 gACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTcc                                                            <  2:1100189/100‑1 (MQ=255)
                                  aCCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCt                                                           >  2:3720737/1‑100 (MQ=255)
                                   ccAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTg                                                          <  2:3348121/100‑1 (MQ=255)
                                   ccAGGTCATCGAGGTAATTGTGCAGGTTATTGCGACTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTg                                                          >  2:880138/1‑100 (MQ=255)
                                    cAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGt                                                         <  1:2139175/100‑1 (MQ=255)
                                     aGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGtt                                                        <  2:1336366/100‑1 (MQ=255)
                                     aGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGtt                                                        <  1:2674897/100‑1 (MQ=255)
                                     aGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGtt                                                        >  2:667131/1‑100 (MQ=255)
                                     aGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGtt                                                        <  1:3754273/100‑1 (MQ=255)
                                     aGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGtt                                                        >  2:2777352/1‑100 (MQ=255)
                                     aGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGtt                                                        >  1:2451268/1‑100 (MQ=255)
                                     aGGTCATCGAGGTAATCGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGtt                                                        <  2:1255117/100‑1 (MQ=255)
                                      ggTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTccagatc                                                       >  1:331430/1‑95 (MQ=255)
                                       gTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTcc                                                      <  2:2170599/100‑1 (MQ=255)
                                       gTCATCGAGGTAATTGTGCAGGATATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTTCCTGTTc                                                       >  2:1480562/1‑100 (MQ=255)
                                        tCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTcca                                                     <  1:1353926/100‑1 (MQ=255)
                                         cATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccag                                                    >  2:177073/1‑100 (MQ=255)
                                         cATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccag                                                    <  1:1613596/100‑1 (MQ=255)
                                          aTCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagc                                                   <  1:1452896/100‑1 (MQ=255)
                                          aTCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagc                                                   >  1:12919/1‑100 (MQ=255)
                                           tCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcg                                                  <  2:1740580/100‑1 (MQ=255)
                                           tCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcg                                                  >  1:4086174/1‑100 (MQ=255)
                                            cGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgc                                                 >  2:1097436/1‑100 (MQ=255)
                                            cGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgc                                                 >  1:3145971/1‑100 (MQ=255)
                                             gAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgcc                                                <  1:306416/100‑1 (MQ=255)
                                              aGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgcca                                               >  2:2317875/1‑100 (MQ=255)
                                               ggTAATTGTGCAGGTTATTGCGATTTGGCAGACCTGGCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgcc                                                <  2:676538/98‑1 (MQ=255)
                                               ggTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgccag                                              >  2:1690731/1‑100 (MQ=255)
                                               ggTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgccag                                              >  1:22875/1‑100 (MQ=255)
                                                 tAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgccagcg                                            >  2:645251/1‑100 (MQ=255)
                                                 tAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgccagcg                                            >  1:3129114/1‑100 (MQ=255)
                                                    ttGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCgg                                         <  1:2140744/100‑1 (MQ=255)
                                                     tgtgCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGc                                        <  2:3223810/100‑1 (MQ=255)
                                                       tgCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCa                                      <  1:3720764/100‑1 (MQ=255)
                                                        gCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTAACTATGCTGACGGCTTTTTTACTGTTCAAGCCCCAACGCGGCCat                                     >  1:1558828/1‑100 (MQ=255)
                                                          aGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCatat                                   >  2:790470/1‑100 (MQ=255)
                                                           ggTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCatat                                   >  2:792691/1‑99 (MQ=255)
                                                           ggTTATTGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATg                                  >  1:3336034/1‑100 (MQ=255)
                                                            gTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTGTTCCTGTTTCAGCGCCAGCGAGGCCATATgc                                 >  1:1606379/1‑100 (MQ=255)
                                                              tATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctg                               <  2:1241250/100‑1 (MQ=255)
                                                                ttGCGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctggc                             <  1:1612111/100‑1 (MQ=255)
                                                                  gtgATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctggctg                           <  1:2078496/98‑1 (MQ=255)
                                                                   cGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGa                          >  1:1642936/1‑100 (MQ=255)
                                                                   cGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGa                          >  1:1652850/1‑100 (MQ=255)
                                                                   cGATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGa                          >  1:793465/1‑100 (MQ=255)
                                                                    gATTTGGCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGat                         >  1:3518121/1‑100 (MQ=255)
                                                                        tGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCt                     <  1:1947249/100‑1 (MQ=255)
                                                                          gCAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGc                   <  1:1757251/100‑1 (MQ=255)
                                                                           cAGACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGcc                  >  1:243468/1‑100 (MQ=255)
                                                                           cAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGcc                  <  1:896025/100‑1 (MQ=255)
                                                                             gACCGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAc                >  1:2848082/1‑100 (MQ=255)
                                                                              aCCGGTCATCGGATCAGATTGGATGAGCTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAcg               <  1:95953/100‑1 (MQ=255)
                                                                               ccGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAcgc              >  2:3797183/1‑100 (MQ=255)
                                                                               ccGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAcgc              >  1:3338700/1‑100 (MQ=255)
                                                                               ccGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAcgc              >  1:4144188/1‑100 (MQ=255)
                                                                                cGGTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAcgcg             <  1:2187866/100‑1 (MQ=255)
                                                                                 ggTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCTAGCGCCAGCGAGGCAATATGCTGGCTGATATCTGCAAcgcg             >  2:944804/1‑99 (MQ=255)
                                                                                 ggTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCCGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCAcgcg             >  2:3363934/1‑99 (MQ=255)
                                                                                 ggTCATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCAGCCACGCGt            >  1:1891363/1‑100 (MQ=255)
                                                                                  gTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGtt           <  2:2451253/100‑1 (MQ=255)
                                                                                  gTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCAGATATCAGCCagcatat           <  2:2820687/100‑7 (MQ=255)
                                                                                    cATCGGATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCg         <  1:1659678/100‑1 (MQ=255)
                                                                                    cATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCg         <  1:2265067/100‑1 (MQ=255)
                                                                                        ggATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCACAGCGGCCATATGCTGGCTGATATCGCCCACGCGTTCGAt       <  2:2573467/98‑1 (MQ=255)
                                                                                        ggATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATaa     <  1:1311707/100‑1 (MQ=255)
                                                                                         gATCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTTCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATaaa    >  2:1153612/1‑100 (MQ=255)
                                                                                           tCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAgg  <  2:2626213/100‑1 (MQ=255)
                                                                                           tCAGATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAgg  <  2:1040062/100‑1 (MQ=255)
                                                                                            cAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTTTTTTTTTTTTTCCAACGAGGCCATATACTGGCTGATATCTGCCACGGGGTCGATaaa    >  2:3577532/1‑97 (MQ=255)
                                                                                              |                                                                                                 
GATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGG  >  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: