Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,323,286 C→T 15.2% intergenic (+35/‑161) 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,2860CT15.2% 237.5 / 27.8 102intergenic (+35/‑161)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 C (35/51);  new base T (5/11);  total (40/62)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e-01

Only 100 of 106 total aligned reads displayed.

ATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTC  >  NC_000913/2323189‑2323385
                                                                                                 |                                                                                                   
aTCAGTCTCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgcgc                                                                                                   <  1:293634/100‑1 (MQ=255)
aTCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgcgc                                                                                                   <  2:2612253/100‑1 (MQ=255)
 tcggTCGCCGTGCATTGATGTATAAAGTTCAGGAATACTGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCa                                                                                                  <  2:1704417/97‑1 (MQ=255)
 tCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCa                                                                                                  <  2:2871378/100‑1 (MQ=255)
       gCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtttt                                                                                             >  1:2215945/1‑99 (MQ=255)
         cGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTg                                                                                          >  1:1034098/1‑100 (MQ=255)
         cGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTg                                                                                          <  2:2677884/100‑1 (MQ=255)
         cGTGCATTGATGTATAAACTCCAGGAATACGGTATAGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTg                                                                                          <  2:441079/100‑1 (MQ=255)
          gTGCATTGATGTATAAACTCCATGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGc                                                                                         <  1:1080956/100‑1 (MQ=255)
          gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCt                                                                                           <  2:825413/98‑1 (MQ=255)
          gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTg                                                                                          <  2:2634314/99‑1 (MQ=255)
           tGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCa                                                                                        <  2:2080608/100‑1 (MQ=25)
           tGCATTGATGTATAAACTCCAGGAAGACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCa                                                                                        <  2:1098637/100‑1 (MQ=37)
            gCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCAt                                                                                       <  1:2340975/100‑1 (MQ=25)
             cATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCAt                                                                                       <  2:2369145/99‑1 (MQ=25)
              aTTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATa                                                                                      >  2:2352874/1‑99 (MQ=255)
              aTTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAg                                                                                     <  1:141208/100‑1 (MQ=255)
              aTTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAg                                                                                     <  1:144900/100‑1 (MQ=255)
              aTTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAg                                                                                     <  1:2013753/100‑1 (MQ=25)
              aTTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAg                                                                                     <  2:1912905/100‑1 (MQ=25)
              aTTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACCTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAg                                                                                     <  2:1651289/100‑1 (MQ=25)
                 gATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAg                                                                                     >  2:3047598/1‑97 (MQ=255)
                 gATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCg                                                                                  >  1:237917/1‑100 (MQ=255)
                  aTGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGc                                                                                 >  1:1906419/1‑100 (MQ=255)
                  aTGTATAAACTCAAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGcaag                                                                                 <  2:1471251/100‑4 (MQ=21)
                   tGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGc                                                                                 >  2:1917504/1‑99 (MQ=255)
                    gTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTAACACAATACGCAATTTTCTGCATAGCCGCt                                                                                >  1:1642753/1‑99 (MQ=255)
                     tataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCa                                                                              >  2:1207782/1‑100 (MQ=255)
                      ntaAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCAt                                                                             >  2:2139069/2‑100 (MQ=255)
                         aaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAATACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtct                                                                          <  1:26593/100‑1 (MQ=255)
                          aCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtctc                                                                         >  2:2338639/1‑100 (MQ=255)
                            tCCAGGAATACGGTATCGATCCGGCGGATGTATAACTCCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                                       <  1:628607/100‑1 (MQ=255)
                            tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                                       <  2:1663798/100‑1 (MQ=255)
                            tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                                       <  2:2344249/100‑1 (MQ=255)
                            tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                                       >  1:1642867/1‑100 (MQ=255)
                            tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                                       <  1:1917527/100‑1 (MQ=255)
                             ccAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCtt                                                                      >  2:2956243/1‑100 (MQ=255)
                             ccAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCtt                                                                      >  1:2176614/1‑100 (MQ=255)
                             ccAGGAATACGGTATCGATCAGGCGGATGTATAACACCAAAACTTCCTATGCAGAAATTTGCACAGTGCGCAAGTTTCTGCATAGCCGCTCATTCTCCtt                                                                      >  2:848996/1‑100 (MQ=255)
                               aGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtat                                                                    <  1:42355/100‑1 (MQ=255)
                                  aaTACGGTATCGATCCGGCGGATGTATAACACCAAATTTTTCTATGTAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATaaa                                                                 <  1:824859/100‑1 (MQ=255)
                                  aaTACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATaaa                                                                 >  2:2957435/1‑100 (MQ=255)
                                    tACGGTATCGATCCTGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTTCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTTTAAATc                                                               >  2:26822/1‑100 (MQ=255)
                                      cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                             <  2:252616/100‑1 (MQ=255)
                                      cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                             <  1:1467670/100‑1 (MQ=255)
                                      cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                             <  1:399998/100‑1 (MQ=255)
                                      cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                             <  2:2506656/100‑1 (MQ=255)
                                       ggTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATcccc                                                            <  2:214827/100‑1 (MQ=255)
                                        gTATCGATCCGGCGGATGTATAACACCAAAGCTTGCTATGCAGAAATTTGCACAATGCGCTATTTTCTGCATAGCCGCTCATTCTACTTATAATTCCCCa                                                           <  2:468081/100‑1 (MQ=255)
                                        gTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCa                                                           >  1:1903517/1‑100 (MQ=255)
                                           tCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATcc                                                        <  1:809396/100‑1 (MQ=255)
                                             gATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCACCATCCaa                                                      >  1:778025/1‑100 (MQ=255)
                                              aTCCGGCGGATGTATAACACCAAAACTTGCTATGCAGTAATTTGCACAGTGCGCAATTTTATAAATAAACGCTCATTATCCTTATAAATCCCCATCCAAt                                                     >  1:915809/1‑100 (MQ=255)
                                                ccGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAttt                                                   <  2:349233/100‑1 (MQ=255)
                                                 cggcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAttt                                                   <  2:2974765/99‑1 (MQ=255)
                                                   gcggATGTATAACACCAAATCTTGCTATGCAGAAATTTGCACAGTGTGCAATTTTCTGAATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATc                                                <  1:1041842/100‑1 (MQ=255)
                                                    cggATGTATAACACCAAAACTTGCTATACAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATcc                                               >  2:586565/1‑100 (MQ=255)
                                                         gTATAACGCCAAAACTTGCTATGCAGATATTTGCACAGTGCGCAATTGTCTGCATAGCCGCTCATTCTCCTTATAAATCCTCATCCAATTTATCCCTTc                                           >  2:204885/1‑99 (MQ=255)
                                                         gTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                                          <  2:2721458/100‑1 (MQ=255)
                                                         gTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                                          >  1:1657182/1‑100 (MQ=255)
                                                          tataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCat                                         >  1:2275664/1‑100 (MQ=255)
                                                           ataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCACTTTTCTGCATAGCCGCTCATTATCCTTATAAATCCACATCCAATATATACATTCata                                        >  1:2817866/1‑100 (MQ=255)
                                                            taACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCatat                                       >  2:468660/1‑100 (MQ=255)
                                                            taACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCatat                                       <  1:2667593/100‑1 (MQ=255)
                                                             aaCACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATAtt                                      >  1:2401009/1‑100 (MQ=255)
                                                              acacCAAACCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTATCCTTATAACTCCCCATCCAATTTATCCCTTCATGTTc                                     >  2:638703/1‑100 (MQ=255)
                                                               cacCAAAACTTGCTATGCAGAAATTTGCACAGTGTGCAATTTTCTGCATAGTCGCTCATTCTCCTTATAAATCCCCTTACAATGTATCCCTTCatatcc                                     <  2:2058890/99‑3 (MQ=255)
                                                                acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                                   >  2:87189/1‑100 (MQ=255)
                                                                acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                                   >  2:2748108/1‑100 (MQ=255)
                                                                  cAAAACTTGCTATTCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                                  <  1:283373/99‑1 (MQ=255)
                                                                   aaaaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTa                                <  1:2748140/100‑1 (MQ=255)
                                                                      aCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGtt                             <  1:2650199/100‑1 (MQ=255)
                                                                       cTTGCTATGTAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTa                            <  1:1207793/100‑1 (MQ=255)
                                                                       cTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTa                            <  1:458028/100‑1 (MQ=255)
                                                                       cTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTa                            <  1:2637582/100‑1 (MQ=255)
                                                                         tGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTaaa                          >  1:1489114/1‑100 (MQ=255)
                                                                          gCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATACCCATACAATTTATCAATTAATATTAAATTAGTTAAAt                         >  1:342683/1‑100 (MQ=12)
                                                                           cTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATa                        <  2:1015369/100‑1 (MQ=255)
                                                                           cTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATa                        <  1:2338666/100‑1 (MQ=255)
                                                                            tATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCTCTTCATATTCAATTAGTTaaaga                        <  2:2882466/99‑3 (MQ=255)
                                                                             aTGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAAc                      <  2:346739/100‑1 (MQ=255)
                                                                              tGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACt                     <  2:2488680/100‑1 (MQ=255)
                                                                               gCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa                    <  2:1411752/100‑1 (MQ=255)
                                                                               gCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa                    >  2:1775173/1‑100 (MQ=255)
                                                                               gCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa                    >  2:439140/1‑100 (MQ=25)
                                                                                cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa                   >  1:2971196/1‑100 (MQ=37)
                                                                                cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa                   >  1:36331/1‑100 (MQ=37)
                                                                                 aGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCATATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaaa                  >  2:1678504/1‑100 (MQ=25)
                                                                                    aaTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATcc               <  1:2221937/100‑1 (MQ=255)
                                                                                    aaTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATcc               <  1:1771071/100‑1 (MQ=255)
                                                                                      tttGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaa             <  2:2090857/100‑1 (MQ=255)
                                                                                       ttGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAATTCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaat            <  1:2584541/100‑1 (MQ=255)
                                                                                       ttGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaat            <  2:917222/100‑1 (MQ=255)
                                                                                       ttGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaat            <  1:2853547/100‑1 (MQ=255)
                                                                                        tGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaata           >  1:3036324/1‑100 (MQ=255)
                                                                                         gCACAGTGCGCAATTTTCTGCATAGCCGCTCATTTTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaataa          <  2:2980796/100‑1 (MQ=255)
                                                                                          cacaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCACTTCATATTCAATTAGTTAAATAACTAAATCCaataat         >  1:872790/1‑100 (MQ=255)
                                                                                           acaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAAtc        <  1:112551/100‑1 (MQ=255)
                                                                                                gcgcAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCAtt   >  2:1233802/1‑100 (MQ=255)
                                                                                                 cgcAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTc  <  2:911159/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
ATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTC  >  NC_000913/2323189‑2323385

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: