Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 951,979 T→C 15.3% G269G (GGT→GGC 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,9790TC15.3% 391.8 / 51.7 158G269G (GGT→GGCcyoBcytochrome bo terminal oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base T (71/62);  new base C (13/11);  total (85/73)
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.73e-01

Only 100 of 173 total aligned reads displayed.

ATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGC  >  NC_002947/951830‑952128
                                                                                                                                                     |                                                                                                                                                     
aTGAATCTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTtggt                                                                                                                                                       <  2:175706/150‑1 (MQ=255)
  gAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTtggtgg                                                                                                                                                     <  1:405070/150‑1 (MQ=255)
  gAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTtggtgg                                                                                                                                                     <  1:99577/150‑1 (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:2905103/150‑1 (MQ=255)
              cATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCgatgat                                                                                                                                         <  2:2905356/150‑1 (MQ=255)
                tGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGt                                                                                                                                       <  2:3467530/150‑1 (MQ=255)
                 gCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTa                                                                                                                                      >  1:984058/1‑150 (MQ=255)
                  ccGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTTCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTAc                                                                                                                                     <  2:3104402/150‑1 (MQ=255)
                  ccGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTAc                                                                                                                                     <  1:1160715/150‑1 (MQ=255)
                  ccGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTAc                                                                                                                                     <  1:985828/150‑1 (MQ=255)
                  ccGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTAc                                                                                                                                     >  2:4883245/1‑150 (MQ=255)
                   cGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACg                                                                                                                                    >  1:70520/1‑150 (MQ=255)
                   cGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACg                                                                                                                                    >  1:500153/1‑150 (MQ=255)
                        ttCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAAc                                                                                                                               <  1:3862139/150‑1 (MQ=255)
                          cACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCt                                                                                                                             >  2:5330419/1‑150 (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)
                                       aCCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTgg                                                                                                                <  1:3678091/150‑1 (MQ=255)
                                       aCCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTgg                                                                                                                >  1:2732816/1‑150 (MQ=255)
                                        ccTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTggg                                                                                                               >  1:3275273/1‑150 (MQ=255)
                                        ccTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTggg                                                                                                               <  1:2641634/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)
                                                    tgtTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAgg                                                                                                   >  1:5157657/1‑150 (MQ=255)
                                                      ctGATCGTGGCTGCGTTCCCGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGtc                                                                                                 <  1:627622/149‑2 (MQ=255)
                                                       tGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGttt                                                                                                <  2:4664259/150‑1 (MQ=255)
                                                        gATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTa                                                                                               >  1:2424045/1‑150 (MQ=255)
                                                         aTCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTAc                                                                                              >  1:1353924/1‑150 (MQ=255)
                                                         aTCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCATCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTAc                                                                                              >  2:1843666/1‑150 (MQ=255)
                                                            gtcgCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACacc                                                                                           <  1:37450/150‑3 (MQ=255)
                                                             gggCTGCGTTCCCGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTatattc                                                                                          <  2:1214651/148‑3 (MQ=255)
                                                                cTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcctga                                                                                       >  1:20954/1‑150 (MQ=255)
                                                                 ttCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcctgat                                                                                      >  2:1059496/1‑150 (MQ=255)
                                                                    gTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcctgatcct                                                                                   >  2:4118558/1‑150 (MQ=255)
                                                                    gTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTAAGTCAACCTGTTCTGGGCTTGGGGTCACCCTGAGGTTTACatcctgatact                                                                                   >  2:2527954/1‑150 (MQ=255)
                                                                    gTTCCCGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATattctggttct                                                                                   >  2:2717358/1‑145 (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                                                            >  2:4547175/1‑150 (MQ=255)
                                                                                           tCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCgg                                                            >  1:4642260/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)
                                                                                             gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                          <  1:1709056/150‑1 (MQ=255)
                                                                                             gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                          <  2:3955881/150‑1 (MQ=255)
                                                                                             gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                          <  2:5265913/150‑1 (MQ=255)
                                                                                              ctctGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGaaa                                                         <  1:5378846/150‑2 (MQ=255)
                                                                                               tctGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGt                                                        >  2:3815332/1‑150 (MQ=255)
                                                                                               tctGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGt                                                        <  1:2065649/150‑1 (MQ=255)
                                                                                                    tgACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                          <  1:3198844/143‑1 (MQ=255)
                                                                                                    tgACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACtt                                                   >  2:4323487/1‑150 (MQ=255)
                                                                                                       cTGTTGACCGTTATCTGGACTTCCACATCTTCACCATCGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGa                                                <  1:5313316/150‑1 (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)
                                                                                                          ttGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACATCGACGt                                             >  1:1151190/1‑150 (MQ=255)
                                                                                                            gACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTtc                                           >  2:353473/1‑150 (MQ=255)
                                                                                                               cGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTctct                                        >  2:1009657/1‑150 (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)
                                                                                                                    tCTGGATTTCCACTTCTTCACCAATGATGGTGGAGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCaa                                   >  2:2521844/1‑150 (MQ=255)
                                                                                                                         atttCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCgtc                              >  2:4544089/3‑148 (MQ=255)
                                                                                                                          cTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCt                             >  2:2440096/1‑150 (MQ=255)
                                                                                                                          cTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCt                             <  1:4991199/150‑1 (MQ=255)
                                                                                                                            tCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGTCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGt                           >  1:5197280/1‑150 (MQ=255)
                                                                                                                             ccACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGtt                          >  2:4194207/1‑150 (MQ=255)
                                                                                                                             ccACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGtt                          >  1:2266067/1‑150 (MQ=255)
                                                                                                                              cacaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTc                         <  2:3087340/150‑1 (MQ=255)
                                                                                                                               acaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGc                                                                      >  2:1713727/1‑104 (MQ=255)
                                                                                                                               acaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCg                        >  2:634892/1‑150 (MQ=255)
                                                                                                                                caTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCgg                       <  2:2854187/150‑1 (MQ=255)
                                                                                                                                caTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCgg                       <  2:1752867/150‑1 (MQ=255)
                                                                                                                                caTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCgg                       <  2:2854200/150‑1 (MQ=255)
                                                                                                                                 ttCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCggt                                                                      >  1:2280614/2‑101 (MQ=255)
                                                                                                                                 ttCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGTTTGGc                      >  1:1897215/2‑150 (MQ=255)
                                                                                                                                   cTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGcca                    <  2:242826/150‑1 (MQ=255)
                                                                                                                                    ttCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccac                   >  2:2196022/1‑150 (MQ=255)
                                                                                                                                     tCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccacc                  <  2:3925161/150‑1 (MQ=255)
                                                                                                                                          aaCGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccacca                 >  1:2162925/1‑146 (MQ=255)
                                                                                                                                                    gtggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACg   >  2:5250499/1‑150 (MQ=255)
                                                                                                                                                    gtggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACg   >  2:2523459/1‑150 (MQ=255)
                                                                                                                                                     tggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGc  >  2:1556058/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGAAACTGATGGACATGCCGATCTTCACCTGGACCTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGC  >  NC_002947/951830‑952128

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: