Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,323,291 T→A 17.4% intergenic (+40/‑156) 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,2910TA17.4% 371.5 / 64.9 149intergenic (+40/‑156)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 (68/55);  new base A (14/12);  total (82/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

Only 100 of 160 total aligned reads displayed.

AGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCA  >  NC_000913/2323192‑2323390
                                                                                                   |                                                                                                   
aGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAt                                                                                                     <  1:1483866/100‑1 (MQ=255)
aGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAt                                                                                                     >  1:688518/1‑100 (MQ=255)
 gTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtt                                                                                                    >  1:3891254/1‑100 (MQ=255)
 gTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAt                                                                                                     >  2:1285168/1‑99 (MQ=255)
   cGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAAtttt                                                                                                  >  1:602243/1‑100 (MQ=255)
    gCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTc                                                                                                 <  2:75946/100‑1 (MQ=255)
      cGTGCATTGATGTATAAACTCCCGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAGTTTGCACAGTGCGCAATTTTCTg                                                                                               >  2:2438339/1‑100 (MQ=255)
       gTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGc                                                                                              <  1:3718254/100‑1 (MQ=25)
        tGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCa                                                                                             >  1:4180619/1‑100 (MQ=255)
         gCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCAt                                                                                            >  2:3887397/1‑100 (MQ=255)
           aTTTATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTAAAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAg                                                                                          <  2:3818330/100‑1 (MQ=21)
           aTTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGCAATTTGCAAAGTGGGCAATTTTCTGCAt                                                                                            >  2:3534604/1‑98 (MQ=255)
           aTTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAg                                                                                          >  1:3788559/1‑100 (MQ=255)
            ttGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGc                                                                                         <  1:3300665/100‑1 (MQ=25)
            ttGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAAAAAATTGCGCACTGTGCAAATTTCTGCATAGc                                                                                         <  1:3977109/100‑1 (MQ=25)
             tGATGTATAAGCACCAGGAATACGGTATCGATCCGGAGGATGTAGAACACCAAAACCTGCGATGCAGAAAATGGCGCACTGTGCAAATTTCTGCATAg                                                                                          <  2:308833/98‑1 (MQ=18)
             tGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGcc                                                                                        >  1:2961259/1‑100 (MQ=255)
             tGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGca                                                                                        <  2:2652930/100‑2 (MQ=25)
                  tataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCa                                                                                   >  1:1518064/1‑100 (MQ=255)
                  tataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCa                                                                                   >  1:1261290/1‑100 (MQ=255)
                   ataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCAt                                                                                  <  2:1422197/100‑1 (MQ=255)
                   ataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCAt                                                                                  <  1:2675281/100‑1 (MQ=255)
                   ataAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTACGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCAt                                                                                  >  2:2675274/1‑100 (MQ=255)
                     aaaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtc                                                                                >  1:1650164/1‑100 (MQ=255)
                     aaaCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtc                                                                                >  1:2406586/1‑100 (MQ=255)
                      naCTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATtct                                                                               >  2:3730996/2‑100 (MQ=255)
                        cTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTcc                                                                             <  2:1483290/100‑1 (MQ=255)
                          ccAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCt                                                                            <  1:516828/99‑1 (MQ=255)
                           cAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTta                                                                          >  1:3571356/1‑100 (MQ=255)
                            aGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtat                                                                         >  1:1236450/1‑100 (MQ=255)
                             ggAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTtata                                                                        <  1:3284673/100‑1 (MQ=255)
                                aTACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAAt                                                                     <  1:668185/100‑1 (MQ=255)
                                aTACGGAATAGATCGGGTGGATGTTTCGCACACAAACTTACTATGCAAAATGTTGCACAGGGCGCTATTTTCTGCATAGCCACtaa                                                                                   <  2:3903108/86‑3 (MQ=255)
                                  aCGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATcc                                                                   <  1:2889381/100‑1 (MQ=255)
                                   cGGTATCGATCGGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTATGCATAGCCGCTCATTCTCCTTATAAAt                                                                     >  2:609472/1‑97 (MQ=255)
                                   cGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATccc                                                                  >  1:991886/1‑100 (MQ=255)
                                    ggTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATcccc                                                                 >  1:815925/1‑100 (MQ=255)
                                     gTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCa                                                                >  1:579316/1‑100 (MQ=255)
                                       aTCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATc                                                              >  2:1062599/1‑100 (MQ=255)
                                         cGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCa                                                            <  1:393225/100‑1 (MQ=255)
                                          gATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCaa                                                           <  2:991886/100‑1 (MQ=255)
                                           aTCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCAGCTCATTCTCCTTATAAATCCCCATCCAAt                                                          >  1:2082482/1‑100 (MQ=255)
                                               ggcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTAt                                                      >  2:2165315/1‑100 (MQ=255)
                                               ggcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTAt                                                      >  2:1490384/1‑100 (MQ=255)
                                               ggcggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTAt                                                      >  1:1353900/1‑100 (MQ=255)
                                                 cggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATcc                                                    <  2:2255820/100‑1 (MQ=255)
                                                 cggATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATcc                                                    <  2:1236449/100‑1 (MQ=255)
                                                   gATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCt                                                  >  2:2523364/1‑100 (MQ=255)
                                                   gATGTATAACACCAAAACTTGCTATGCAGAAATTCGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCCATTTATCCCt                                                  <  1:2928286/100‑1 (MQ=255)
                                                      gTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                                               <  2:2497773/100‑1 (MQ=255)
                                                      gTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCa                                               <  2:1228390/100‑1 (MQ=255)
                                                       tataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCat                                              <  2:3177535/100‑1 (MQ=255)
                                                       tataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCat                                              >  2:1658172/1‑100 (MQ=255)
                                                       tataACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCat                                              >  2:3298030/1‑100 (MQ=255)
                                                           acacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTc                                          <  1:1995793/100‑1 (MQ=255)
                                                            cacCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCa                                         <  2:3597828/100‑1 (MQ=255)
                                                             acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                                        >  1:95942/1‑100 (MQ=255)
                                                             acCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCaa                                        >  2:1803120/1‑100 (MQ=255)
                                                              ccAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                                       >  1:367244/1‑100 (MQ=255)
                                                              ccAAAACTTGCTATGCAGAAATTTGCACAGTGCACAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAt                                       >  2:4001018/1‑100 (MQ=255)
                                                               cAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAtt                                      >  1:849961/1‑100 (MQ=255)
                                                               cAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAtt                                      <  1:2054923/100‑1 (MQ=255)
                                                               cAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAAtt                                      <  1:121517/100‑1 (MQ=255)
                                                                  aaCTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGt                                   <  1:883/100‑1 (MQ=255)
                                                                    cTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTa                                 >  1:232628/1‑100 (MQ=255)
                                                                    cTTGATATGTAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCATTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGt                                   <  1:3158421/98‑1 (MQ=255)
                                                                      tGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTaaa                               <  2:481636/100‑1 (MQ=255)
                                                                       gCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAAt                              <  1:3842369/100‑1 (MQ=255)
                                                                        cTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCTCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAg                                    <  2:719967/93‑1 (MQ=255)
                                                                        cTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATa                             >  1:268106/1‑100 (MQ=25)
                                                                        cTATGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATa                             >  1:914887/1‑100 (MQ=25)
                                                                           tGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACt                          <  1:3379193/100‑1 (MQ=255)
                                                                           tGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACt                          >  1:3472335/1‑100 (MQ=255)
                                                                           tGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACt                          >  2:3068883/1‑100 (MQ=37)
                                                                           tGCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACt                          >  1:4109732/1‑100 (MQ=37)
                                                                            gCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa                         >  2:193834/1‑100 (MQ=255)
                                                                            gCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa                         >  2:33739/1‑100 (MQ=255)
                                                                            gCAGAAAATTGCGCACTGTGCAAATTTCTGCCTAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa                         >  1:1895799/1‑100 (MQ=37)
                                                                            gCAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTa                         >  1:1086528/1‑100 (MQ=25)
                                                                             cAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCACCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa                        >  2:1836556/1‑100 (MQ=255)
                                                                             cAGAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa                        >  1:395521/1‑100 (MQ=37)
                                                                             cAGAAAATTGCGCACTGTGCAAATTTCTACATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTaa                        >  2:860231/1‑100 (MQ=37)
                                                                              aGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCTCAACCCATTCATCCCTTCATATTCAATTAGTTAAATAACTaa                        <  2:2739779/99‑1 (MQ=255)
                                                                               gAAATTTGCACCGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATACAATATATCCCTTCATATTCAATTAGTTAAATAACTAAAt                      <  2:437909/100‑1 (MQ=255)
                                                                               gAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt                      >  1:3148143/1‑100 (MQ=255)
                                                                               gAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt                      >  1:4078157/1‑100 (MQ=255)
                                                                               gAAAATTGCGCACTGTGCAAATTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAAt                      >  1:2231319/1‑100 (MQ=37)
                                                                                 aaTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATcc                    <  2:2613415/100‑1 (MQ=255)
                                                                                 aaTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATcc                    <  1:352031/100‑1 (MQ=255)
                                                                                  aTTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCa                   <  1:3372681/100‑1 (MQ=255)
                                                                                     tGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaata                >  1:4176280/1‑100 (MQ=255)
                                                                                       cacaGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCaataat              >  2:3240943/1‑100 (MQ=255)
                                                                                             gcgcAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCAtt        <  1:835203/100‑1 (MQ=255)
                                                                                             gcgcAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCAt         >  2:1365290/1‑99 (MQ=255)
                                                                                              cgcAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTc       >  2:3806292/1‑100 (MQ=255)
                                                                                                cAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTAAAATAACTAAATCCAATAATCTCATTCTg     >  1:1556310/1‑100 (MQ=255)
                                                                                                  aTTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGc   <  2:4180590/100‑1 (MQ=255)
                                                                                                  aTTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGc   >  1:1128038/1‑100 (MQ=255)
                                                                                                   ttttCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCa  <  2:352192/100‑1 (MQ=255)
                                                                                                   ttttCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCa  >  1:3189462/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AGTCGCCGTGCATTGATGTATAAACTCCAGGAATACGGTATCGATCCGGCGGATGTATAACACCAAAACTTGCTATGCAGAAATTTGCACAGTGCGCAATTTTCTGCATAGCCGCTCATTCTCCTTATAAATCCCCATCCAATTTATCCCTTCATATTCAATTAGTTAAATAACTAAATCCAATAATCTCATTCTGGCA  >  NC_000913/2323192‑2323390

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: