Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,323,275 T→A 18.1% 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,2750TA18.1% 469.5 / 79.7 182intergenic (+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 (83/65);  new base A (16/17);  total (99/83)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

TTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAA  >  NC_000913/2323176‑2323374
                                                                                                   |                                                                                                   
tttAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAt                                                                                                     >  1:1321042/1‑100 (MQ=255)
   aaTGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTg                                                                                                  >  1:4342238/1‑100 (MQ=255)
    aTGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGc                                                                                                 <  2:2133219/100‑1 (MQ=255)
       cTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGcaca                                                                                              >  2:4315235/1‑100 (MQ=255)
        tGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAg                                                                                             >  2:3047602/1‑100 (MQ=255)
           gCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgc                                                                                          >  2:190649/1‑100 (MQ=255)
             aTCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTgcgc                                                                                        <  2:2433052/100‑1 (MQ=255)
              tCAGTCGCCGTGCATTGATGTATAAACTCCAGAAATACGGTATCGCTCCGGCGGATGTATAACACCAAACCTTGCTATGCAGAAAATTGCGCACTGTGCa                                                                                       <  2:3885042/100‑1 (MQ=37)
               cAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCaa                                                                                      >  2:1412838/1‑100 (MQ=255)
               cAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCaa                                                                                      >  1:3757300/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)
                       ttGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGc                                                                              <  1:1095150/99‑1 (MQ=255)
                       gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGc                                                                              <  1:1695458/100‑1 (MQ=37)
                        tGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCa                                                                             <  2:4021385/100‑1 (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:4227855/100‑1 (MQ=25)
                            ttGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGc                                                                         <  2:3362274/100‑1 (MQ=25)
                             tGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGcc                                                                        >  2:181097/1‑100 (MQ=255)
                             tGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAACTTGCACAGTGCGCAATTTTCTGCATAGcc                                                                        <  2:2565345/100‑1 (MQ=255)
                              gATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCg                                                                       >  2:2489065/1‑100 (MQ=255)
                                tGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGc                                                                      <  1:1910935/99‑1 (MQ=255)
                                tGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCt                                                                     <  1:1038279/100‑1 (MQ=255)
                                tGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCt                                                                     <  2:1859006/100‑1 (MQ=255)
                                   ataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATATCACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCa                                                                   >  2:2828715/1‑99 (MQ=255)
                                     aaaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtc                                                                >  1:222557/1‑100 (MQ=255)
                                      aaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtct                                                               <  2:971847/100‑1 (MQ=255)
                                        cTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                             >  2:1076493/1‑100 (MQ=255)
                                         tCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                            <  2:2944017/100‑1 (MQ=255)
                                          ccAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCtt                                                           >  1:3690661/1‑100 (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)
                                               aaTACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGTATAGCCGCTCATTCTCCTTATaaa                                                      <  2:1348384/100‑1 (MQ=255)
                                                   cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                  <  2:3890819/100‑1 (MQ=255)
                                                   cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                  <  2:2603907/100‑1 (MQ=255)
                                                     gTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCa                                                <  2:573180/100‑1 (MQ=255)
                                                       aTCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATc                                              <  2:3145723/100‑1 (MQ=255)
                                                        tCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATcc                                             >  2:4247652/1‑100 (MQ=255)
                                                         cGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCa                                            >  2:2710708/1‑100 (MQ=255)
                                                          gATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCaa                                           >  1:843340/1‑100 (MQ=255)
                                                          gATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCaa                                           >  2:1521163/1‑100 (MQ=255)
                                                          gATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCaa                                           >  1:3268925/1‑100 (MQ=255)
                                                           aTCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTACACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCGTATAAATCCCCATCCAAt                                          <  1:389063/100‑1 (MQ=255)
                                                            tCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAtt                                         <  1:1409818/100‑1 (MQ=255)
                                                             ccGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAAtt                                         <  1:181097/99‑1 (MQ=255)
                                                                   gATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCt                                  >  2:722705/1‑100 (MQ=255)
                                                                     tGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATACCCATCCAATTTATCCCTTc                                >  1:3398518/1‑100 (MQ=255)
                                                                      gTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                               <  1:4236783/100‑1 (MQ=255)
                                                                      gTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAACTTATCCCTTCa                               <  2:2080754/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)
                                                                           acacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTc                          <  1:3022784/100‑1 (MQ=255)
                                                                            cacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCa                         >  2:546120/1‑100 (MQ=255)
                                                                            cacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCa                         >  2:1691419/1‑100 (MQ=255)
                                                                            cacCAAAACTTCCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCa                         <  2:1132206/100‑1 (MQ=255)
                                                                             acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                        >  2:2230224/1‑100 (MQ=255)
                                                                             acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                        >  2:2158834/1‑100 (MQ=255)
                                                                             acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                        >  2:3568763/1‑100 (MQ=255)
                                                                              ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                       >  2:2262563/1‑100 (MQ=255)
                                                                              ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                       >  1:3818420/1‑100 (MQ=255)
                                                                              ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                       >  2:2863718/1‑100 (MQ=255)
                                                                               cAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAtt                      >  1:341303/1‑100 (MQ=255)
                                                                                aaaaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAtt                      <  1:2444328/99‑1 (MQ=255)
                                                                                   aCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGtt                  >  2:1370512/1‑100 (MQ=255)
                                                                                    cTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGtt                  <  1:3488460/99‑1 (MQ=255)
                                                                                    cTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTa                 <  1:1386826/100‑1 (MQ=255)
                                                                                     cggCTATGCAGAAATTTGCACAGTGCTCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTaa                <  1:1218056/98‑1 (MQ=255)
                                                                                      ggCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTaaa               >  2:814084/2‑100 (MQ=25)
                                                                                       gCTATGCAGAAATTTGCACAGTGCGCACTTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATCTTCAATTAGt                   >  1:2487840/1‑95 (MQ=255)
                                                                                       gCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAAt              >  2:3663887/1‑100 (MQ=255)
                                                                                         tATGCAGAAATTTGCACAGTTCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTGCATATTCAATTAGTTAAATa             <  1:4192920/99‑1 (MQ=255)
                                                                                         tATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCGATCCACTTTATCCCCTCATATTCAATTAGTTAAATaa            >  1:2971193/1‑100 (MQ=21)
                                                                                          aTGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAAc           >  1:1545066/1‑100 (MQ=255)
                                                                                           tGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACt          >  2:764239/1‑100 (MQ=255)
                                                                                            gCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCTCCATCCAATTTATCCCTTCATATTCAATTAGTTCAATAACTa         >  1:507336/1‑100 (MQ=25)
                                                                                            gCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTAATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa         >  2:918746/1‑100 (MQ=25)
                                                                                            gCAGAAAATTGCGCACTGTGAAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCACCATCCAATTTACCCATTCATATTCAATTAGTCAAATAACTa         >  1:2073822/1‑100 (MQ=21)
                                                                                             cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa        >  2:1389705/1‑100 (MQ=255)
                                                                                             cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa        >  1:2151082/1‑100 (MQ=255)
                                                                                             cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa        >  2:1649091/1‑100 (MQ=255)
                                                                                             cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa        >  2:1491938/1‑100 (MQ=37)
                                                                                             cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa        >  2:290628/1‑100 (MQ=37)
                                                                                             cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa        >  1:779679/1‑100 (MQ=37)
                                                                                              aGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaaa       <  1:1054913/100‑1 (MQ=255)
                                                                                               gAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaaa       <  2:3724493/99‑1 (MQ=255)
                                                                                               gAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt      >  2:3252124/1‑100 (MQ=255)
                                                                                                aaaTTTGCACAGTTCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt      <  1:1953914/99‑1 (MQ=255)
                                                                                                   tttGCACAGTGCGCAATTTTCTGCATAGCCTCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt      <  1:608145/96‑1 (MQ=255)
                                                                                                   tttGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaa  <  1:1513439/100‑1 (MQ=255)
                                                                                                   tttGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTAACCCTTCATATTCAATTAGCTAAATAACTAAATCCaa  >  2:3856554/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TTTAATGCTGGGCATCAGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAA  >  NC_000913/2323176‑2323374

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: