Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 951,959 A→T 13.6% I263F (ATC→TTC)  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,9590AT13.6% 419.8 / 49.3 154I263F (ATC→TTC) cyoBcytochrome bo terminal oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base A (67/66);  new base T (11/10);  total (78/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.29e-01

Only 100 of 165 total aligned reads displayed.

GAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGG  >  NC_002947/951811‑952107
                                                                                                                                                    |                                                                                                                                                    
gaAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACAt                                                                                                                                                     <  1:2963584/150‑1 (MQ=255)
gaAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACAt                                                                                                                                                     >  2:3949903/1‑150 (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)
                ggCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCtt                                                                                                                                     >  1:1826440/1‑150 (MQ=255)
                   aTGAATCTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTtggt                                                                                                                                  <  2:175706/150‑1 (MQ=255)
                     gAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTtggtgg                                                                                                                                <  1:99577/150‑1 (MQ=255)
                     gAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTtggtgg                                                                                                                                <  1:405070/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)
                          tGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAAcc                                                                                                                           >  1:5310111/1‑150 (MQ=255)
                            aTGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCg                                                                                                                         >  1:3467530/1‑150 (MQ=255)
                              ggACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCgat                                                                                                                       <  1:4862694/150‑1 (MQ=255)
                                aCATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCgatga                                                                                                                     >  1:1752867/1‑150 (MQ=255)
                                 cATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCgatgat                                                                                                                    <  2:2905356/150‑1 (MQ=255)
                                 cATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCgatgat                                                                                                                    <  2:2905250/150‑1 (MQ=255)
                                   tGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGt                                                                                                                  <  2:3467530/150‑1 (MQ=255)
                                     ccGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCACCGAGCTTGGTGGGAACCCGATGATGTAc                                                                                                                <  1:985480/150‑1 (MQ=255)
                                     ccGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTAc                                                                                                                <  1:1160715/150‑1 (MQ=255)
                                     ccGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTAc                                                                                                                >  2:4883245/1‑150 (MQ=255)
                                     ccGATCTGCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAAGCCGATGATGTAc                                                                                                                >  2:1729771/1‑150 (MQ=255)
                                      cGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACg                                                                                                               >  1:500153/1‑150 (MQ=255)
                                      cGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACg                                                                                                               >  1:2474289/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)
                                               ccTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTtggt                                                                                                                                  >  1:4329133/1‑122 (MQ=255)
                                                cTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGtt                                                                                                     <  1:991302/150‑1 (MQ=255)
                                                  ggACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCt                                                                                                   <  1:3377061/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)
                                                          aCCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTgg                                                                                           >  1:2732816/1‑150 (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)
                                                                          tGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGttt                                                                           <  2:4664244/150‑1 (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:20954/1‑150 (MQ=255)
                                                                                   cTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcctga                                                                  <  1:3242196/150‑1 (MQ=255)
                                                                                       gTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTAAGTCAACCTGTTCTGGGCTTGGGGTCACCCTGAGGTTTACatcctgatact                                                              >  2:2527954/1‑150 (MQ=255)
                                                                                       gTTCCCGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATattctggttct                                                              >  2:2717358/1‑145 (MQ=255)
                                                                                            cGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGcctt                                                         >  2:1709085/1‑148 (MQ=255)
                                                                                                  tGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCg                                                   >  1:1668633/1‑150 (MQ=255)
                                                                                                      tgctgcACTCGTTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCgt                                               >  1:4539813/1‑150 (MQ=255)
                                                                                                      tgctgcACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCgt                                               <  2:1965923/150‑1 (MQ=255)
                                                                                                        tgacgcTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATCTTCTGGTTCTGCCTTGCTTCGgtgttt                                             >  2:1241703/6‑148 (MQ=255)
                                                                                                          gcACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTtc                                           <  1:2975723/150‑1 (MQ=255)
                                                                                                            gcTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTtctc                                         <  1:2586198/149‑1 (MQ=255)
                                                                                                              tCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCgg                                       >  1:2117276/1‑150 (MQ=255)
                                                                                                              tCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCgg                                       >  2:4547175/1‑150 (MQ=255)
                                                                                                                gCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGaa                                     <  2:5029/150‑1 (MQ=255)
                                                                                                                gCTCTGCTGACTGTTGAGCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                     <  2:2814423/150‑1 (MQ=255)
                                                                                                                gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCACCGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                     <  2:353416/150‑1 (MQ=255)
                                                                                                                gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                     <  1:1709056/150‑1 (MQ=255)
                                                                                                                gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                     <  2:3955881/150‑1 (MQ=255)
                                                                                                                 ctctGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGaaa                                    <  1:5378846/150‑2 (MQ=255)
                                                                                                                       tgACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                     <  1:3198844/143‑1 (MQ=255)
                                                                                                                       tgACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                     >  2:3198844/1‑143 (MQ=255)
                                                                                                                       tgACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACtt                              >  2:4323487/1‑150 (MQ=255)
                                                                                                                         aCTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCg                            >  2:2485391/1‑150 (MQ=255)
                                                                                                                          cTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGa                           <  1:5313335/150‑1 (MQ=255)
                                                                                                                           tGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGAc                          >  2:5098894/1‑150 (MQ=255)
                                                                                                                           tGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGAc                          >  2:3264248/1‑150 (MQ=255)
                                                                                                                           tGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGAc                          <  2:736376/150‑1 (MQ=255)
                                                                                                                           tGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGAc                          >  1:4041971/1‑150 (MQ=255)
                                                                                                                           tGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGAc                          >  1:380949/1‑150 (MQ=255)
                                                                                                                             ttGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGt                        >  1:1145185/1‑150 (MQ=255)
                                                                                                                               gACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTtc                      >  2:353473/1‑150 (MQ=255)
                                                                                                                                  cGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTctct                   <  1:5446111/150‑1 (MQ=255)
                                                                                                                                   gTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTg                  <  1:4166372/150‑1 (MQ=255)
                                                                                                                                   gTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTg                  >  2:3621622/1‑150 (MQ=255)
                                                                                                                                     tATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGc                <  1:1009657/150‑1 (MQ=255)
                                                                                                                                       tCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCaa              >  2:3864144/1‑150 (MQ=255)
                                                                                                                                       tCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCaa              >  2:2433166/1‑150 (MQ=255)
                                                                                                                                       tCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCaa              >  2:2522080/1‑150 (MQ=255)
                                                                                                                                       tCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCaa              >  2:1013570/1‑150 (MQ=255)
                                                                                                                                         tGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCaagc            <  2:1109616/150‑3 (MQ=255)
                                                                                                                                           gACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGc          <  1:3621622/150‑1 (MQ=255)
                                                                                                                                             cTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCt        <  1:4991199/150‑1 (MQ=255)
                                                                                                                                               tCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGTCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGt      >  1:5197280/1‑150 (MQ=255)
                                                                                                                                               tCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGt      <  2:1826440/150‑1 (MQ=255)
                                                                                                                                                ccACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGtt     >  2:4194207/1‑150 (MQ=255)
                                                                                                                                                ccACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGtt     >  1:2266067/1‑150 (MQ=255)
                                                                                                                                                 cactTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGttt    <  1:4674435/150‑2 (MQ=255)
                                                                                                                                                 cacaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTc    <  2:3087340/150‑1 (MQ=255)
                                                                                                                                                  acaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGc                                                 <  1:1713727/104‑1 (MQ=255)
                                                                                                                                                  acaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCg   >  2:634892/1‑150 (MQ=255)
                                                                                                                                                   caTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCgg  <  2:2854187/150‑1 (MQ=255)
                                                                                                                                                    aTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                     <  1:951448/114‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAAGATGCGCGCACCTGGCATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGG  >  NC_002947/951811‑952107

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: