Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,323,286 C→T 17.8% 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,2860CT17.8% 436.6 / 68.9 184intergenic (+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 (79/72);  new base T (16/17);  total (95/89)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 189 total aligned reads displayed.

GCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTC  >  NC_000913/2323187‑2323385
                                                                                                   |                                                                                                   
 cATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGTGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgcg                                                                                                    >  2:2908475/1‑100 (MQ=255)
  aTCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgcgc                                                                                                   <  2:2433052/100‑1 (MQ=255)
   tCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCa                                                                                                  <  2:1361696/100‑1 (MQ=255)
    cAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCaa                                                                                                 >  2:1412838/1‑100 (MQ=255)
          ccGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGAGGTATAACACCAAAACTTGCTATGCAGAAATTTGCACGGTGCGCAAtt                                                                                               >  1:4408714/1‑96 (MQ=255)
           cGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTg                                                                                          >  2:41234/1‑100 (MQ=255)
           cGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTg                                                                                          <  1:3374888/100‑1 (MQ=37)
            gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGc                                                                                         <  1:225229/100‑1 (MQ=37)
            gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGc                                                                                         <  1:1665562/100‑1 (MQ=37)
            gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGc                                                                                         <  1:1695458/100‑1 (MQ=37)
             tGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCa                                                                                        <  2:4021385/100‑1 (MQ=255)
             tGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCa                                                                                        >  2:1095153/1‑100 (MQ=255)
             tGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCa                                                                                        <  2:258344/100‑1 (MQ=25)
              gCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCAt                                                                                       <  1:2844824/100‑1 (MQ=25)
               cATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATa                                                                                      <  1:493473/100‑1 (MQ=255)
                aTTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAg                                                                                     <  2:3342477/100‑1 (MQ=25)
                aTTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAg                                                                                     <  2:4227855/100‑1 (MQ=25)
                 ttGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGc                                                                                    >  1:4187230/1‑100 (MQ=255)
                 ttGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGc                                                                                    <  2:3365243/100‑1 (MQ=25)
                 ttGATGTATAAACTCCAGGAATACGGTAGCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGc                                                                                    <  1:1380853/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)
                          aaaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtc                                                                           <  1:1953477/100‑1 (MQ=255)
                           aaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtct                                                                          <  2:971847/100‑1 (MQ=255)
                           aaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtct                                                                          >  1:1160077/1‑100 (MQ=255)
                             cTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                        >  2:1076493/1‑100 (MQ=255)
                             cTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                        >  2:4236815/1‑100 (MQ=255)
                              tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                                       <  2:2944081/100‑1 (MQ=255)
                                cAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTTCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                                       <  1:2714760/98‑1 (MQ=255)
                                cAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTta                                                                     <  1:1681807/100‑1 (MQ=255)
                                cAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTta                                                                     <  1:2884193/100‑1 (MQ=255)
                                 aGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtat                                                                    <  2:2875688/100‑1 (MQ=255)
                                   gAATACGGTATCGATCCAGCTGATGTATAACACCAAATCTTGCTATGCAGAAATTTGCACAGTGCGGAATTTTCTGCATAGCCGCTCATTCTCCTta                                                                     <  1:704795/97‑1 (MQ=255)
                                        cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                             <  2:2603907/100‑1 (MQ=255)
                                        cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                             <  2:3890819/100‑1 (MQ=255)
                                        cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATcc                                                              <  1:1768062/99‑1 (MQ=255)
                                            aTCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATc                                                         <  2:3145723/100‑1 (MQ=255)
                                            aTCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATc                                                         >  2:2751594/1‑100 (MQ=255)
                                            aTCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATc                                                         >  2:4288827/1‑100 (MQ=255)
                                              cGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCa                                                       >  2:2710708/1‑100 (MQ=255)
                                              cGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCa                                                       >  1:3693862/1‑100 (MQ=255)
                                              agATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCa                                                       >  1:3228814/2‑100 (MQ=255)
                                               gATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCaa                                                      >  1:843340/1‑100 (MQ=255)
                                               gATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCaa                                                      >  2:1521163/1‑100 (MQ=255)
                                                aTCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTACACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCGTATAAATCCCCATCCAAt                                                     <  1:389063/100‑1 (MQ=255)
                                                 tCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAtt                                                    <  1:1409818/100‑1 (MQ=255)
                                                      cggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATcc                                               >  2:72033/1‑100 (MQ=255)
                                                       ggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATccc                                              <  1:4222601/100‑1 (MQ=255)
                                                          tGTATAACACCAAAACTTGCTACGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTc                                           >  2:73283/1‑100 (MQ=255)
                                                            tataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCat                                         >  2:415537/1‑100 (MQ=255)
                                                              taACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCatat                                       <  2:1400608/100‑1 (MQ=255)
                                                                 cacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCa                                    >  2:1691419/1‑100 (MQ=255)
                                                                 cacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCa                                    >  2:1842933/1‑100 (MQ=255)
                                                                 cacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCa                                    >  2:546120/1‑100 (MQ=255)
                                                                  acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                                   >  2:3568763/1‑100 (MQ=255)
                                                                  acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                                   >  2:2158834/1‑100 (MQ=255)
                                                                  acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                                   >  2:2230224/1‑100 (MQ=255)
                                                                  acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                                   >  2:596128/1‑100 (MQ=255)
                                                                  acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCATTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                                   >  1:185106/1‑100 (MQ=255)
                                                                   ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                                  >  2:2863718/1‑100 (MQ=255)
                                                                   ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                                  >  1:3818420/1‑100 (MQ=255)
                                                                   ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                                  >  2:2100571/1‑100 (MQ=255)
                                                                   ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                                  >  2:2262563/1‑100 (MQ=255)
                                                                   ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                                  >  2:1335867/1‑100 (MQ=255)
                                                                    cAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAtt                                 <  2:2028735/100‑1 (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:3401840/1‑100 (MQ=255)
                                                                          cggCTATGCAGAAATTTGCACAGTGCTCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTaa                           <  1:1218056/98‑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)
                                                                            gCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAGTTTATCCCTTCATATTCAATTAGTTAAAt                         <  1:2420610/100‑1 (MQ=255)
                                                                            gCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAAt                         >  2:3663887/1‑100 (MQ=255)
                                                                              tATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCGATCCACTTTATCCCCTCATATTCAATTAGTTAAATaa                       >  1:2971193/1‑100 (MQ=21)
                                                                                tGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACt                     >  2:2986864/1‑100 (MQ=25)
                                                                                 gCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCTCCATCCAATTTATCCCTTCATATTCAATTAGTTCAATAACTa                    >  1:507336/1‑100 (MQ=25)
                                                                                 gCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTAATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa                    >  2:918746/1‑100 (MQ=25)
                                                                                 gCAGAAAATTGCGCACTGTGAAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCACCATCCAATTTACCCATTCATATTCAATTAGTCAAATAACTa                    >  1:2073822/1‑100 (MQ=21)
                                                                                  cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa                   >  2:1649091/1‑100 (MQ=255)
                                                                                  cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa                   >  1:2151082/1‑100 (MQ=255)
                                                                                  cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa                   >  2:3871983/1‑100 (MQ=255)
                                                                                  cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa                   >  2:1389705/1‑100 (MQ=255)
                                                                                  cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa                   >  1:779679/1‑100 (MQ=37)
                                                                                  cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa                   >  2:1491938/1‑100 (MQ=37)
                                                                                    gAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaaa                  <  2:3724493/99‑1 (MQ=255)
                                                                                      aaTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATcc               <  1:1064427/100‑1 (MQ=255)
                                                                                        tttGCACAGTGCGCAATTTTCTGCATAGCCTCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt                 <  1:608145/96‑1 (MQ=255)
                                                                                        tttGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaa             <  1:1513439/100‑1 (MQ=255)
                                                                                         ttGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaat            >  1:1961062/1‑100 (MQ=255)
                                                                                           gCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaataa          >  2:1548423/1‑100 (MQ=255)
                                                                                            cacaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaataat         <  1:104805/100‑1 (MQ=255)
                                                                                             acaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAAtc        <  1:3681907/100‑1 (MQ=255)
                                                                                              caGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAAtct       <  2:3480729/100‑1 (MQ=255)
                                                                                                   cgcAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTc  >  1:4324264/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTC  >  NC_000913/2323187‑2323385

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