Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 951,917 G→A 5.2% A249T (GCA→ACA) ‡ cyoB → cytochrome bo terminal oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947951,9170GA 5.2% 421.1 / 8.4 133A249T (GCA→ACA) ‡cyoBcytochrome bo terminal oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base G (56/70);  new base A (3/4);  total (59/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 137 total aligned reads displayed.

GGGTACGACGCTTACTGGCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTC  >  NC_002947/951769‑952066
                                                                                                                                                    |                                                                                                                                                     
gggTACGACGCTTACTGGCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACtgctgc                                                                                                                                                      <  1:4545557/150‑1 (MQ=255)
   tACGACGCTTACTGGCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACt                                                                                                                                                   <  2:350541/150‑1 (MQ=255)
    acgacgCTTACTGGCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTc                                                                                                                                                  <  1:2928327/150‑1 (MQ=255)
         gCTTACTGGCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGctct                                                                                                                                             <  2:4330849/150‑1 (MQ=255)
         gCTTACTGGCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGctct                                                                                                                                             >  1:5265913/1‑150 (MQ=255)
           ttACTGGCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTctgc                                                                                                                                           <  1:3353165/150‑1 (MQ=255)
               tGGCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGAc                                                                                                                                       <  1:1767604/150‑1 (MQ=255)
               tGGCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGAc                                                                                                                                       >  1:2854200/1‑150 (MQ=255)
               tGGCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGAc                                                                                                                                       >  1:2854187/1‑150 (MQ=255)
               tGGCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGAc                                                                                                                                       >  1:2369501/1‑150 (MQ=255)
                 gCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTg                                                                                                                                     <  1:3219013/150‑1 (MQ=255)
                   gTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGtt                                                                                                                                   <  1:103438/150‑1 (MQ=255)
                   gTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGtt                                                                                                                                   <  2:1334749/150‑1 (MQ=255)
                   gTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGtt                                                                                                                                   <  1:3897258/150‑1 (MQ=255)
                     gAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGa                                                                                                                                 <  2:595068/150‑1 (MQ=255)
                       aCTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGAcc                                                                                                                               <  1:4151607/150‑1 (MQ=255)
                        cTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCg                                                                                                                              <  2:1543363/150‑1 (MQ=255)
                         ttCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGt                                                                                                                             >  2:4166372/1‑150 (MQ=255)
                          tCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGtt                                                                                                                            <  1:1582966/150‑1 (MQ=255)
                          tCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGtt                                                                                                                            <  2:58220/150‑1 (MQ=255)
                            cTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCAGGGCCAACGTGTTGATCGTCGCTTCGTTACCGATCCTGACTGCAGCACTCGCTCTGCTGACTGTTGACCGTTAt                                                                                                                          >  1:2814423/1‑150 (MQ=255)
                              cGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCt                                                                                                                        <  2:5446286/150‑1 (MQ=255)
                              cGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCt                                                                                                                        <  2:5446170/150‑1 (MQ=255)
                              cGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCt                                                                                                                        <  1:1614816/150‑1 (MQ=255)
                                  aCCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGAc                                                                                                                    >  2:4991199/1‑150 (MQ=255)
                                  aCCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGAc                                                                                                                    >  2:217120/1‑150 (MQ=255)
                                      tgatgaAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTcc                                                                                                                <  2:3737846/150‑1 (MQ=255)
                                          gaAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACAt                                                                                                            >  2:3949903/1‑150 (MQ=255)
                                           aaGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATc                                                                                                           <  2:1238625/150‑1 (MQ=255)
                                            aGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCt                                                                                                          <  1:1723382/150‑1 (MQ=255)
                                              aTGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTc                                                                                                        <  1:1422315/150‑1 (MQ=255)
                                              aTGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTc                                                                                                        <  1:1422364/150‑1 (MQ=255)
                                                  gcgcACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCa                                                                                                    <  1:3712651/150‑1 (MQ=255)
                                                    gcACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAAc                                                                                                  <  2:4727175/150‑1 (MQ=255)
                                                         tGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACAAACGAGCt                                                                                             >  1:4213284/1‑150 (MQ=255)
                                                          ggCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCtt                                                                                            >  1:1826440/1‑150 (MQ=255)
                                                             aTGAATCTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTtggt                                                                                         <  2:175706/150‑1 (MQ=255)
                                                             aTGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCGCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTtggt                                                                                         <  1:1895929/150‑1 (MQ=255)
                                                               gAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTtggtgg                                                                                       <  1:4500605/150‑1 (MQ=255)
                                                               gAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTtggtgg                                                                                       <  1:405070/150‑1 (MQ=255)
                                                               gAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTtggtgg                                                                                       <  1:3949903/150‑1 (MQ=255)
                                                                 aaCTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGa                                                                                     >  1:67557/1‑150 (MQ=255)
                                                                    tGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCAc                                                                                                      >  2:1694618/1‑130 (MQ=255)
                                                                    tGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCAc                                                                                                      <  1:1694618/130‑1 (MQ=255)
                                                                        ggACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCgat                                                                              <  1:4862694/150‑1 (MQ=255)
                                                                           cATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCgatgat                                                                           <  2:2905356/150‑1 (MQ=255)
                                                                           cATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCgatgat                                                                           <  2:2905103/150‑1 (MQ=255)
                                                                             tGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGt                                                                         <  1:3821957/150‑1 (MQ=255)
                                                                             tGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGt                                                                         <  2:3467530/150‑1 (MQ=255)
                                                                              gCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTa                                                                        >  1:5166102/1‑150 (MQ=255)
                                                                              gCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTa                                                                        >  1:984058/1‑150 (MQ=255)
                                                                               ccGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTTCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTAc                                                                       <  2:3104402/150‑1 (MQ=255)
                                                                               ccGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTAc                                                                       >  2:4883245/1‑150 (MQ=255)
                                                                               ccGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTAc                                                                       <  1:985828/150‑1 (MQ=255)
                                                                               ccGATCTGCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAAGCCGATGATGTAc                                                                       >  2:1729771/1‑150 (MQ=255)
                                                                                cGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACg                                                                      >  1:70520/1‑150 (MQ=255)
                                                                                cGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACg                                                                      >  1:500110/1‑150 (MQ=255)
                                                                                    cTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCaa                                                                  <  2:2856287/150‑1 (MQ=255)
                                                                                     ttCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAAc                                                                 <  1:3862139/150‑1 (MQ=255)
                                                                                       cACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCt                                                               >  2:5330419/1‑150 (MQ=255)
                                                                                         ccTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTtggt                                                                                         >  1:4329133/1‑122 (MQ=255)
                                                                                         ccTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTtggt                                                                                         <  2:4329133/122‑1 (MQ=255)
                                                                                          cTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGtt                                                            <  1:991302/150‑1 (MQ=255)
                                                                                            ggACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCt                                                          <  1:3377061/150‑1 (MQ=255)
                                                                                            ggACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCt                                                          <  2:2712530/150‑1 (MQ=255)
                                                                                              aCCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTgg                                                        <  2:3482942/150‑1 (MQ=255)
                                                                                                cTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGc                                                      <  1:807239/150‑1 (MQ=255)
                                                                                                  gCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGt                                                    >  2:2071333/1‑150 (MQ=255)
                                                                                                  gCACCTGGGACAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGATATCTGGACTTCCACATCTTAACCAACGAGCTTGGTGGGAACCCTATGATGTACGTCAACCGGTTCTGTGCGt                                                    >  2:2377431/1‑150 (MQ=255)
                                                                                                    aCCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTgg                                                  >  1:2732816/1‑150 (MQ=255)
                                                                                                    aCCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTgg                                                  <  1:3678091/150‑1 (MQ=255)
                                                                                                    aCCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTgg                                                  >  2:2955297/1‑150 (MQ=255)
                                                                                                     ccTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTggg                                                 <  1:2641634/150‑1 (MQ=255)
                                                                                                          gCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCAc                                            <  1:5330419/150‑1 (MQ=255)
                                                                                                           ccAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCAcc                                           <  2:4798993/150‑1 (MQ=255)
                                                                                                            cAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCAccc                                          >  2:4641779/1‑150 (MQ=255)
                                                                                                               cGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGa                                       <  1:1109981/150‑1 (MQ=255)
                                                                                                                gtgtTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAg                                      >  1:1011561/1‑150 (MQ=255)
                                                                                                                 tgtTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAgg                                     >  1:5157657/1‑150 (MQ=255)
                                                                                                                   ctGATCGTGGCTGCGTTCCCGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGtc                                   <  2:3646590/149‑2 (MQ=255)
                                                                                                                   ctGATCGTGGCTGCGTTCCCGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGtc                                   <  1:627622/149‑2 (MQ=255)
                                                                                                                    tGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGttt                                  <  2:4664259/150‑1 (MQ=255)
                                                                                                                     gATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTa                                 >  1:2424045/1‑150 (MQ=255)
                                                                                                                     gATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCAATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTa                                 >  2:2735733/1‑150 (MQ=255)
                                                                                                                      aTCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTAc                                >  1:1353924/1‑150 (MQ=255)
                                                                                                                      aTCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTAc                                >  2:1843969/1‑150 (MQ=255)
                                                                                                                         gtcgCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACacc                             <  1:37450/150‑3 (MQ=255)
                                                                                                                          tggCTGCGTTCCCGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTatattc                            >  1:1999457/3‑148 (MQ=255)
                                                                                                                          gggCTGCGTTCCCGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTatattc                            <  2:1214651/148‑3 (MQ=255)
                                                                                                                             cTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcctga                         <  1:3242196/150‑1 (MQ=255)
                                                                                                                             cTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcctga                         >  1:20954/1‑150 (MQ=255)
                                                                                                                              ttCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcctgat                        >  2:1059496/1‑150 (MQ=255)
                                                                                                                                cGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGa                                                                                                                     >  2:2103323/1‑55 (MQ=255)
                                                                                                                                cGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGa                                                                                                                     <  1:2103323/55‑1 (MQ=255)
                                                                                                                                 gTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcctgatcct                     >  2:4118558/1‑150 (MQ=255)
                                                                                                                                 gTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTAAGTCAACCTGTTCTGGGCTTGGGGTCACCCTGAGGTTTACatcctgatact                     >  2:2527954/1‑150 (MQ=255)
                                                                                                                                            tGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCg          >  1:1668633/1‑150 (MQ=255)
                                                                                                                                                tgctgcACTCGTTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCgt      >  1:4539813/1‑150 (MQ=255)
                                                                                                                                                tgctgcACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCgt      <  2:1965923/150‑1 (MQ=255)
                                                                                                                                                    gcACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTtc  <  1:2975723/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGGTACGACGCTTACTGGCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTC  >  NC_002947/951769‑952066

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