Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 952,014 C→T 8.3% A281V (GCG→GTG) ‡ 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,0140CT 8.3% 928.6 / 51.9 311A281V (GCG→GTG) ‡cyoBcytochrome bo terminal oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (153/132);  new base T (15/11);  total (168/143)
Fisher's exact test for biased strand distribution p-value = 8.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 316 total aligned reads displayed.

CTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTAT  >  NC_002947/951865‑952163
                                                                                                                                                     |                                                                                                                                                     
  gCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGt                                                                                                                                                     >  1:5291418/1‑150 (MQ=255)
    aCCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTgg                                                                                                                                                   <  1:184199/150‑1 (MQ=255)
          tcgaaCGTCCTGATCGTGGCTGCGTTCCCGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGacac                                                                                                                                             >  1:1195963/4‑150 (MQ=255)
           ccAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCAcc                                                                                                                                            <  2:4851047/150‑1 (MQ=255)
            cAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCAccc                                                                                                                                           >  1:1925066/1‑150 (MQ=255)
             aaCGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCt                                                                                                                                          <  2:8854102/150‑1 (MQ=255)
             aaCGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCt                                                                                                                                          <  1:8291735/150‑1 (MQ=255)
               cGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGa                                                                                                                                        <  1:375731/150‑1 (MQ=255)
                 tgtTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAgg                                                                                                                                      <  2:125806/150‑1 (MQ=255)
                    tGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGttt                                                                                                                                   <  1:5702670/150‑1 (MQ=255)
                      aTCGTGGCTGCGTTCCCGATCCTGACCGTGACGCTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTat                                                                                                                                 <  1:5983972/150‑2 (MQ=255)
                        cgtcgCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACat                                                                                                                               <  1:4797259/150‑1 (MQ=255)
                         gtcgCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatc                                                                                                                              <  1:2800332/150‑1 (MQ=255)
                          tcgCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcc                                                                                                                             >  1:2940116/1‑150 (MQ=255)
                              ttCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcctgat                                                                                                                         >  1:2371795/1‑150 (MQ=255)
                              ttCGTTCCCGATCCTGACTGCTGCACTCGCTCAGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTCTTCTGGGCGTGGGCTCACCCTGAGGTTTACatcctgat                                                                                                                         >  1:8728909/1‑150 (MQ=255)
                                   tCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGACCGCGTTGGTGTACGTCAACCTGTTCTGGGCGTGGGGTCTGCCTGAGGTTGGCCTGCTGAGCCTGc                                                                                                                    >  2:7569998/1‑150 (MQ=255)
                                      cGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCgg                                                                                                                 >  1:2754149/1‑150 (MQ=255)
                                      cGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCgg                                                                                                                 >  1:8084524/1‑150 (MQ=255)
                                      cGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCgg                                                                                                                 >  2:5733600/1‑150 (MQ=255)
                                            tGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCg                                                                                                           <  1:315532/150‑1 (MQ=255)
                                             gACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcc                                                                                                                             >  2:2118485/1‑131 (MQ=255)
                                             gACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACatcc                                                                                                                             <  1:2118485/131‑1 (MQ=255)
                                             gACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCgg                                                                                                          >  1:8465022/1‑150 (MQ=255)
                                              aCTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGc                                                                                                         <  1:2183406/150‑1 (MQ=255)
                                                tgctgcACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCgt                                                                                                       <  1:7916727/150‑1 (MQ=255)
                                                 gctgcACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCgtg                                                                                                      >  1:1798915/1‑150 (MQ=255)
                                                  tgacgcTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGgtgttt                                                                                                     >  1:6955801/6‑148 (MQ=255)
                                                  tgacgcTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGgtgttt                                                                                                     <  2:1830407/145‑3 (MQ=255)
                                                  tgacgcTCGCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGgtgttt                                                                                                     >  2:3705928/6‑148 (MQ=255)
                                                   tgcACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGtt                                                                                                    <  1:3466373/150‑1 (MQ=255)
                                                     cACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTtct                                                                                                  >  1:7151703/1‑150 (MQ=255)
                                                      aCTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTtctc                                                                                                 <  2:4287660/150‑1 (MQ=255)
                                                      aCTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTtctc                                                                                                 <  1:3479222/150‑1 (MQ=255)
                                                       cTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCg                                                                                                >  1:9131732/1‑150 (MQ=255)
                                                       cTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCg                                                                                                <  2:7045349/150‑1 (MQ=255)
                                                        tCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCgg                                                                                               <  1:7469054/150‑1 (MQ=255)
                                                         cGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGa                                                                                              <  2:7441693/150‑1 (MQ=255)
                                                          gCTCTGCTGTCTCTTGATCGTTATCTGGATTTCCACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGaa                                                                                             >  1:7086764/1‑150 (MQ=255)
                                                          gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                                                             <  1:306845/150‑1 (MQ=255)
                                                          gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                                                             <  2:1521966/150‑1 (MQ=255)
                                                          gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                                                             <  1:1476902/150‑1 (MQ=255)
                                                          gCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGaa                                                                                             <  2:8654710/150‑1 (MQ=255)
                                                              tgctgACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTaa                                                                                         >  2:3850796/1‑150 (MQ=255)
                                                              tgctgACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTaa                                                                                         <  1:3818615/150‑1 (MQ=255)
                                                              tgctgACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTaa                                                                                         >  2:914337/1‑150 (MQ=255)
                                                                 tgACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACtt                                                                                      >  1:2429993/1‑150 (MQ=255)
                                                                  gACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTc                                                                                     <  1:8604444/150‑1 (MQ=255)
                                                                   aCTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCg                                                                                    >  1:1747970/1‑150 (MQ=255)
                                                                       ttGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGt                                                                                >  1:4309648/1‑150 (MQ=255)
                                                                         gACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCATCCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTtc                                                                              >  2:2581616/1‑150 (MQ=255)
                                                                          aCCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTtct                                                                             >  1:8176111/1‑150 (MQ=255)
                                                                                  cTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCaaa                                                                     >  1:7268620/1‑150 (MQ=255)
                                                                                   tGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAAc                                                                    >  2:4120387/1‑150 (MQ=255)
                                                                                   tGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAAc                                                                    >  1:4097250/1‑150 (MQ=255)
                                                                                   tGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAAc                                                                    >  1:5383169/1‑150 (MQ=255)
                                                                                    ggACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACg                                                                   >  2:6737069/1‑150 (MQ=255)
                                                                                       cTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCt                                                                >  1:368390/1‑150 (MQ=255)
                                                                                       cTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCt                                                                <  1:6086960/150‑1 (MQ=255)
                                                                                        ttCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTg                                                               >  1:5226293/1‑150 (MQ=255)
                                                                                          ccACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGtt                                                             >  1:5823767/1‑150 (MQ=255)
                                                                                          ccACTTCTTCACCAATGATGGTGGCGGCAATCCGATGATGTACGTCAACCTCATCTGGATCTGGGGACACCCGGAAGTCTATATTCTGGTTCTGCCTTGCTTCGGTGTTTTCTCCGAAATCATTGCGACCTTTTCGGGCAAGCGTCTGtt                                                             >  1:6491334/1‑150 (MQ=255)
                                                                                          ccACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGtt                                                             >  1:1162908/1‑150 (MQ=255)
                                                                                          ccACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGtt                                                             >  1:3426442/1‑150 (MQ=255)
                                                                                           cacaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTc                                                            >  1:4448746/1‑150 (MQ=255)
                                                                                            acaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCg                                                           >  1:662962/1‑150 (MQ=255)
                                                                                            acaTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCg                                                           >  2:3761515/1‑150 (MQ=255)
                                                                                              aTCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGc                                                         <  2:7287500/150‑1 (MQ=255)
                                                                                               tCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGcc                                                        >  1:8537966/1‑150 (MQ=255)
                                                                                               tCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGcc                                                        >  1:7959952/1‑150 (MQ=255)
                                                                                                 ttCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccac                                                      >  1:1295457/1‑150 (MQ=255)
                                                                                                  tCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccacc                                                     >  2:2259552/1‑150 (MQ=255)
                                                                                                   cACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGccacca                                                    >  2:8977048/1‑150 (MQ=255)
                                                                                                      cAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTc                                                 >  2:6903519/1‑150 (MQ=255)
                                                                                                      cAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTc                                                 >  1:6589700/1‑150 (MQ=255)
                                                                                                      cAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTc                                                 <  1:3004351/150‑1 (MQ=255)
                                                                                                        aCGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCga                                               <  2:6221986/150‑1 (MQ=255)
                                                                                                         cGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCgat                                              <  2:3525917/150‑1 (MQ=255)
                                                                                                             cTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGAtt                                          <  1:944035/150‑1 (MQ=255)
                                                                                                              ttGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGa                                                                                   >  2:8103425/1‑108 (MQ=255)
                                                                                                               tggtggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTa                                        <  2:1706642/150‑1 (MQ=255)
                                                                                                                  tggGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGc                                     >  2:8932995/1‑150 (MQ=255)
                                                                                                                      aaCCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCg                                 <  1:7566480/150‑1 (MQ=255)
                                                                                                                      aaCCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCg                                 <  2:6052981/150‑1 (MQ=255)
                                                                                                                         ccGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGgc                              <  2:2715796/150‑1 (MQ=255)
                                                                                                                               atgTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATc                        <  1:3971492/150‑1 (MQ=255)
                                                                                                                               atgTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATc                        <  2:4777118/150‑1 (MQ=255)
                                                                                                                                  tACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGcc                     <  2:7830188/150‑1 (MQ=255)
                                                                                                                                   aCGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCa                    <  1:5011209/150‑1 (MQ=255)
                                                                                                                                    cGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCAt                   <  2:8271566/150‑1 (MQ=255)
                                                                                                                                      tCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAATCGCCTGTTCGGCGACCACTCGTTGATTTACGCATCGGGCgcgatcgcgatcg                 >  2:5629713/1‑149 (MQ=255)
                                                                                                                                         aCCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCg              >  1:5534322/1‑150 (MQ=255)
                                                                                                                                          ccTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCCGATCCTGCCGGCCTTCGGCGTCTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCgg             <  2:3056942/150‑1 (MQ=255)
                                                                                                                                           cTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGc            >  1:3603606/1‑150 (MQ=255)
                                                                                                                                            tGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCt           <  2:4635811/150‑1 (MQ=255)
                                                                                                                                                cTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGc       <  2:855871/150‑1 (MQ=255)
                                                                                                                                                 tGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCTACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCt      >  1:5142767/1‑150 (MQ=255)
                                                                                                                                                 tGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCt      >  2:6586583/1‑150 (MQ=255)
                                                                                                                                                 tGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCt      >  2:5032995/1‑150 (MQ=255)
                                                                                                                                                   ggCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGt    >  1:2241700/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGCACCTGGGCCAACGTGTTGATCGTCGCTTCGTTCCCGATCCTGACTGCTGCACTCGCTCTGCTGACTGTTGACCGTTATCTGGACTTCCACATCTTCACCAACGAGCTTGGTGGGAACCCGATGATGTACGTCAACCTGTTCTGGGCGTGGGGTCACCCTGAGGTTTACATCCTGATCCTGCCGGCCTTCGGCGTGTTCTCGGAAGTAACTTCGACGTTCTCTGGCAAACGCCTGTTCGGCCACCACTCGATGATTTACGCATCGGGCGCGATCGCCATCCTCGGCTTTGCTGTAT  >  NC_002947/951865‑952163

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: