Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,270,455 C→T 15.2% G357G (GGC→GGT PP_4643 → xanthine/uracil permease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,270,4550CT15.2% 529.4 / 75.9 217G357G (GGC→GGTPP_4643xanthine/uracil permease family protein
Reads supporting (aligned to +/- strand):  ref base C (95/88);  new base T (16/17);  total (112/105)
Fisher's exact test for biased strand distribution p-value = 8.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 223 total aligned reads displayed.

AGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCA  >  NC_002947/5270306‑5270604
                                                                                                                                                     |                                                                                                                                                     
     gTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGATGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGgcgc                                                                                                                                                  >  2:17810/1‑150 (MQ=255)
            cAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATc                                                                                                                                           <  2:1753075/150‑1 (MQ=255)
            cAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATc                                                                                                                                           <  2:4891844/150‑1 (MQ=255)
            cAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATc                                                                                                                                           >  2:1134182/1‑150 (MQ=255)
              gAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGt                                                                                                                                         <  1:2746598/150‑1 (MQ=255)
                aCGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGc                                                                                                                                       <  2:2359592/150‑1 (MQ=255)
                     ggCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTACCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTc                                                                                                                                  <  1:3370524/150‑1 (MQ=255)
                      gCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCg                                                                                                                                 >  1:2849575/1‑150 (MQ=255)
                           gtggCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCAcc                                                                                                                            >  2:4645100/1‑150 (MQ=255)
                              gCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTg                                                                                                                         >  1:1550165/1‑150 (MQ=255)
                              gCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTg                                                                                                                         <  1:2302776/150‑1 (MQ=255)
                              gCCGTAACCGGCGTGAAAAGTCGTTATGTGGTGGCAACGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGtt                                                                                                                         >  2:2532160/1‑149 (MQ=255)
                               ccGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTgg                                                                                                                        <  2:243286/150‑1 (MQ=255)
                               ccGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTgg                                                                                                                        <  2:2973586/150‑1 (MQ=255)
                                 gTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTggcg                                                                                                                      <  2:1922105/150‑1 (MQ=255)
                                 gTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTggcg                                                                                                                      <  2:1204812/150‑1 (MQ=255)
                                     ccGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCa                                                                                                                  <  1:5020288/150‑1 (MQ=255)
                                      cGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAg                                                                                                                 <  2:2167485/150‑1 (MQ=255)
                                       ggggTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGc                                                                                                                >  2:2177345/1‑150 (MQ=255)
                                          gTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGt                                                                                                             <  2:1204710/150‑1 (MQ=255)
                                           tCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTa                                                                                                            >  2:3593689/1‑150 (MQ=255)
                                           tCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTa                                                                                                            >  1:4237566/1‑150 (MQ=255)
                                            cAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAt                                                                                                           <  2:2245417/150‑1 (MQ=255)
                                            cAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTAt                                                                                                           >  1:5398619/1‑150 (MQ=255)
                                             aaGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATc                                                                                                          <  2:1550165/150‑1 (MQ=255)
                                             aaGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATc                                                                                                          <  1:3999023/150‑1 (MQ=255)
                                              agagCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATcc                                                                                                         <  1:3344071/150‑1 (MQ=255)
                                              agagCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATcc                                                                                                         >  1:1977743/1‑150 (MQ=255)
                                              agagCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATcc                                                                                                         <  1:4815527/150‑1 (MQ=255)
                                                  ccGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCAc                                                                                                     <  1:4434460/150‑1 (MQ=255)
                                                   cGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCAcc                                                                                                    <  1:1624562/150‑1 (MQ=255)
                                                     cTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCt                                                                                                  <  1:133250/150‑1 (MQ=255)
                                                         gtggtgGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTcc                                                                                              >  1:4782939/1‑150 (MQ=255)
                                                               gcaaCGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTc                                                                                        <  2:3043287/147‑1 (MQ=255)
                                                                cTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTCTCGTGCTGTTCGGCCCCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCa                                                                                       >  2:4453578/1‑150 (MQ=255)
                                                                 tACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAg                                                                                      >  1:1698080/1‑150 (MQ=255)
                                                                 tACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAg                                                                                      >  2:2088506/1‑150 (MQ=255)
                                                                 tACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAg                                                                                      >  1:926926/1‑150 (MQ=255)
                                                                 tACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAg                                                                                      >  2:4279429/1‑150 (MQ=255)
                                                                 tACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAg                                                                                      >  2:3205599/1‑150 (MQ=255)
                                                                 tACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAg                                                                                      >  2:2538283/1‑150 (MQ=255)
                                                                 aaCGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTcga                                                                                      >  2:322020/2‑148 (MQ=255)
                                                                 aaCGGGCGGTCTTATTCTGATCACACTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTcga                                                                                      >  2:2345711/2‑148 (MQ=255)
                                                                     ggcggcGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTAc                                                                                  <  1:1938109/150‑1 (MQ=255)
                                                                     ggcggcGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTAc                                                                                  <  1:1528156/150‑1 (MQ=255)
                                                                           gTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGaac                                                                            <  1:1736176/150‑1 (MQ=255)
                                                                            tGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGaaca                                                                           <  2:2787202/150‑1 (MQ=255)
                                                                             gATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGaacaa                                                                          <  1:4402628/150‑1 (MQ=255)
                                                                                tCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGt                                                                       >  1:2199422/1‑150 (MQ=255)
                                                                                   ggtggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCaa                                                                    >  2:3316111/1‑150 (MQ=255)
                                                                                    gtggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAAc                                                                   <  2:1685298/150‑1 (MQ=255)
                                                                                      ggtACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCt                                                                 <  2:3423873/150‑1 (MQ=255)
                                                                                      cacaCTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCt                                                                 <  2:4802975/147‑1 (MQ=255)
                                                                                       gtACTCGGCTTGCTGCCAGTCATGAGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTg                                                                >  1:3059547/1‑150 (MQ=255)
                                                                                        tACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGa                                                               <  2:4409283/150‑1 (MQ=255)
                                                                                         aCTCGGTCTTCTGCCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGat                                                              >  1:344999/1‑150 (MQ=255)
                                                                                             ggtcttctgcCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGatcatt                                                          >  2:4034952/7‑149 (MQ=255)
                                                                                             ggCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGatcatc                                                          >  2:4670280/1‑150 (MQ=255)
                                                                                               tcttctgcCTGTAATGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCAttgt                                                        <  2:3857568/146‑1 (MQ=255)
                                                                                                   ctgcCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTcgcg                                                    >  2:669827/1‑150 (MQ=255)
                                                                                                   ctgcCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAACGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTcgcg                                                    <  2:1043963/150‑1 (MQ=255)
                                                                                                    tgcCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCgg                                                   >  2:4989341/1‑150 (MQ=255)
                                                                                                      ccAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGcc                                                 >  1:1251514/1‑150 (MQ=255)
                                                                                                      ccAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGcc                                                 >  1:1434578/1‑150 (MQ=255)
                                                                                                      ccAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGcc                                                 >  2:556453/1‑150 (MQ=255)
                                                                                                        aGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTc                                               >  1:432505/1‑150 (MQ=255)
                                                                                                            aTGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTg                                           <  2:4101509/150‑1 (MQ=255)
                                                                                                             tGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTgg                                          >  1:2103889/1‑150 (MQ=255)
                                                                                                             tGGGACGGGTGATTGCCTGTGTTCCGCCAGCAGTGCTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTggc                                                     <  1:1021270/139‑1 (MQ=255)
                                                                                                              gggCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTg                                           >  2:4549435/1‑148 (MQ=255)
                                                                                                              gggCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTg                                           <  1:4549435/148‑1 (MQ=255)
                                                                                                               ggCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGc                                        <  1:4279429/150‑1 (MQ=255)
                                                                                                               ggCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGAGCATCGTCGCGGCCTCGCTGGGc                                        >  2:3616627/1‑150 (MQ=255)
                                                                                                                gCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCt                                       >  1:5159594/1‑150 (MQ=255)
                                                                                                                 ccGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCgg                                                   <  1:5426771/137‑1 (MQ=255)
                                                                                                                  cGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTc                                     <  1:3316111/150‑1 (MQ=255)
                                                                                                                       gATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCAt                                <  1:802928/150‑1 (MQ=255)
                                                                                                                         ttGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGa                              >  1:624933/1‑150 (MQ=255)
                                                                                                                            ccgccgTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATcc                           >  1:3698887/1‑150 (MQ=255)
                                                                                                                             cgccgTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATccc                          >  2:4608486/1‑150 (MQ=255)
                                                                                                                                      tACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGc                 >  1:3114670/1‑150 (MQ=255)
                                                                                                                                         agcagtgcTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCAcc              <  2:76353/143‑1 (MQ=255)
                                                                                                                                          ccGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCg                                                                                                        <  1:4083068/59‑1 (MQ=255)
                                                                                                                                          ccGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCg                                                                                                        >  2:1331500/1‑59 (MQ=255)
                                                                                                                                          ccGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTcgc                                                     >  1:2479557/1‑110 (MQ=255)
                                                                                                                                          ccGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTcgc                                                     <  2:2479557/110‑1 (MQ=255)
                                                                                                                                          ccGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCg             >  1:1158400/1‑150 (MQ=255)
                                                                                                                                          ccGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCg             >  1:2238454/1‑150 (MQ=255)
                                                                                                                                             gtgcTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCAccagaa          >  1:618160/4‑146 (MQ=255)
                                                                                                                                             gTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGAcc          >  2:3641285/1‑150 (MQ=255)
                                                                                                                                              tgcTGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCATTGTGGCTACCTCCATCGGCTTCGGCATGATCCCGATTGCAGCAccagaat         >  2:140412/3‑145 (MQ=255)
                                                                                                                                              tACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCt         >  1:3603723/1‑150 (MQ=255)
                                                                                                                                              tACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCt         >  1:3002067/1‑150 (MQ=255)
                                                                                                                                               aCTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCtt        <  2:5398619/150‑1 (MQ=255)
                                                                                                                                                 tGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCAttgtgg                                                      <  1:4968971/102‑3 (MQ=255)
                                                                                                                                                 tGGGTGGTGCCGGTATCGTGCTGTTCGGAACGGTTGCCGCAAGCGGCATCCGCACACTTTCAAAGGTCGACTATCAGAACAACATCAATCTGATCAttgtgg                                                      >  2:4968971/1‑100 (MQ=255)
                                                                                                                                                 tGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCt      >  1:2018464/1‑150 (MQ=255)
                                                                                                                                                   ggcggcGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTac    >  2:3321868/1‑150 (MQ=255)
                                                                                                                                                    gcggcgCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTacc   <  2:294109/150‑1 (MQ=255)
                                                                                                                                                     cggcgCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTacca  <  1:4645100/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCGCGTTCGCCCAGAACGTCGGCCTGGTGGCCGTGACCGGGGTCAAGAGCCGCTACGTGGTGGCTACGGGCGGCGTGATTCTGGTGGTACTCGGCTTGCTGCCAGTCATGGGCCGGGTGATTGCCGCCGTGCCTACCCCGGTACTGGGCGGCGCCGGTATCGTGCTGTTCGGCACCGTGGCGGCCAGCGGTATCCGCACCCTGTCCAAGGTCAGCTACAAGAACAACGTCAACCTGATCATCGTCGCGGCCTCGCTGGGCTTCGGCATGATCCCGATTGCCGCACCGACCTTCTACCA  >  NC_002947/5270306‑5270604

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: