Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,720,543 C→A 59.8% R197R (CGC→CGA pka → protein lysine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,720,5430CA59.8% 104.9 / 222.3 180R197R (CGC→CGApkaprotein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (36/36);  new base A (59/48);  total (96/84)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

Only 100 of 185 total aligned reads displayed.

GACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGAT  >  NC_000913/2720445‑2720641
                                                                                                  |                                                                                                  
gACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCAcgcg                                                                                                   <  2:1659144/100‑1 (MQ=255)
gACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCAcgcg                                                                                                   >  1:382861/1‑100 (MQ=255)
  cTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGAc                                                                                                 >  2:1437546/1‑100 (MQ=255)
    gggCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAg                                                                                               <  1:2600521/100‑1 (MQ=255)
     ggCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACgaga                                                                                                  <  2:324467/96‑1 (MQ=255)
     cgcgcAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGt                                                                                              >  2:1763458/2‑100 (MQ=255)
       cgcAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTaa                                                                                            >  1:3456809/1‑100 (MQ=255)
         cAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTaaaa                                                                                          >  2:1850852/1‑100 (MQ=255)
         cAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTaaaa                                                                                          >  2:4038589/1‑100 (MQ=255)
         cAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTaaaa                                                                                          >  1:807967/1‑100 (MQ=255)
           aCAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGCCAGTAAAAcc                                                                                        >  2:2457475/1‑100 (MQ=255)
            cAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCa                                                                                       >  2:1305058/1‑100 (MQ=255)
              gCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAgc                                                                                     >  2:896481/1‑100 (MQ=255)
              gCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAgc                                                                                     >  1:2666578/1‑100 (MQ=255)
                gTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAgcgc                                                                                   >  2:1294906/1‑100 (MQ=255)
                  aaGATGTGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCa                                                                                 <  2:3198381/100‑1 (MQ=255)
                     aTGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATcc                                                                              <  2:728550/100‑1 (MQ=255)
                     aTGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATcc                                                                              <  1:3283970/100‑1 (MQ=255)
                         gCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGct                                                                          >  1:3652788/1‑100 (MQ=255)
                         gCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGct                                                                          >  2:1127005/1‑100 (MQ=255)
                         gCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGct                                                                          >  2:1172723/1‑100 (MQ=255)
                         gCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGct                                                                          >  1:3575078/1‑100 (MQ=255)
                           ttttCCTACTTTATTCCGCTCGAAGACAGCCTGGATATCGAAAATGATTAATTGCTTAACTATCTGTCACTAAACATTAAAACCAACGGCATCCta                                                                            >  2:4122564/1‑95 (MQ=255)
                               ccTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTAtct                                                                    <  2:2335720/100‑1 (MQ=255)
                                cTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGGATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCtctctctc                                                                   >  2:1672229/1‑100 (MQ=255)
                                cTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTAtctc                                                                   >  1:3125620/1‑100 (MQ=255)
                                cTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTAtctc                                                                   >  2:55049/1‑100 (MQ=255)
                                 tACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCg                                                                  >  1:2977058/1‑100 (MQ=255)
                                 tACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACAAGCGCCATCCTGCTATATGTCg                                                                  >  2:3519304/1‑100 (MQ=255)
                                 tACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCg                                                                  >  1:2504332/1‑100 (MQ=255)
                                    tttATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAAc                                                               <  1:3804181/100‑1 (MQ=255)
                                     ttATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACa                                                              <  2:923883/100‑1 (MQ=255)
                                     ttATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACa                                                              >  2:1899241/1‑100 (MQ=255)
                                     ttATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACa                                                              >  2:515690/1‑100 (MQ=255)
                                      tATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAg                                                             <  1:1371126/100‑1 (MQ=255)
                                      tATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAg                                                             <  1:3697726/100‑1 (MQ=255)
                                      tATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAg                                                             <  1:2286865/100‑1 (MQ=255)
                                        ttGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGtt                                                           >  1:959301/1‑100 (MQ=255)
                                         tGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTa                                                          >  1:917286/1‑100 (MQ=255)
                                         tGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTa                                                          <  1:841214/100‑1 (MQ=255)
                                           cgcTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAg                                                        >  2:406582/1‑100 (MQ=255)
                                            gcTCGGCGACAGCCTTGATATCGACGTGGATGAATTGCTTGACTATCTGGCACGCGTCAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAg                                                        <  2:1494234/99‑1 (MQ=255)
                                             cTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAgcg                                                      >  1:455198/1‑100 (MQ=255)
                                              tCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAgcga                                                     >  2:135793/1‑100 (MQ=255)
                                               cGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAgcgac                                                    >  1:1569777/1‑100 (MQ=255)
                                                 gCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAgcgacgc                                                  >  2:3273323/1‑100 (MQ=255)
                                                    aCAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAgcgacgcg                                                 >  2:1502387/1‑98 (MQ=255)
                                                      aGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAgcgacgcgcgac                                             <  1:1815410/100‑1 (MQ=255)
                                                       gCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAgcgacgcgcgacg                                            >  2:2368546/1‑100 (MQ=255)
                                                       gCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAgcgacgcgcgacg                                            <  1:3208313/100‑1 (MQ=255)
                                                       gCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAgcgacgcgcgacg                                            >  2:3375283/1‑100 (MQ=255)
                                                       gCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAgcgacgcgcgacg                                            >  2:607495/1‑100 (MQ=255)
                                                        ccTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAgcgacgcgcgacgc                                           >  1:1434566/1‑100 (MQ=255)
                                                        ccTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAgcgacgcgcgacgc                                           >  1:533785/1‑100 (MQ=255)
                                                           ggATATCGACGTTGATGAATTGCTTGACTCTCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCttt                                        <  2:2478559/100‑1 (MQ=255)
                                                            gATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGCCAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTg                                       <  1:420168/100‑1 (MQ=255)
                                                            gATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTg                                       <  1:1760134/100‑1 (MQ=255)
                                                            gATATCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTg                                       >  2:3144844/1‑100 (MQ=255)
                                                             atatCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGTCGCGCGTCGCttaac                                      <  1:2825537/100‑4 (MQ=255)
                                                             atatCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGt                                      <  1:406585/100‑1 (MQ=255)
                                                             atatCGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGt                                      <  1:2460606/100‑1 (MQ=255)
                                                              tatCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGtt                                     >  2:817167/1‑100 (MQ=255)
                                                              tatCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGtt                                     <  1:2423406/100‑1 (MQ=255)
                                                                 cGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTcg                                  <  2:1474826/100‑1 (MQ=255)
                                                                 cGACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTcg                                  >  1:2652741/1‑100 (MQ=255)
                                                                  gACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTcgg                                 >  1:1349361/1‑100 (MQ=255)
                                                                  gACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTcg                                  <  2:3775145/99‑1 (MQ=255)
                                                                  gACGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGTCGCGCGTCGCttaactgttc                                 <  2:2686572/100‑9 (MQ=38)
                                                                   aCGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTcggc                                >  1:150405/1‑100 (MQ=255)
                                                                    cGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTCAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTcggcg                               >  2:2207045/1‑100 (MQ=255)
                                                                    cGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTcggcg                               <  2:3639696/100‑1 (MQ=255)
                                                                    cGTTGATGAATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTcggcg                               >  2:2253896/1‑100 (MQ=255)
                                                                       tgatgaATTGCTTGGCTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCGCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGcc                            <  1:3190792/100‑1 (MQ=255)
                                                                       tgatgaATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGcc                            >  1:3407305/1‑100 (MQ=255)
                                                                        gatgaATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGccc                           >  2:2993961/1‑100 (MQ=255)
                                                                          tgaATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGt                         <  2:1085762/100‑1 (MQ=255)
                                                                           gaATTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTa                        >  1:3993706/1‑100 (MQ=255)
                                                                            aaTTGCTTGACTATCTGGCATGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGCACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAg                       <  2:3474856/100‑1 (MQ=255)
                                                                            aaTTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAg                       <  1:1666226/100‑1 (MQ=255)
                                                                            aaTTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAg                       >  1:2630758/1‑100 (MQ=255)
                                                                             aTTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGt                      >  1:3356320/1‑100 (MQ=255)
                                                                             aTTGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCACCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGt                      <  1:638159/100‑1 (MQ=255)
                                                                              ttGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTg                     >  1:2700098/1‑100 (MQ=255)
                                                                               tGCTTGACTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGc                    <  1:2012057/100‑1 (MQ=255)
                                                                                 cTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCt                  >  1:1038496/1‑100 (MQ=255)
                                                                                  ttGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTc                 <  2:451283/100‑1 (MQ=255)
                                                                                  ttGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTc                 <  1:2368681/100‑1 (MQ=255)
                                                                                    gACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTAACTCGAACAGTTACGCGACGCGAGACGCTTTGTTTCGGACGCACGTAGTGCCTcgc               >  1:841808/1‑100 (MQ=255)
                                                                                      cTATCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGt             <  1:773877/100‑1 (MQ=255)
                                                                                      cTATCTGGCACGAGAAAGTAAAACCGGCGCCATCCTGCGCGGTCGGGAACAGTGGAGCGACGCGTGGCGCTTTGTTTCGGCGGCCCGTAGTGCCTcgcgg             <  2:3738936/100‑2 (MQ=255)
                                                                                         tCTGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGtaat          <  2:1970500/100‑1 (MQ=255)
                                                                                          cTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGtaata         <  1:1198749/100‑1 (MQ=255)
                                                                                           tGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAAAAGTTAGGCGACGCGCGACGCTTTGATTCGCCGGGCCGTAgcgc                    >  2:3837958/1‑88 (MQ=255)
                                                                                           tGGCACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGtaataa        >  1:1910935/1‑100 (MQ=255)
                                                                                            ggCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATaaa       >  2:1921135/1‑100 (MQ=255)
                                                                                              cACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAAcc     <  1:2871986/100‑1 (MQ=255)
                                                                                              cACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAAcc     >  2:263765/1‑100 (MQ=255)
                                                                                              cACGAGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAAcc     <  2:2177339/100‑1 (MQ=255)
                                                                                                cgcgACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGa   >  2:1611312/1‑100 (MQ=255)
                                                                                                  cgACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCggt  >  1:1486193/1‑97 (MQ=255)
                                                                                                  |                                                                                                  
GACTGGGCGCAACAGCGTAAGATGGGCTTTTCCTACTTTATTGCGCTCGGCGACAGCCTGGATATCGACGTTGATGAATTGCTTGACTATCTGGCACGCGACAGTAAAACCAGCGCCATCCTGCTCTATCTCGAACAGTTAAGCGACGCGCGACGCTTTGTTTCGGCGGCCCGTAGTGCCTCGCGTAATAAACCGAT  >  NC_000913/2720445‑2720641

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: