Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 952,021 T→A 15.3% G283G (GGT→GGA cyoB → cytochrome bo terminal oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947952,0210TA15.3% 430.5 / 61.3 163G283G (GGT→GGAcyoBcytochrome bo terminal oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base T (86/52);  new base A (14/11);  total (100/63)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 164 total aligned reads displayed.

GCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCA  >  NC_002947/951875‑952170
                                                                                                                                                  |                                                                                                                                                     
gCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCAc                                                                                                                                                    <  1:5330419/150‑1 (MQ=255)
 ccAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCAcc                                                                                                                                                   <  2:4798993/150‑1 (MQ=255)
  cAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCAccc                                                                                                                                                  <  2:3925392/150‑1 (MQ=255)
  cAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCAccc                                                                                                                                                  >  2:4641779/1‑150 (MQ=255)
     cGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGa                                                                                                                                               >  2:2101582/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)
         ctGATCGTGGCTGCGTTCCCGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGtc                                                                                                                                           <  1:627622/149‑2 (MQ=255)
         ctGATCGTGGCTGCGTTCCCGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGtc                                                                                                                                           <  2:3646590/149‑2 (MQ=255)
          tGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGttt                                                                                                                                          <  2:4664244/150‑1 (MQ=255)
           gATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTa                                                                                                                                         >  1:2424045/1‑150 (MQ=255)
           gATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCAATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTa                                                                                                                                         >  2:2735733/1‑150 (MQ=255)
            aTCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTAc                                                                                                                                        >  2:1843969/1‑150 (MQ=255)
                gggCTGCGTTCCCGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTatattc                                                                                                                                    <  2:1214651/148‑3 (MQ=255)
                   cTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcctga                                                                                                                                 >  1:20954/1‑150 (MQ=255)
                   cTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcctga                                                                                                                                 <  1:3242196/150‑1 (MQ=255)
                    ttCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcctgat                                                                                                                                >  2:1059496/1‑150 (MQ=255)
                       gTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTAAGTCAACCTGTTCTGGGCTTGGGGTCACCCTGAGGTTTACatcctgatact                                                                                                                             >  2:2527954/1‑150 (MQ=255)
                            cGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGcctt                                                                                                                        >  2:1709085/1‑148 (MQ=255)
                                  tGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCg                                                                                                                  >  1:1668633/1‑150 (MQ=255)
                                    aCTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGc                                                                                                                >  2:3573632/1‑150 (MQ=255)
                                      tgctgcACTCGTTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCgt                                                                                                              >  1:4539813/1‑150 (MQ=255)
                                            gcTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTtctc                                                                                                        <  1:2586198/149‑1 (MQ=255)
                                              tCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCgg                                                                                                      >  2:4547175/1‑150 (MQ=255)
                                              tCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCgg                                                                                                      >  1:2117276/1‑150 (MQ=255)
                                               cGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGa                                                                                                     >  1:1713225/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                                                                                                    <  2:3955881/150‑1 (MQ=255)
                                                gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                                                                    <  2:5265913/150‑1 (MQ=255)
                                                gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                                                                    <  1:1709056/150‑1 (MQ=255)
                                                  tctGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGt                                                                                                  <  1:2065649/150‑1 (MQ=255)
                                                       tgACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                                                                    <  1:3198844/143‑1 (MQ=255)
                                                       tgACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                                                                    >  2:3198844/1‑143 (MQ=255)
                                                         aCTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCg                                                                                           >  2:2485391/1‑150 (MQ=255)
                                                          cTGTTGACCGTTATCTGGACTTCCACATCTTCACCATCGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGa                                                                                          <  1:5313316/150‑1 (MQ=255)
                                                           tGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGAc                                                                                         >  2:5098894/1‑150 (MQ=255)
                                                           tGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGAc                                                                                         <  2:736376/150‑1 (MQ=255)
                                                           tGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGAc                                                                                         >  1:380949/1‑150 (MQ=255)
                                                             ttGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGt                                                                                       >  1:1145185/1‑150 (MQ=255)
                                                                aCCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTtct                                                                                    >  1:2212542/1‑150 (MQ=255)
                                                                  cGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTctct                                                                                  <  1:5446111/150‑1 (MQ=255)
                                                                     tATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGc                                                                               <  2:4564146/150‑1 (MQ=255)
                                                                     tATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGc                                                                               <  1:1009657/150‑1 (MQ=255)
                                                                       tCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCaa                                                                             >  2:2433166/1‑150 (MQ=255)
                                                                       tCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCaa                                                                             >  2:3864144/1‑150 (MQ=255)
                                                                       tCTGGATTTCCACTTCTTCACCAATGATGGTGGAGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCaa                                                                             >  2:2521844/1‑150 (MQ=255)
                                                                       tCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCaa                                                                             >  2:1013570/1‑150 (MQ=255)
                                                                           gACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGc                                                                         <  1:3621622/150‑1 (MQ=255)
                                                                             cTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCt                                                                       <  1:4991199/150‑1 (MQ=255)
                                                                              ttCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTg                                                                      >  2:2036190/1‑150 (MQ=255)
                                                                              ttCCACATCTGCCCCAACGAGCTTGGTGGGAACCCGAGGATGTACGTCAACCTGTTCGGGGCGCGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTg                                                                      <  2:4213284/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)
                                                                                 cactTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGttt                                                                   <  1:4674435/150‑2 (MQ=255)
                                                                                 cacaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTc                                                                   <  2:3087340/150‑1 (MQ=255)
                                                                                  acaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCg                                                                  >  2:634892/1‑150 (MQ=255)
                                                                                   caTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCgg                                                                 <  2:2854187/150‑1 (MQ=255)
                                                                                   caTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCgg                                                                 <  2:1752867/150‑1 (MQ=255)
                                                                                    ttCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCggt                                                                                                                >  1:2280614/2‑101 (MQ=255)
                                                                                    ttCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGTTTGGc                                                                >  1:1897215/2‑150 (MQ=255)
                                                                                    aTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                                                                    <  1:951448/114‑1 (MQ=255)
                                                                                       ttCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccac                                                             >  1:1398490/1‑150 (MQ=255)
                                                                                          aCCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCcaccac                                                          >  2:4040984/1‑150 (MQ=255)
                                                                                            cAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTc                                                        >  2:3918321/1‑150 (MQ=255)
                                                                                             aaCGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccacca                                                           <  2:2162925/146‑1 (MQ=255)
                                                                                             aaCGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccacca                                                           >  1:2162925/1‑146 (MQ=255)
                                                                                                     tggtggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTa                                               >  1:1220028/1‑150 (MQ=255)
                                                                                                      ggtggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTAc                                              <  2:2369501/150‑1 (MQ=255)
                                                                                                       gtggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACg                                             >  2:2523459/1‑150 (MQ=255)
                                                                                                       gtggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACg                                             >  2:5250499/1‑150 (MQ=255)
                                                                                                        tggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGc                                            >  2:1556058/1‑150 (MQ=255)
                                                                                                         gggAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCa                                           >  2:4893045/1‑150 (MQ=255)
                                                                                                          ggAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCAt                                          >  2:43499/1‑150 (MQ=255)
                                                                                                               ccGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgc                                     >  2:3892218/1‑150 (MQ=255)
                                                                                                                cGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAATAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcg                                    >  1:1407906/1‑150 (MQ=255)
                                                                                                                 gatgatGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcgc                                   >  1:218509/1‑150 (MQ=255)
                                                                                                                 gatgatGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcgc                                   >  1:2205247/1‑150 (MQ=255)
                                                                                                                 gatgatGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcgc                                   >  1:2874739/1‑150 (MQ=255)
                                                                                                                 gatgatGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcgc                                   >  2:933319/1‑150 (MQ=255)
                                                                                                                 gatgatGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgcgc                                   >  1:2874928/1‑150 (MQ=255)
                                                                                                                   tgatgTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGa                                 <  2:5166102/150‑1 (MQ=255)
                                                                                                                     atgTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATc                               <  2:3275273/150‑1 (MQ=255)
                                                                                                                        tACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGcc                            >  1:4139168/1‑150 (MQ=255)
                                                                                                                           gTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATc                         >  2:3760605/1‑150 (MQ=255)
                                                                                                                           gTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATc                         >  2:2602451/1‑150 (MQ=255)
                                                                                                                                   gTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAATAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCtt                 <  2:1407906/150‑1 (MQ=255)
                                                                                                                                    ttCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCttt                >  1:5135841/1‑150 (MQ=255)
                                                                                                                                       tGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCt             >  2:2424429/1‑150 (MQ=255)
                                                                                                                                        gggCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGAtt                                                 <  1:3471645/113‑1 (MQ=255)
                                                                                                                                        gggCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGAtt                                                 >  2:3471645/1‑113 (MQ=255)
                                                                                                                                        gggCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTGACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTg            >  2:2406986/1‑150 (MQ=255)
                                                                                                                                         ggCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGt           >  1:3064242/1‑150 (MQ=255)
                                                                                                                                         ggCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGt           >  1:4504641/1‑150 (MQ=255)
                                                                                                                                             tGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATgg       <  1:2602451/150‑1 (MQ=255)
                                                                                                                                              ggggTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGc      <  1:3974883/150‑1 (MQ=255)
                                                                                                                                                 gTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGc   >  2:1890281/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTATGGCTGCA  >  NC_002947/951875‑952170

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: