Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,187,851 G→T 5.7% G115V (GGC→GTC)  PP_4562 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,187,8510GT 5.7% 1060.9 / 38.0 316G115V (GGC→GTC) PP_4562hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (143/154);  new base T (15/3);  total (158/158)
Fisher's exact test for biased strand distribution p-value = 6.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 319 total aligned reads displayed.

TGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGA  >  NC_002947/5187702‑5188000
                                                                                                                                                     |                                                                                                                                                     
tGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTgg                                                                                                                                                       <  2:3639647/150‑1 (MQ=255)
 gAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGc                                                                                                                                                      <  1:6896253/150‑1 (MQ=255)
   gCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGa                                                                                                                                                    >  1:2499615/1‑150 (MQ=255)
      gcagcaCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGagca                                                                                                                                                 <  1:6043527/150‑1 (MQ=255)
          caCCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTg                                                                                                                                             <  2:7873578/150‑1 (MQ=255)
             cTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGAc                                                                                                                                          <  1:1566165/150‑1 (MQ=255)
               gAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACtt                                                                                                                                        >  1:1109835/1‑150 (MQ=255)
                 aGGTGCCCACCAAGTTCCTCACGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGCCCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTcc                                                                                                                                      <  1:5060286/150‑1 (MQ=255)
                   gTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGc                                                                                                                                    >  2:4065670/1‑150 (MQ=255)
                        caccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACtt                                                                                                                               >  2:4259588/1‑150 (MQ=255)
                        caccaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACtt                                                                                                                               >  2:7209253/1‑150 (MQ=255)
                         accaATTTCCTCAAGGCTCACCTTGAGCATATGGTAGTGCGCAAGCCTGGCAAGCGCCGCTGCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCTCAACTTc                                                                                                                              <  2:5798572/150‑1 (MQ=255)
                         accaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTc                                                                                                                              >  1:2227835/1‑150 (MQ=255)
                         accaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTc                                                                                                                              <  1:1359058/150‑1 (MQ=255)
                         accaATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTc                                                                                                                              >  2:8036629/1‑150 (MQ=255)
                           ccaTTTGCTCAAGGTGAGCCTTGAGGAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGt                                                                                                                            >  1:2642257/3‑150 (MQ=255)
                           caATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGt                                                                                                                            >  1:57125/1‑150 (MQ=255)
                           caATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGt                                                                                                                            <  2:312735/150‑1 (MQ=255)
                           caATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGt                                                                                                                            >  1:101550/1‑150 (MQ=255)
                            aaTTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTa                                                                                                                           >  1:2179711/1‑150 (MQ=255)
                            aaTTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTa                                                                                                                           <  1:8662173/150‑1 (MQ=255)
                               ttCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAg                                                                                                                        <  2:455101/150‑1 (MQ=255)
                                  cTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCg                                                                                                                     >  1:4672508/1‑150 (MQ=255)
                                  cTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCg                                                                                                                     >  2:5044355/1‑150 (MQ=255)
                                   tCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGTTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGt                                                                                                                    >  1:3995677/1‑150 (MQ=255)
                                   tCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGt                                                                                                                    >  2:2412403/1‑150 (MQ=255)
                                   tCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGt                                                                                                                    >  1:8398665/1‑150 (MQ=255)
                                   tCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGt                                                                                                                    >  1:3273444/1‑150 (MQ=255)
                                      aggtgagccTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAg                                                                                                                 >  2:4364316/8‑150 (MQ=255)
                                      aGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAg                                                                                                                 <  2:7313067/150‑1 (MQ=255)
                                         cTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGcc                                                                                                              >  1:7151709/1‑150 (MQ=255)
                                            aCCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTc                                                                                                           <  1:5918754/150‑1 (MQ=255)
                                              cTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAg                                                                                                         >  2:5453897/1‑150 (MQ=255)
                                               ttGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAgg                                                                                                        <  2:1058156/150‑1 (MQ=255)
                                                 gAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCt                                                                                                      >  2:7591319/1‑150 (MQ=255)
                                                   gCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCttt                                                                                                    <  1:2576761/150‑1 (MQ=255)
                                                    cAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTtc                                                                                                   <  1:4217428/150‑1 (MQ=255)
                                                      aaTGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTtctc                                                                                                 <  1:5472053/150‑1 (MQ=255)
                                                      aaTGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTtctc                                                                                                 <  1:8728030/150‑1 (MQ=255)
                                                       aTGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCg                                                                                                >  2:3917257/1‑150 (MQ=255)
                                                        tGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGa                                                                                               <  1:3747506/150‑1 (MQ=255)
                                                        tGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGa                                                                                               <  1:7648238/150‑1 (MQ=255)
                                                        tGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGa                                                                                               <  1:1885698/150‑1 (MQ=255)
                                                         ggAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGAt                                                                                              <  2:1130250/150‑1 (MQ=255)
                                                          gAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATg                                                                                             >  1:2384019/1‑150 (MQ=255)
                                                           aaGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGc                                                                                            >  1:5694298/1‑150 (MQ=255)
                                                             gTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGa                                                                                          <  1:985266/150‑1 (MQ=255)
                                                              tGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGAt                                                                                         >  2:6217406/1‑150 (MQ=255)
                                                               gcgcACGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATc                                                                                        <  1:3033740/150‑1 (MQ=255)
                                                               gcgcAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGACGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATc                                                                                        <  2:2696671/150‑1 (MQ=255)
                                                                 gcAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGt                                                                                      >  1:3659064/1‑150 (MQ=255)
                                                                      ccTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATgg                                                                                 <  2:5597763/150‑1 (MQ=255)
                                                                             aGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGa                                                                          >  1:6586846/1‑150 (MQ=255)
                                                                              gcgcCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGAt                                                                         <  2:4723894/150‑1 (MQ=255)
                                                                              gcgcCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGAt                                                                         >  2:2798578/1‑150 (MQ=255)
                                                                               cgccgcTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATCCATGACGATc                                                                        <  2:7006418/150‑1 (MQ=255)
                                                                                gccgcTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCa                                                                       >  2:646381/1‑150 (MQ=255)
                                                                                 ccgcTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAg                                                                      <  1:3957460/150‑1 (MQ=255)
                                                                                   gcTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAggg                                                                    <  2:7057120/150‑2 (MQ=255)
                                                                                    cTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGcc                                                                   <  2:6036726/150‑1 (MQ=255)
                                                                                            tCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGa                                                           >  2:8498623/1‑150 (MQ=255)
                                                                                            tCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGa                                                           <  1:5261363/150‑1 (MQ=255)
                                                                                             cAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGaa                                                          <  2:6801502/150‑1 (MQ=255)
                                                                                                  tGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAAcc                                                     >  1:6389116/1‑150 (MQ=255)
                                                                                                     cccTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGc                                                  <  2:35541/150‑1 (MQ=255)
                                                                                                               cGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCctt                                        <  1:8036629/150‑1 (MQ=255)
                                                                                                                  ccGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCcttctt                                     <  2:1729617/150‑1 (MQ=255)
                                                                                                                  ccGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCcttctt                                     <  1:7951094/150‑1 (MQ=255)
                                                                                                                  ccGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCcttctt                                     <  2:3283709/150‑1 (MQ=255)
                                                                                                                   cGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTg                                    >  2:7814256/1‑150 (MQ=255)
                                                                                                                   cGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTg                                    >  1:889766/1‑150 (MQ=255)
                                                                                                                    gggACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCATCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGa                                   <  2:126687/150‑1 (MQ=255)
                                                                                                                     ggACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGaa                                  <  1:4005824/150‑1 (MQ=255)
                                                                                                                     ggACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGaa                                  <  1:2046497/150‑1 (MQ=255)
                                                                                                                      gACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCGTCTTGAAg                                 >  1:7455476/1‑150 (MQ=255)
                                                                                                                         tCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTg                              >  1:93873/1‑150 (MQ=255)
                                                                                                                            cagctgctCGCCACCCGGACGGGTGTCGCGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGccc                           >  2:5342725/3‑150 (MQ=255)
                                                                                                                             tgctgcGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCTTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCt                          >  1:5393562/1‑150 (MQ=255)
                                                                                                                               ctgcTCGCCACCCGGACGGGTGTCGCGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTg                        >  2:5248231/1‑150 (MQ=255)
                                                                                                                               ctgcGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTg                        <  1:4882527/150‑1 (MQ=255)
                                                                                                                               ctgcGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTg                        <  1:3489529/150‑1 (MQ=255)
                                                                                                                                tgcTCGCCACCCGGACGGGTGTCGCGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGc                       >  1:6439473/1‑150 (MQ=255)
                                                                                                                                    cgACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGAc                   <  2:1568362/150‑1 (MQ=255)
                                                                                                                                       cacCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGtt                <  2:1835262/150‑1 (MQ=255)
                                                                                                                                       cacCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGtt                <  2:8301907/150‑1 (MQ=255)
                                                                                                                                        acCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTg               <  1:4430673/150‑1 (MQ=255)
                                                                                                                                        acCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTg               >  2:41589/1‑150 (MQ=255)
                                                                                                                                          ccGTCCGGTTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGTCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGaa             <  1:7209474/150‑1 (MQ=255)
                                                                                                                                          ccGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGaa             >  2:1522843/1‑150 (MQ=255)
                                                                                                                                          ccGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGaa             <  2:3041703/150‑1 (MQ=255)
                                                                                                                                           cGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAAc            >  1:182521/1‑150 (MQ=255)
                                                                                                                                           cGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAAc            <  2:4214008/150‑1 (MQ=255)
                                                                                                                                            gTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACa           <  2:8746440/150‑1 (MQ=255)
                                                                                                                                            gTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACa           <  1:8565281/150‑1 (MQ=255)
                                                                                                                                              ccGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGc         >  2:1376241/1‑150 (MQ=255)
                                                                                                                                                gtgTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCtt       >  2:628724/3‑150 (MQ=255)
                                                                                                                                                  gTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGc     >  2:2321550/1‑150 (MQ=255)
                                                                                                                                                   tGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCt    >  1:6697645/1‑150 (MQ=255)
                                                                                                                                                   tGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCt    <  1:1289276/150‑1 (MQ=255)
                                                                                                                                                     gCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGa  <  2:2954239/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAGCTGCAGCACCTGAAGGTGCCCACCAATTTCCTCAAGGCTCACCTTGAGCAAATGGAAGTGCGCAAGCCTGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGA  >  NC_002947/5187702‑5188000

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: