Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 951,919 A→G 5.1% A249A (GCA→GCG) ‡ 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,9190AG 5.1% 439.7 / 7.7 137A249A (GCA→GCG) ‡cyoBcytochrome bo terminal oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base A (58/71);  new base G (3/4);  total (61/76)
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.72e-01

Only 100 of 137 total aligned reads displayed.

TACGACGCTTACTGGCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCT  >  NC_002947/951772‑952067
                                                                                                                                                   |                                                                                                                                                    
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:2854200/1‑150 (MQ=255)
            tGGCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGAc                                                                                                                                        >  1:2854187/1‑150 (MQ=255)
            tGGCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGAc                                                                                                                                        <  1:1767604/150‑1 (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)
                  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                                                                                                                             <  2:58220/150‑1 (MQ=255)
                       tCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGtt                                                                                                                             <  1:1582966/150‑1 (MQ=255)
                         cTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCAGGGCCAACGTGTTGATCGTCGCTTCGTTACCGATCCTGACTGCAGCACTCGCTCTGCTGACTGTTGACCGTTAt                                                                                                                           >  1:2814423/1‑150 (MQ=255)
                          tCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATc                                                                                                                          >  1:3925392/1‑150 (MQ=255)
                           cGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCt                                                                                                                         <  2:5446170/150‑1 (MQ=255)
                           cGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCt                                                                                                                         <  1:1614816/150‑1 (MQ=255)
                           cGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCt                                                                                                                         <  2:5446286/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:1422364/150‑1 (MQ=255)
                                           aTGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTc                                                                                                         <  1:1422315/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:3949903/150‑1 (MQ=255)
                                                            gAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTtggtgg                                                                                        <  1:4500605/150‑1 (MQ=255)
                                                              aaCTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGa                                                                                      >  1:67557/1‑150 (MQ=255)
                                                                 tGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCAc                                                                                                       <  1:1694618/130‑1 (MQ=255)
                                                                 tGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCAc                                                                                                       >  2:1694618/1‑130 (MQ=255)
                                                                   aTGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCg                                                                                 >  1:3467530/1‑150 (MQ=255)
                                                                     ggACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCgat                                                                               <  1:4862694/150‑1 (MQ=255)
                                                                        cATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCgatgat                                                                            <  2:2905103/150‑1 (MQ=255)
                                                                        cATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCgatgat                                                                            <  2:2905356/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:500110/1‑150 (MQ=255)
                                                                             cGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACg                                                                       >  1:70520/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                                                                                          <  2:4329133/122‑1 (MQ=255)
                                                                                       cTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGtt                                                             <  1:991302/150‑1 (MQ=255)
                                                                                         ggACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCt                                                           <  2:2712530/150‑1 (MQ=255)
                                                                                         ggACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCt                                                           <  1:3377061/150‑1 (MQ=255)
                                                                                           aCCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTgg                                                         <  2:3482942/150‑1 (MQ=255)
                                                                                             cTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGc                                                       <  1:807239/150‑1 (MQ=255)
                                                                                              tGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCg                                                      >  1:4424228/1‑150 (MQ=255)
                                                                                               gCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGt                                                     >  2:2071333/1‑150 (MQ=255)
                                                                                               gCACCTGGGACAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGATATCTGGACTTCCACATCTTAACCAACGAGCTTGGTGGGAACCCTATGATGTACGTCAACCGGTTCTGTGCGt                                                     >  2:2377431/1‑150 (MQ=255)
                                                                                                cACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTg                                                    <  1:4076295/150‑1 (MQ=255)
                                                                                                 aCCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTgg                                                   >  2:2955297/1‑150 (MQ=255)
                                                                                                 aCCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTgg                                                   >  1:2732816/1‑150 (MQ=255)
                                                                                                 aCCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTgg                                                   <  1:3678091/150‑1 (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)
                                                                                                             gtgtTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAg                                       >  1:3761555/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)
                                                                                                                       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)
                                                                                                                                                   |                                                                                                                                                    
TACGACGCTTACTGGCGTGAACTTCCTCGTCACCGTGATGAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCT  >  NC_002947/951772‑952067

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: