Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,323,275 T→A 16.5% intergenic (+24/‑172) atoC → / → atoD fused response regulator of ato operon, in two‑component system with AtoS: response regulator/sigma54 interaction protein/acetyl‑CoA:acetoacetyl‑CoA transferase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,323,2750TA16.5% 404.9 / 63.4 158intergenic (+24/‑172)atoC/atoDfused response regulator of ato operon, in two‑component system with AtoS: response regulator/sigma54 interaction protein/acetyl‑CoA:acetoacetyl‑CoA transferase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (70/60);  new base A (14/12);  total (84/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

Only 100 of 159 total aligned reads displayed.

TTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCA  >  NC_000913/2323177‑2323373
                                                                                                  |                                                                                                  
ttAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAtt                                                                                                   <  1:3613598/100‑1 (MQ=255)
   aTGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTGCAGGAATACGGTCTTGATACGTCGGATGTATAACACCAACACTTGCTATGCAGAAATTTg                                                                                                 >  1:3736870/1‑99 (MQ=255)
   aTGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGc                                                                                                <  1:492957/100‑1 (MQ=255)
   aTGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGc                                                                                                >  1:3084058/1‑100 (MQ=255)
     gCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGcac                                                                                              <  1:3081766/100‑1 (MQ=255)
       tGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAg                                                                                            <  1:1288418/100‑1 (MQ=255)
         ggCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTg                                                                                          <  1:1963842/100‑1 (MQ=255)
           cATCTGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCTGATGTATAACACCAAAACTTGCTATGCAGAAATTTGAAAAGTgc                                                                                         >  2:778336/1‑99 (MQ=255)
               aGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAt                                                                                    <  1:1483866/100‑1 (MQ=255)
                gTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtt                                                                                   <  1:3973531/100‑1 (MQ=255)
                gTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAt                                                                                    >  2:1285168/1‑99 (MQ=255)
                  cGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtttt                                                                                 >  1:3176030/1‑100 (MQ=255)
                   gCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTc                                                                                >  1:1131258/1‑100 (MQ=255)
                     cGTGCATTGATGTATAAACTCCCGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAGTTTGCACAGTGCGCAATTTTCTg                                                                              >  2:2438339/1‑100 (MQ=255)
                     cGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTg                                                                              <  2:1601458/100‑1 (MQ=37)
                      gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGc                                                                             <  2:2303814/100‑1 (MQ=37)
                      gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGc                                                                             <  1:3449386/100‑1 (MQ=25)
                       tGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCa                                                                            <  1:383894/100‑1 (MQ=255)
                       tGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCa                                                                            >  1:4180619/1‑100 (MQ=255)
                        gCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCAt                                                                           <  2:1839691/100‑1 (MQ=25)
                          aTTTATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTAAAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAg                                                                         <  2:3818330/100‑1 (MQ=21)
                           ttGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGc                                                                        >  2:2675241/1‑100 (MQ=255)
                            tGATGTATAAGCACCAGGAATACGGTATCGATCCGGAGGATGTAGAACACCAAAACCTGCGATGCAGAAAATGGCGCACTGTGCAAATTTCTGCATAg                                                                         <  2:308833/98‑1 (MQ=18)
                            tGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGcc                                                                       <  2:2621457/100‑1 (MQ=255)
                            tGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGca                                                                       <  1:525689/100‑2 (MQ=25)
                             gATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGcaa                                                                      <  2:3534577/100‑3 (MQ=21)
                                  ataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCAt                                                                 <  2:1422197/100‑1 (MQ=255)
                                  ataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTACGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCAt                                                                 >  2:2675274/1‑100 (MQ=255)
                                    aaaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtc                                                               >  1:1650164/1‑100 (MQ=255)
                                          cAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTta                                                         >  1:3571356/1‑100 (MQ=255)
                                           aGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtat                                                        >  1:401342/1‑100 (MQ=255)
                                           aGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtat                                                        >  1:1236450/1‑100 (MQ=255)
                                           aGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtat                                                        <  1:819213/100‑1 (MQ=255)
                                            ggAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtata                                                       <  1:3284673/100‑1 (MQ=255)
                                            ggAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtata                                                       >  1:2172867/1‑100 (MQ=255)
                                              aacaCGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtata                                                       <  1:779913/95‑1 (MQ=255)
                                              aaTACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATaaa                                                     <  1:3489912/100‑1 (MQ=255)
                                               aTACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAAt                                                    <  1:668185/100‑1 (MQ=255)
                                               aTACGGAATAGATCGGGTGGATGTTTCGCACACAAACTTACTATGCAAAATGTTGCACAGGGCGCTATTTTCTGCATAGCCACtaa                                                                  <  2:3903108/86‑3 (MQ=255)
                                                  cGGTATCGATCGGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTATGCATAGCCGCTCATTCTCCTTATAAAt                                                    >  2:609472/1‑97 (MQ=255)
                                                  cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                 >  1:991886/1‑100 (MQ=255)
                                                  cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                 >  1:3719941/1‑100 (MQ=255)
                                                   ggTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATcccc                                                >  2:513944/1‑100 (MQ=255)
                                                    gTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCa                                               >  1:579316/1‑100 (MQ=255)
                                                      aTCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATc                                             >  2:1062599/1‑100 (MQ=255)
                                                       tCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATcc                                            <  1:1253603/100‑1 (MQ=255)
                                                         gATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCAGAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCaa                                          >  2:3884732/1‑100 (MQ=255)
                                                         gATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCaa                                          <  2:991886/100‑1 (MQ=255)
                                                         gATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCaa                                          >  2:2179399/1‑100 (MQ=255)
                                                         gATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCaa                                          >  2:3184798/1‑100 (MQ=255)
                                                          aTCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCAGCTCATTCTCCTTATAAATCCCCATCCAAt                                         >  1:2082482/1‑100 (MQ=255)
                                                            ccGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAttt                                       >  1:651510/1‑100 (MQ=255)
                                                              ggcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTAt                                     >  2:1490384/1‑100 (MQ=255)
                                                              ggcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTAt                                     >  2:2165315/1‑100 (MQ=255)
                                                              ggcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTAt                                     >  1:4163114/1‑100 (MQ=255)
                                                                 ggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATccc                                  <  2:4104453/100‑1 (MQ=255)
                                                                 ggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATccc                                  <  1:1303898/100‑1 (MQ=255)
                                                                  gATGTATAACACCAAAACTTGCTATGCAGAAATTCGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCCATTTATCCCt                                 <  1:2928286/100‑1 (MQ=255)
                                                                     gTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                              <  2:1228390/100‑1 (MQ=255)
                                                                     gTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                              <  1:2923657/100‑1 (MQ=255)
                                                                     gTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                              <  2:1797428/100‑1 (MQ=255)
                                                                      tataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCat                             >  2:3298030/1‑100 (MQ=255)
                                                                      tataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCat                             <  2:3177535/100‑1 (MQ=255)
                                                                      tataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCat                             <  1:2307953/100‑1 (MQ=255)
                                                                          acacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTc                         <  2:2862167/100‑1 (MQ=255)
                                                                           cacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCa                        >  2:741165/1‑100 (MQ=255)
                                                                           cacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCa                        <  2:3597828/100‑1 (MQ=255)
                                                                            acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                       >  1:95942/1‑100 (MQ=255)
                                                                            acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                       >  2:121517/1‑100 (MQ=255)
                                                                            acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                       >  2:1803120/1‑100 (MQ=255)
                                                                             ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCTCATCCAATTTATCCCTTCATATTCAAt                      >  2:2498885/1‑100 (MQ=255)
                                                                             ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                      >  1:367244/1‑100 (MQ=255)
                                                                              cAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAtt                     >  1:849961/1‑100 (MQ=255)
                                                                                 aaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGt                  <  1:883/100‑1 (MQ=255)
                                                                                   cTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTa                >  1:3357113/1‑100 (MQ=255)
                                                                                   cTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTa                <  1:1645593/100‑1 (MQ=255)
                                                                                   cTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTa                >  1:232628/1‑100 (MQ=255)
                                                                                     tGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTaaa              <  2:481636/100‑1 (MQ=255)
                                                                                     ggCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTaaa              >  1:2626974/2‑100 (MQ=25)
                                                                                      gCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAAt             <  1:3842369/100‑1 (MQ=255)
                                                                                       cTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCTCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAg                   <  2:719967/93‑1 (MQ=255)
                                                                                       cTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAAtn            <  2:3307241/100‑2 (MQ=255)
                                                                                       cTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATa            >  1:268106/1‑100 (MQ=25)
                                                                                       cTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATa            >  1:914887/1‑100 (MQ=25)
                                                                                         aTGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAAc          >  1:3930619/1‑100 (MQ=25)
                                                                                          tGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACt         >  1:2827271/1‑100 (MQ=37)
                                                                                          tGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACt         >  2:3068883/1‑100 (MQ=37)
                                                                                           gCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa        >  2:3570780/1‑100 (MQ=255)
                                                                                           gCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa        >  1:204662/1‑100 (MQ=255)
                                                                                           gCAGAAAATTGCGCACTGTGCAAATTTCTGCCTAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa        >  1:1895799/1‑100 (MQ=37)
                                                                                           gCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa        >  1:1086528/1‑100 (MQ=25)
                                                                                            cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa       >  2:1267485/1‑100 (MQ=37)
                                                                                            cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa       >  1:395521/1‑100 (MQ=37)
                                                                                            cAGAAAATTGCGCACTGTGCAAATTTCTACATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa       >  2:860231/1‑100 (MQ=37)
                                                                                              gAAATTTGCACCGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATACAATATATCCCTTCATATTCAATTAGTTAAATAACTAAAt     <  2:437909/100‑1 (MQ=255)
                                                                                              gAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt     >  1:4078157/1‑100 (MQ=255)
                                                                                              gAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt     >  1:3148143/1‑100 (MQ=255)
                                                                                                aaTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATcc   <  1:352031/100‑1 (MQ=255)
                                                                                                aaTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATcc   <  2:2613415/100‑1 (MQ=255)
                                                                                                 aTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCa  <  1:3372681/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
TTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCA  >  NC_000913/2323177‑2323373

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: