Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,323,291 T→A 8.6% intergenic (+40/‑156) 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,2910TA 8.6% 326.5 / 16.5 93intergenic (+40/‑156)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 (42/43);  new base A (5/3);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCAA  >  NC_000913/2323143‑2323434
                                                                                                                                                    |                                                                                                                                               
ggAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAt                                                                                                                                                 >  2:1311476/1‑149 (MQ=255)
          gaagGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATa                                                                                                                                       >  1:329156/1‑149 (MQ=255)
             ggAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtttt                                                                                                                                              >  4:1210588/1‑139 (MQ=255)
             ggAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtttt                                                                                                                                              <  3:1210588/139‑1 (MQ=255)
             ggAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGcc                                                                                                                                    >  2:354794/1‑149 (MQ=255)
             ggAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGcc                                                                                                                                    >  2:504550/1‑149 (MQ=255)
                aaCCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCt                                                                                                                                 >  4:804232/1‑149 (MQ=255)
                 aCCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGCTGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGcc                                                                                                                                    <  1:1433630/145‑1 (MQ=255)
                 aCCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGcc                                                                                                                                    >  2:1433630/1‑145 (MQ=255)
                 aCCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTc                                                                                                                                >  1:1105369/1‑149 (MQ=255)
                   cGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCAt                                                                                                                              >  4:195823/1‑149 (MQ=255)
                   cGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCAt                                                                                                                              >  1:654072/1‑149 (MQ=255)
                          gCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCATAGTGCGCAATTTTCTGCATAGCCGCTCATCCTCCtt                                                                                                                       >  3:973327/1‑149 (MQ=255)
                           cACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTta                                                                                                                      <  1:504550/149‑1 (MQ=255)
                            aCTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCAt                                                                                                                                        <  1:1407656/130‑1 (MQ=38)
                            aCTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCAt                                                                                                                                        >  2:1407656/1‑130 (MQ=38)
                               gCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATaaa                                                                                                                  <  3:195823/149‑1 (MQ=255)
                                cTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                                                                         >  2:1165565/1‑141 (MQ=255)
                                cTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                                                                         <  1:1165565/141‑1 (MQ=255)
                                 tttAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATc                                                                                                                >  4:377279/1‑149 (MQ=255)
                                  ttAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATcc                                                                                                               <  2:654072/149‑1 (MQ=255)
                                   tAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                                                                              <  2:329156/149‑1 (MQ=255)
                                     aTGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCa                                                                                                            <  3:804232/149‑1 (MQ=255)
                                     aTGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCa                                                                                                            <  1:354794/149‑1 (MQ=255)
                                         tGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCa                                                                                                        <  4:973327/149‑1 (MQ=255)
                                          gggCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCaa                                                                                                       >  4:1451090/1‑149 (MQ=255)
                                                  gTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                                                                         >  3:1105135/1‑123 (MQ=255)
                                                  gTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                                                                         <  4:1105135/123‑1 (MQ=255)
                                                      ccGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                                                                                           <  2:1105369/149‑1 (MQ=255)
                                                         tGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCatat                                                                                        >  4:12650/1‑149 (MQ=255)
                                                                   tataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                                                                         <  3:1371199/106‑1 (MQ=255)
                                                                   tataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                                                                         >  4:1371199/1‑106 (MQ=255)
                                                                   tataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATAtt                                                                                       >  2:715781/1‑140 (MQ=255)
                                                                   tataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATAtt                                                                                       <  1:715781/140‑1 (MQ=255)
                                                                       aaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTc                                                                                            <  2:1291009/131‑1 (MQ=255)
                                                                       aaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTc                                                                                            >  1:1291009/1‑131 (MQ=255)
                                                                         cTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATaa                                                                        <  3:377279/149‑1 (MQ=255)
                                                                          tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCtt                                                                                                                       >  1:506851/1‑101 (MQ=255)
                                                                          tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCtt                                                                                                                       <  2:506851/101‑1 (MQ=255)
                                                                          tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                                                                                   <  3:540056/137‑1 (MQ=255)
                                                                          tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                                                                                   >  4:540056/1‑137 (MQ=255)
                                                                             aGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa                                                                    >  2:990240/1‑149 (MQ=255)
                                                                              ggAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGt                                                                               <  2:420219/137‑1 (MQ=255)
                                                                              ggAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGt                                                                               >  1:420219/1‑137 (MQ=255)
                                                                               gAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt                                                                  <  1:990240/149‑1 (MQ=255)
                                                                                aaTACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATcc                                                                                                         <  2:1267162/109‑1 (MQ=255)
                                                                                aaTACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATcc                                                                                                         >  1:1267162/1‑109 (MQ=255)
                                                                                     ggTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATcc                                                                                                >  2:856966/1‑113 (MQ=255)
                                                                                     ggTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATcc                                                                                                <  1:856966/113‑1 (MQ=255)
                                                                                     ggTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaata                                                            <  3:1451090/149‑1 (MQ=255)
                                                                                      gTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaa                                                              <  1:105861/146‑1 (MQ=255)
                                                                                      gTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaa                                                              >  2:105861/1‑146 (MQ=255)
                                                                                       tATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTAt                                                                                                  <  1:1286710/109‑1 (MQ=255)
                                                                                       tATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTAt                                                                                                  >  2:1286710/1‑109 (MQ=255)
                                                                                       tATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTACAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaataat                                                          <  1:1311476/149‑1 (MQ=255)
                                                                                        aTCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGt                                                                               >  2:278744/1‑127 (MQ=255)
                                                                                        aTCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGt                                                                               <  1:278744/127‑1 (MQ=255)
                                                                                            aTCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATAtt                                                                                       <  3:1151904/115‑1 (MQ=255)
                                                                                            aTCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATAtt                                                                                       >  4:1151904/1‑115 (MQ=255)
                                                                                                ggcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTg                                                 >  3:737545/1‑149 (MQ=255)
                                                                                                   ggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACt                                                                      <  1:592798/125‑1 (MQ=255)
                                                                                                   ggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACt                                                                      >  2:592798/1‑125 (MQ=255)
                                                                                                         ataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTcccc                                        <  4:737545/149‑1 (MQ=255)
                                                                                                           aaCACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCtt                                      >  3:1060381/1‑149 (MQ=255)
                                                                                                              acCAAAACTTGCTATGCAGAAGTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCt                                   <  4:1060381/149‑1 (MQ=255)
                                                                                                              acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTcccc                                        >  3:983500/1‑144 (MQ=255)
                                                                                                              acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTcccc                                        <  4:983500/144‑1 (MQ=255)
                                                                                                               ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGcc                                                                                                                                    <  1:1257485/51‑1 (MQ=255)
                                                                                                               ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGcc                                                                                                                                    >  2:1257485/1‑51 (MQ=255)
                                                                                                                cAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaat                                                             <  2:5170/121‑1 (MQ=255)
                                                                                                                cAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaat                                                             >  1:5170/1‑121 (MQ=255)
                                                                                                                cAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAAtc                                                         >  3:772686/1‑125 (MQ=255)
                                                                                                                cAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAAtc                                                         <  4:772686/125‑1 (MQ=255)
                                                                                                                      ttGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATc                                                                 >  3:551447/1‑111 (MQ=255)
                                                                                                                      ttGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATc                                                                 <  4:551447/111‑1 (MQ=255)
                                                                                                                      ttGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATc                                                                 <  3:1315038/111‑1 (MQ=255)
                                                                                                                      ttGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATc                                                                 >  4:1315038/1‑111 (MQ=255)
                                                                                                                          tATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATAtt                                                                                       <  4:1337448/85‑1 (MQ=21)
                                                                                                                          tATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATAtt                                                                                       >  3:1337448/1‑85 (MQ=25)
                                                                                                                          tATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTgg                                                <  2:947283/124‑1 (MQ=38)
                                                                                                                          tATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTg                       >  3:738204/1‑149 (MQ=39)
                                                                                                                          tATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCAATCTgg                                                >  1:947283/1‑124 (MQ=37)
                                                                                                                              cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAttt                                                                                                    <  3:807521/68‑1 (MQ=255)
                                                                                                                              cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAttt                                                                                                    >  4:807521/1‑68 (MQ=255)
                                                                                                                               aGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATcc                                                                                                <  4:1282463/71‑1 (MQ=255)
                                                                                                                               aGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATcc                                                                                                >  3:1282463/1‑71 (MQ=255)
                                                                                                                                 aaaaTTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCcac                >  2:602244/1‑149 (MQ=255)
                                                                                                                                       gCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGgt          <  3:12650/149‑1 (MQ=255)
                                                                                                                                        cacaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGgtg         >  4:554259/1‑149 (MQ=255)
                                                                                                                                         acaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGc                                    <  1:385947/121‑1 (MQ=255)
                                                                                                                                         acaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGc                                    >  2:385947/1‑121 (MQ=255)
                                                                                                                                         acaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGgtgt        >  1:389329/1‑149 (MQ=255)
                                                                                                                                         acaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGgtgt        <  3:554259/149‑1 (MQ=255)
                                                                                                                                               cgcAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCaa  <  2:389329/149‑1 (MQ=255)
                                                                                                                                                  aaTTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGgt          >  4:1103808/1‑138 (MQ=255)
                                                                                                                                                  aaTTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGgt          <  3:1103808/138‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCACAACGGTGTATGCAA  >  NC_000913/2323143‑2323434

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: