Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,696,120 A→C 13.7% K26T (AAG→ACG)  PP_1487 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,696,1200AC13.7% 845.6 / 13.5 313K26T (AAG→ACG) PP_1487hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (187/71);  new base C (2/52);  total (189/124)
Fisher's exact test for biased strand distribution p-value = 1.99e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 320 total aligned reads displayed.

CAAATCCCTGCATGGCGTGTGCTCCTCCAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGCCTCGCCGC  >  NC_002947/1695971‑1696267
                                                                                                                                                     |                                                                                                                                                   
cAAATCCCTGCATGGCGTGTGCTCCTCCAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCaa                                                                                                                                                     <  1:4488364/150‑1 (MQ=255)
cAAATCCCTGCATGGCGTGTGCTCCTCCAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCaa                                                                                                                                                     >  2:2894558/1‑150 (MQ=255)
cAAATCCCTGCATGGCGTGTGCTCCTCCAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCaa                                                                                                                                                     >  2:2977486/1‑150 (MQ=255)
  aaTCCCTGCATGGCGTGTGCTCCTCCAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGc                                                                                                                                                   <  2:1711167/150‑1 (MQ=255)
   aTCCCTGCATGGCGTGTGCTCCTCCAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGcc                                                                                                                                                  >  1:4859904/1‑150 (MQ=255)
     cccTGCATGGCGTGTGCTCCTCCAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTg                                                                                                                                                >  2:2439836/1‑150 (MQ=255)
      ccTGCATGGCGTGTGCTCCTCCAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTgg                                                                                                                                               <  1:5314119/150‑1 (MQ=255)
       cTGCATGGCGTGTGCTCCTCCAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGc                                                                                                                                              >  1:4783094/1‑150 (MQ=255)
             ggCGTGTGCTCCTCCAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCt                                                                                                                                        >  2:3301326/1‑150 (MQ=255)
              gCGTGTGCTCCTCCAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCtt                                                                                                                                       >  1:6487450/1‑150 (MQ=255)
                gtgtGCTCCTCCAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTcg                                                                                                                                     >  2:1474224/1‑150 (MQ=255)
                gtgtGCTCCTCCAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTcg                                                                                                                                     <  2:4374287/150‑1 (MQ=255)
                    gCTCCTCCAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACACGTGGCTTTCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCGCGCCTGGCAGCCCTTcgggcg                                                                                                                                 >  1:5647684/1‑150 (MQ=255)
                      tcctccAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTcgggcggg                                                                                                                               >  1:5131011/1‑150 (MQ=255)
                        ctccAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCgggcgggtt                                                                                                                             >  2:5626364/1‑150 (MQ=255)
                        ctccAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCgggcgggtt                                                                                                                             >  1:81103/1‑150 (MQ=255)
                         tccAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAGAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCgggcggggtg                                                                                                                            >  2:6862712/1‑150 (MQ=255)
                          ccAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCgggcgggttgg                                                                                                                           >  2:2147020/1‑150 (MQ=255)
                           cAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCgggcgggttggg                                                                                                                          >  1:388269/1‑150 (MQ=255)
                           cAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCgggcgggttggg                                                                                                                          >  1:428895/1‑150 (MQ=255)
                            aaGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCgggcgggttgggc                                                                                                                         >  2:5046827/1‑150 (MQ=255)
                            aaGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCgggcgggttgggc                                                                                                                         >  1:5367053/1‑150 (MQ=255)
                              gCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCgggcgggttgggcgg                                                                                                                       <  1:5409131/150‑1 (MQ=255)
                                    ctctGGGGTTTTAATCTTCGTACGTGGAGACACGCGACTGAACACCCGGGGCTTGCTCGATCAACTCCACAAATCCGGCCCGCAACTGCTCGAAAAGCAGGGGGGCGCCACAAAGCCTGGCAGCCATTCGGGCGGGTTGGGCGGgttgtt                                                                                                                 <  2:3663407/150‑1 (MQ=255)
                                    ctctGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGgttgtt                                                                                                                 <  2:2770864/150‑1 (MQ=255)
                                       tGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTc                                                                                                              >  1:3826803/1‑150 (MQ=255)
                                         gCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAAGCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTGGGGCGGGTTGGGCGGGTTGTTGTCCg                                                                                                            <  2:3535344/150‑1 (MQ=255)
                                          cGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCgg                                                                                                           >  1:1088848/1‑150 (MQ=255)
                                           ttttGAATCTTCTTACGTGTACACACGCTAATCAACACCCGTGGATTGCTCTATAAACTCCTCAAATCCGCCAAGAAACTGATCGAAAAGATGGGGGGCGCAACAAAGCCTGGCAGCAATTCGGGCGTTTTGGGCGGGTTGTTGTCCggg                                                                                                          <  2:6606762/149‑1 (MQ=255)
                                               taaTCTTCTTAGGTGGAGACACGCTATTCAACACCGGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCTAAAAGCTGGGGGGCGCCACCACGCGGGGCAGCCCTTGGGGGGGTTTGGGGGGGTTGTTGTCCGGGGCAg                                                                                                      <  1:3238537/149‑1 (MQ=255)
                                               gAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAg                                                                                                      >  1:137703/1‑150 (MQ=255)
                                               gAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAg                                                                                                      >  2:4002973/1‑150 (MQ=255)
                                               gAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAg                                                                                                      >  1:2840923/1‑150 (MQ=255)
                                                aaTCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAgg                                                                                                     >  2:6710129/1‑150 (MQ=255)
                                                    ttcttACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGggc                                                                                                 <  2:1115667/150‑1 (MQ=255)
                                                     tcttCCGTGGAGACACGCTAATGAACCCCCGTGGCTTGCTCGATCAACTCCCCAACTCCGGCCAGCAACTGCTCGAAAAGCCGGGGGGCGCCACCACGCCGGGCAGCCCTTCGGGGGGGTGGGGCGGGTTGTTGTCCGGGGCAAGGggcg                                                                                                <  1:6638724/150‑1 (MQ=255)
                                                       ttACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGggcggc                                                                                              >  1:4627132/1‑150 (MQ=255)
                                                            tGGAGACACGCTAATTAACACCCGTGGCTTGCTCGATCAACTCCTCAAACCCGGCCAGAAACTGCTCGAAAAGCTGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGGGGGTTGGGGGGGTTGTTGTCCGGGGCAGGGGGCGGCctgct                                                                                         <  2:3201963/150‑1 (MQ=255)
                                                             ggAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCctgctg                                                                                        <  1:6813120/150‑1 (MQ=255)
                                                             ggAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCctgctg                                                                                        >  2:827521/1‑150 (MQ=255)
                                                                  caggcgATTGAACACCCGTGGCTGGCTCGATCAACTCCTCAAATCCGGCCAGCAAGTGCTCTAAAAGCTGGGGGGCGCCACCACGCCTGGCAGCCCTTGGGGCGGTTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCtgggtgg                                                                                   <  2:5827704/147‑1 (MQ=255)
                                                                   acGCGAATGAACACCCTTGACTTGCTCGATAAACTCCTAAAATCCGGCCAGCAACTGCTCTAAAAGCAGGGGGGCGCCACCACGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCtgggtggg                                                                                  <  1:2052723/150‑1 (MQ=255)
                                                                   acGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCtgggtggg                                                                                  >  2:3831273/1‑150 (MQ=255)
                                                                   acGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCtgggtggg                                                                                  >  1:1756718/1‑150 (MQ=255)
                                                                   acGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCtgggtggg                                                                                  >  2:1020832/1‑150 (MQ=255)
                                                                   acGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCtgggtggg                                                                                  >  1:1238782/1‑150 (MQ=255)
                                                                    cgcgAATGAACACCGGTGGGTTGCCCGACCAACCCCCCAAACCCGGCCCCCACCTGCGCCAAACGCGGGGGGGCCCCACCACCCCGGGCAGCCCTTCGGGGGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTgggg                                                                                 <  2:1242068/150‑1 (MQ=255)
                                                                    cgcgAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTgggg                                                                                 <  2:1669696/150‑1 (MQ=255)
                                                                     gctatttaaCACCCGTGGCTTGCTAAAACAACCCCCAAAATCCGCCCCCCAACTGCCCTAAAACCCGGGGGGCGCCACCAAGCCTGGCAGCCCTTGGGGGGGTTTGGGGGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTggggg                                                                                <  2:3037941/143‑1 (MQ=255)
                                                                     gcgAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTggggg                                                                                >  2:3594308/1‑150 (MQ=255)
                                                                       gAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTgggg                                                                                 <  1:1096805/147‑1 (MQ=255)
                                                                       gAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGgcg                                                                              >  1:371440/1‑150 (MQ=255)
                                                                       gAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGgcg                                                                              >  2:2811893/1‑150 (MQ=255)
                                                                          tGAACCCCCGGGGCTTGCTCGAACAACTCCTCAAAGCCGGCCAGCACCGGCGCGAAAAGCTGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCAAGGTgggggtgatc                                                                           <  2:4624049/150‑6 (MQ=255)
                                                                          tGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGctc                                                                           >  2:5261771/1‑150 (MQ=255)
                                                                           gaccacCCGGGGCTTGCACGAACAACTCCCCAAACCCGGCCAGCAACGGCTCGAAACGCTGGGGGGCGCCACCACGCCTGGCCGCCCTGCGGGGGGGTGGGGGGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGctct                                                                          <  2:407234/147‑1 (MQ=255)
                                                                            acccccggTGGGTTGCTCGAACAACTCCTAAACGCCGCCCCGCAACTGCTCGAAAAGCCGGGGGGCGCCACCACGCCGGGCCGCCCTTGGGGGGGGTTGGGGGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTg                                                                         <  2:3428194/143‑1 (MQ=255)
                                                                            aaCACCCGTGGCTTGCTCGATCAACTCCCCAAATCCGGCCAGCAACTGCTCTAAAAGCCGGGGGGCGCCACCAAGCCTGGCCGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTg                                                                         <  1:1560719/150‑1 (MQ=255)
                                                                             acacCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTgg                                                                        >  1:2291530/1‑150 (MQ=255)
                                                                                   gTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGt                                                                  >  2:3937630/1‑150 (MQ=255)
                                                                                    tGGTTTGCTCTATAAACTCCTAAAATCGGGCCAGAAAGTGCTCTAAATGCTGGGGGGCGCCACAACGCCGGGCAGCCCTTGGGGGGGTTTGGGGGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGtt                                                                 <  1:4271268/150‑1 (MQ=255)
                                                                                       gtggcgcGACCAACGCCCCAAATCCGGCCCGCAACTGCTCGAAACGCGGGGGGGCGCCACCACGCCGGGCCGCCGTGGGGGGGGGTTGGGGGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCt                                                              <  2:2575250/144‑1 (MQ=255)
                                                                                        tggCTCTAACAACTCCACAAACCCGGCCACCAAGTGCTCCAAAAGCGGGGGGGCGCCACCACGCGTGGCAGCCCTGGGGGGGGTTGGGGGGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTg                                                             <  1:6046966/148‑1 (MQ=255)
                                                                                            tCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTa                                                         >  2:2507006/1‑150 (MQ=255)
                                                                                            tCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTa                                                         <  2:4569825/150‑1 (MQ=255)
                                                                                               aacAACCCCTCAAACCCGGCCAGCACCGGCTCTAAAAGCCGGGGGGCGCCACCACGCCTGGCCGCCCTTGGGGGGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCa                                                      <  2:6441263/148‑1 (MQ=255)
                                                                                                   aCTCCCCAAACCCGGCCCGCACGGGCTCGAAACGCGGGGGGGCGCCACCACGCGTGGCCGCCCTTGGGGCGGGTTGGGGGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCaagaa                                                  <  2:3184040/150‑1 (MQ=255)
                                                                                                    ctcctcAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCCGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCaagaag                                                 >  2:4293209/1‑150 (MQ=255)
                                                                                                     tcctcAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGTGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCGGGGTAGCAAGAAgg                                                >  2:6100665/1‑150 (MQ=255)
                                                                                                     tcctcAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAgg                                                >  2:4821375/1‑150 (MQ=255)
                                                                                                      cctcAACTCCGGCCTGCAACTGCTCGAAAAGCTGGGGGGCGCCACCAAGCGGGGCAGCCCTTGGGGCGGGTTGGGGGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGc                                               <  1:2147020/150‑1 (MQ=255)
                                                                                                       ctcAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGcc                                              >  2:1655221/1‑150 (MQ=255)
                                                                                                       ctcAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGcc                                              >  1:4300705/1‑150 (MQ=255)
                                                                                                       ctcAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGcc                                              >  1:5740669/1‑150 (MQ=255)
                                                                                                        tcAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGccc                                             <  2:5661459/150‑1 (MQ=255)
                                                                                                        tcAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGccc                                             >  2:605870/1‑150 (MQ=255)
                                                                                                        tcAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGccc                                             >  2:6283916/1‑150 (MQ=255)
                                                                                                        tcAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGccc                                             >  1:914392/1‑150 (MQ=255)
                                                                                                         cAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCg                                            >  1:600508/1‑150 (MQ=255)
                                                                                                            agccGGCCTCCACCTGCTCCAAAAGCTGGGGGGCGCAACAACGCGTGGCAGCCCTTGGGGGGGGTTGGGGGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCaa                                         <  1:4061264/148‑1 (MQ=255)
                                                                                                            aTCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCaa                                         >  2:4931859/1‑150 (MQ=255)
                                                                                                             tCCGGCCAGCAACGGCTCGAAAAGCCGGGGGGCGCCACCCAGCCGGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAg                                        <  2:2062447/150‑1 (MQ=255)
                                                                                                                    cgCAACGGCGCTAAACGCGGGGGGGCGCCCCCCCGCGGGGCCGCCGTGGGGGGGGGTGGGGGGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTg                                 <  2:1402090/149‑1 (MQ=255)
                                                                                                                    agaaaCTGCTCTAAAAGCTGGGGGGCGCCACCACGCCGGGCAGCCCTTGGGGGGGTTTGGGGGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTg                                 <  1:1430469/147‑1 (MQ=255)
                                                                                                                     gCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTgg                                >  2:1766416/1‑150 (MQ=255)
                                                                                                                       aaCTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCa                              <  2:3317936/150‑1 (MQ=255)
                                                                                                                           gCTCTAAAAGCCGGGGGGCGCCACCACGCCTGGCAGCCCTCGGGGGGGGTTGGGGGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGc                          <  1:4193910/150‑1 (MQ=255)
                                                                                                                            ctctAAATGCTGGGGGGGGCCACCATGCCTGGCTGCCCTTGGGGGGGGTTGGGGGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGcc                         <  1:4437958/150‑1 (MQ=255)
                                                                                                                            cTCGAAAGGCTGGGGGGCGCCACCACGCCGGGCCGCCGTTCGGGGGGGTGGGGGGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGcc                         <  2:6742161/150‑1 (MQ=255)
                                                                                                                             tCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGccc                        >  1:2013844/1‑150 (MQ=255)
                                                                                                                             tCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGccc                        >  1:5850962/1‑150 (MQ=255)
                                                                                                                             tCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGccc                        >  1:5092060/1‑150 (MQ=255)
                                                                                                                               gAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTg                      >  1:3540631/1‑150 (MQ=255)
                                                                                                                               gAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTg                      >  1:2532448/1‑150 (MQ=255)
                                                                                                                                 aaaGCCGGGGGGCGCCACCACGCCGGGCTGCCCTTGGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGAc                    <  1:4699308/150‑1 (MQ=255)
                                                                                                                                  aaGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACt                   >  1:2303659/1‑150 (MQ=255)
                                                                                                                                        gggggCGCCACCACGCCGGGCCGCCCTTGGGGGGGGTGGGGGGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGTAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTAcggc             <  2:6026867/150‑1 (MQ=255)
                                                                                                                                           ggCGCCACCACGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACggcggc          <  2:4692376/150‑1 (MQ=255)
                                                                                                                                              gcCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGCCTc       >  2:3113695/1‑150 (MQ=255)
                                                                                                                                                caccaAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGT‑GGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGCCTcgcc    >  2:2379102/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAAATCCCTGCATGGCGTGTGCTCCTCCAAGCCTTACTCTGGCGTTTGAATCTTCTTACGTGGAGACACGCGAATGAACACCCGTGGCTTGCTCGATCAACTCCTCAAATCCGGCCAGCAACTGCTCGAAAAGCAGGGGGGCGCCACCAAGCCTGGCAGCCCTTCGGGCGGGTTGGGCGGGTTGTTGTCCGGGGCAGGGGGCGGCCTGCTGGGTGGGGGCGCTCTGGGCCTGTTGCTGGGTAGCAAGAAGGCCCGCAAGTACGGTGGCAAGGCCCTGACTTACGGCGGCCTCGCCGC  >  NC_002947/1695971‑1696267

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑