Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,323,280 A→G 18.3% intergenic (+29/‑167) 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,2800AG18.3% 398.0 / 79.7 183intergenic (+29/‑167)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 A (81/68);  new base G (16/17);  total (97/86)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

Only 100 of 188 total aligned reads displayed.

TGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATC  >  NC_000913/2323181‑2323379
                                                                                                   |                                                                                                   
  cTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGcaca                                                                                                   >  2:4315235/1‑100 (MQ=255)
      gCATCAGTCGCCGTGCATTGATGTATAATCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgc                                                                                               <  1:1296358/100‑1 (MQ=255)
        aTCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgcgc                                                                                             <  2:2433052/100‑1 (MQ=255)
         tCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCa                                                                                            <  2:1361696/100‑1 (MQ=255)
          cAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCaa                                                                                           >  2:1412838/1‑100 (MQ=255)
               gCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTc                                                                                      <  2:2949877/100‑1 (MQ=255)
                ccGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCt                                                                                     <  1:1676161/100‑1 (MQ=255)
                ccGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGAGGTATAACACCAAAACTTGCTATGCAGAAATTTGCACGGTGCGCAAtt                                                                                         >  1:4408714/1‑96 (MQ=255)
                 cGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTg                                                                                    >  2:41234/1‑100 (MQ=255)
                 cGTGCATTGATGTATAAACTCCAGGAAGACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTg                                                                                    <  1:2328874/100‑1 (MQ=255)
                  gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGc                                                                                   <  1:1695458/100‑1 (MQ=37)
                  gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGc                                                                                   <  2:1440018/100‑1 (MQ=25)
                  gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGc                                                                                   <  1:1665562/100‑1 (MQ=37)
                   tGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCa                                                                                  >  2:1095153/1‑100 (MQ=255)
                   tGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCa                                                                                  <  2:4021385/100‑1 (MQ=255)
                    gCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCAt                                                                                 <  1:2844824/100‑1 (MQ=25)
                     cATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATa                                                                                <  1:493473/100‑1 (MQ=255)
                      aTTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAg                                                                               <  2:1514298/100‑1 (MQ=25)
                      aTTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAg                                                                               <  2:4227855/100‑1 (MQ=25)
                       ttGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGc                                                                              >  1:4187230/1‑100 (MQ=255)
                       ttGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGc                                                                              <  2:3362274/100‑1 (MQ=25)
                        tGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAACTTGCACAGTGCGCAATTTTCTGCATAGcc                                                                             <  2:2565345/100‑1 (MQ=255)
                         gATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCg                                                                            >  2:2489065/1‑100 (MQ=255)
                           tGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCt                                                                          <  2:190102/100‑1 (MQ=255)
                           tGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCt                                                                          <  1:1038279/100‑1 (MQ=255)
                            gTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTc                                                                         >  2:1949752/1‑100 (MQ=255)
                              ataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCAt                                                                       >  2:4403319/1‑100 (MQ=255)
                                aaaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtc                                                                     >  1:222557/1‑100 (MQ=255)
                                aaaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtc                                                                     <  1:1953477/100‑1 (MQ=255)
                                 aaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtct                                                                    <  2:971847/100‑1 (MQ=255)
                                 aaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtct                                                                    >  1:1160077/1‑100 (MQ=255)
                                   cTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                  <  1:805518/100‑1 (MQ=255)
                                    tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                                 <  2:2944017/100‑1 (MQ=255)
                                      cAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTTCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                                 <  1:2714760/98‑1 (MQ=255)
                                      cAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTta                                                               <  1:2884193/100‑1 (MQ=255)
                                        ggAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtata                                                             <  2:20263/100‑1 (MQ=255)
                                        ggAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtata                                                             >  2:681179/1‑100 (MQ=255)
                                        ggAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtata                                                             >  2:3022802/1‑100 (MQ=255)
                                         gAATACGGTATCGATCCAGCTGATGTATAACACCAAATCTTGCTATGCAGAAATTTGCACAGTGCGGAATTTTCTGCATAGCCGCTCATTCTCCTta                                                               <  1:704795/97‑1 (MQ=255)
                                          aaTACTGTATCGATCCGGCGGATGTATACCACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtat                                                              <  1:873653/97‑1 (MQ=255)
                                              cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                       <  2:3890819/100‑1 (MQ=255)
                                              cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATcc                                                        <  1:1768062/99‑1 (MQ=255)
                                                gTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCa                                                     <  2:573180/100‑1 (MQ=255)
                                                  aTCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATc                                                   >  2:2751594/1‑100 (MQ=255)
                                                   tCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATcc                                                  >  2:4247652/1‑100 (MQ=255)
                                                    cGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCa                                                 >  1:3693862/1‑100 (MQ=255)
                                                    cGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTATGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCa                                                 >  2:3088388/1‑100 (MQ=255)
                                                    agATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCa                                                 >  1:3228814/2‑100 (MQ=255)
                                                     gATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCaa                                                >  1:3268925/1‑100 (MQ=255)
                                                      aTCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTACACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCGTATAAATCCCCATCCAAt                                               <  1:389063/100‑1 (MQ=255)
                                                       tCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAtt                                              <  1:1409818/100‑1 (MQ=255)
                                                           gcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATc                                          >  2:2158028/1‑100 (MQ=255)
                                                             ggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATccc                                        <  1:4222601/100‑1 (MQ=255)
                                                              gATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCt                                       >  2:722705/1‑100 (MQ=255)
                                                                tGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATACCCATCCAATTTATCCCTTc                                     >  1:3398518/1‑100 (MQ=255)
                                                                 gTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAACTTATCCCTTCa                                    <  2:2080754/100‑1 (MQ=255)
                                                                 gTATAACACCAAAACTTGCTACGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                                    <  1:73283/100‑1 (MQ=255)
                                                                  tataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCat                                   >  2:415537/1‑100 (MQ=255)
                                                                  tataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCat                                   >  2:3402505/1‑100 (MQ=255)
                                                                    taACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCatat                                 <  2:1400608/100‑1 (MQ=255)
                                                                     aaCACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTCTCCCTTCatat                                 >  2:2569075/1‑99 (MQ=255)
                                                                     aaCACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATAtt                                >  2:2160339/1‑100 (MQ=255)
                                                                      acacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTc                               >  2:2602022/1‑100 (MQ=255)
                                                                       cacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCa                              >  2:1842933/1‑100 (MQ=255)
                                                                       cacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCa                              >  2:1691419/1‑100 (MQ=255)
                                                                        acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCATTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                             >  1:185106/1‑100 (MQ=255)
                                                                         ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                            >  2:2262563/1‑100 (MQ=255)
                                                                         ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                            >  2:1335867/1‑100 (MQ=255)
                                                                         ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                            >  2:2863718/1‑100 (MQ=255)
                                                                         ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                            >  1:3818420/1‑100 (MQ=255)
                                                                          cAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAtt                           >  1:341303/1‑100 (MQ=255)
                                                                           aaaaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAtt                           <  1:2444328/99‑1 (MQ=255)
                                                                               cTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGtt                       <  1:3488460/99‑1 (MQ=255)
                                                                               cTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTa                      <  1:1386826/100‑1 (MQ=255)
                                                                                 tGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTaaa                    <  1:4268083/100‑1 (MQ=255)
                                                                                 ggCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTaaa                    >  2:814084/2‑100 (MQ=25)
                                                                                 ggCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATAATCCATATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTaaa                    >  2:3023300/2‑100 (MQ=21)
                                                                                  gCTATGCAGAAATTTGCACAGTGCGCACTTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATCTTCAATTAGt                        >  1:2487840/1‑95 (MQ=255)
                                                                                  gCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAAt                   >  2:2643602/1‑100 (MQ=25)
                                                                                  gCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAAt                   >  2:1538638/1‑100 (MQ=25)
                                                                                    tATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCGATCCACTTTATCCCCTCATATTCAATTAGTTAAATaa                 >  1:2971193/1‑100 (MQ=21)
                                                                                      tGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACt               >  2:3698283/1‑100 (MQ=255)
                                                                                      tGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACt               >  2:764239/1‑100 (MQ=255)
                                                                                       gCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCTCCATCCAATTTATCCCTTCATATTCAATTAGTTCAATAACTa              >  1:507336/1‑100 (MQ=25)
                                                                                       gCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa              >  2:3612049/1‑100 (MQ=37)
                                                                                       gCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTAATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa              >  2:918746/1‑100 (MQ=25)
                                                                                       gCAGAAAATTGCGCACTGTGAAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCACCATCCAATTTACCCATTCATATTCAATTAGTCAAATAACTa              >  1:2073822/1‑100 (MQ=21)
                                                                                        cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa             >  2:3871983/1‑100 (MQ=255)
                                                                                        cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa             >  1:2151082/1‑100 (MQ=255)
                                                                                        cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa             >  2:1649091/1‑100 (MQ=255)
                                                                                        cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa             >  2:1389705/1‑100 (MQ=255)
                                                                                        cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa             >  2:290628/1‑100 (MQ=37)
                                                                                          gAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt           >  2:3252124/1‑100 (MQ=255)
                                                                                            aaTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATcc         <  1:1064427/100‑1 (MQ=255)
                                                                                              tttGCACAGTGCGCAATTTTCTGCATAGCCTCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt           <  1:608145/96‑1 (MQ=255)
                                                                                              tttGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaa       <  1:1513439/100‑1 (MQ=255)
                                                                                              tttGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTAACCCTTCATATTCAATTAGCTAAATAACTAAATCCaa       >  2:3856554/1‑100 (MQ=255)
                                                                                               ttGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaat      >  1:1961062/1‑100 (MQ=255)
                                                                                                  cacaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaataat   <  1:104805/100‑1 (MQ=255)
                                                                                                   acaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAAtc  <  1:3681907/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATC  >  NC_000913/2323181‑2323379

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