Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,323,275 T→A 7.9% 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,2750TA 7.9% 313.2 / 13.6 88intergenic (+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 (41/40);  new base A (4/3);  total (45/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

CATTATGGAAGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCAC  >  NC_000913/2323128‑2323420
                                                                                                                                                   |                                                                                                                                                 
catTATGGAAGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAtt                                                                                                                                                  <  1:1063844/149‑1 (MQ=255)
               ggAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAt                                                                                                                                   >  2:1311476/1‑149 (MQ=255)
                         gaagGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATa                                                                                                                         >  1:329156/1‑149 (MQ=255)
                            ggAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtttt                                                                                                                                <  3:1210588/139‑1 (MQ=255)
                            ggAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtttt                                                                                                                                >  4:1210588/1‑139 (MQ=255)
                            ggAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGcc                                                                                                                      >  2:504550/1‑149 (MQ=255)
                            ggAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGcc                                                                                                                      >  2:354794/1‑149 (MQ=255)
                              aaaCCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAttt                                                                                                                                                 >  4:74257/1‑120 (MQ=255)
                              aaaCCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAttt                                                                                                                                                 <  3:74257/120‑1 (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                                                                                                                          >  2:1407656/1‑130 (MQ=38)
                                           aCTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCAt                                                                                                                          <  1:1407656/130‑1 (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)
                                                        tGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTg                                                                                                                                                >  2:743224/1‑95 (MQ=255)
                                                        tGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTg                                                                                                                                                <  1:743224/95‑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                                                                         <  1:715781/140‑1 (MQ=255)
                                                                                  tataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATAtt                                                                         >  2:715781/1‑140 (MQ=255)
                                                                                      aaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTc                                                                              >  1:1291009/1‑131 (MQ=255)
                                                                                      aaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTc                                                                              <  2:1291009/131‑1 (MQ=255)
                                                                                        cTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATaa                                                          <  3:377279/149‑1 (MQ=255)
                                                                                         tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCtt                                                                                                         <  2:506851/101‑1 (MQ=255)
                                                                                         tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCtt                                                                                                         >  1:506851/1‑101 (MQ=255)
                                                                                         tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                                                                     >  4:540056/1‑137 (MQ=255)
                                                                                         tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                                                                     <  3:540056/137‑1 (MQ=255)
                                                                                            aGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa                                                      >  2:990240/1‑149 (MQ=255)
                                                                                             ggAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGt                                                                 >  1:420219/1‑137 (MQ=255)
                                                                                             ggAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGt                                                                 <  2:420219/137‑1 (MQ=255)
                                                                                              gAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt                                                    <  1:990240/149‑1 (MQ=255)
                                                                                               aaTACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATcc                                                                                           <  2:1267162/109‑1 (MQ=255)
                                                                                               aaTACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATcc                                                                                           >  1:1267162/1‑109 (MQ=255)
                                                                                                    ggTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATcc                                                                                  <  1:856966/113‑1 (MQ=255)
                                                                                                    ggTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATcc                                                                                  >  2:856966/1‑113 (MQ=255)
                                                                                                    ggTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaata                                              <  3:1451090/149‑1 (MQ=255)
                                                                                                     gTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaa                                                >  2:105861/1‑146 (MQ=255)
                                                                                                     gTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaa                                                <  1:105861/146‑1 (MQ=255)
                                                                                                      tATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTAt                                                                                    >  2:1286710/1‑109 (MQ=255)
                                                                                                      tATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTAt                                                                                    <  1:1286710/109‑1 (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                          <  4:983500/144‑1 (MQ=255)
                                                                                                                             acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTcccc                          >  3:983500/1‑144 (MQ=255)
                                                                                                                              ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGcc                                                                                                                      >  2:1257485/1‑51 (MQ=255)
                                                                                                                              ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGcc                                                                                                                      <  1:1257485/51‑1 (MQ=255)
                                                                                                                               cAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaat                                               >  1:5170/1‑121 (MQ=255)
                                                                                                                               cAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaat                                               <  2:5170/121‑1 (MQ=255)
                                                                                                                               cAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAAtc                                           >  3:772686/1‑125 (MQ=255)
                                                                                                                               cAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAAtc                                           <  4:772686/125‑1 (MQ=255)
                                                                                                                                     ttGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATc                                                   >  4:1315038/1‑111 (MQ=255)
                                                                                                                                     ttGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATc                                                   <  3:1315038/111‑1 (MQ=255)
                                                                                                                                     ttGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATc                                                   <  4:551447/111‑1 (MQ=255)
                                                                                                                                     ttGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATc                                                   >  3:551447/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                                                                                  >  3:1282463/1‑71 (MQ=255)
                                                                                                                                              aGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATcc                                                                                  <  4:1282463/71‑1 (MQ=255)
                                                                                                                                                aaaaTTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCcac  >  2:602244/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CATTATGGAAGTGCTGGAACAACAAGAAGGAAACCGAACCCGCACTGCTTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCACTCCCCTTGCTATTGCCTGACTGTACCCAC  >  NC_000913/2323128‑2323420

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: