Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 85,081 T→G 8.5% intergenic (+59/+71) aroE → / ← PP_0075 shikimate 5‑dehydrogenase/choline sulfate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294785,0810TG 8.5% 594.7 / 39.6 223intergenic (+59/+71)aroE/PP_0075shikimate 5‑dehydrogenase/choline sulfate transporter
Reads supporting (aligned to +/- strand):  ref base T (85/111);  new base G (10/9);  total (96/127)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-02

Only 100 of 236 total aligned reads displayed.

GCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCG  >  NC_002947/84933‑85229
                                                                                                                                                     |                                                                                                                                                    
gcagCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtc                                                                                                                                                     >  2:4376518/1‑150 (MQ=255)
  agCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtcgt                                                                                                                                                   >  2:4839876/1‑150 (MQ=255)
  agCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTAcgtcgt                                                                                                                                                   >  1:5240521/1‑150 (MQ=255)
  agCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTACCCCCAAGTAcgtcgt                                                                                                                                                   >  1:3449792/1‑150 (MQ=255)
  agCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAACGCCGTTACCACAAAATAactcgt                                                                                                                                                   >  1:936818/1‑150 (MQ=255)
   gCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTa                                                                                                                                                  <  1:1771367/150‑1 (MQ=255)
      gAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGc                                                                                                                                               >  2:1281170/1‑150 (MQ=255)
       aaGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGAGTAGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGcc                                                                                                                                              >  2:4121590/1‑150 (MQ=255)
            ttcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATcc                                                                                                                                         >  2:2922267/1‑150 (MQ=255)
            ttcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTacgacgtac                                                                                                                                          >  2:2285662/1‑147 (MQ=37)
             tcttTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCct                                                                                                                                        <  2:4957113/150‑1 (MQ=255)
               tttATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCctct                                                                                                                                      <  2:3499633/150‑1 (MQ=255)
                 tATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                    <  2:1743425/150‑1 (MQ=255)
                 tATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACCCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                    <  2:1378391/150‑1 (MQ=255)
                 gaTCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgt                                                                                                                                    <  2:3662046/149‑1 (MQ=255)
                  agcgggcgtggcgaGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCtgtg                                                                                                                                   <  2:5500879/146‑1 (MQ=255)
                     tggcgtggcgTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGa                                                                                                                                <  1:426430/150‑1 (MQ=255)
                       gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGaac                                                                                                                              <  1:2285662/150‑3 (MQ=25)
                       gcgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGaac                                                                                                                              <  1:84540/150‑3 (MQ=21)
                        cgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACgg                                                                                                                            <  2:273559/150‑1 (MQ=17)
                        cgtggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTTGCAGCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACgg                                                                                                                            <  1:3663128/150‑1 (MQ=17)
                          tggcgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCCATCGCGACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGc                                                                                                                           <  1:3101926/150‑1 (MQ=21)
                             cgtgAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCtt                                                                                                                        <  2:1436106/150‑1 (MQ=255)
                                gAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCg                                                                                                                     <  1:456855/150‑1 (MQ=255)
                                 aGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGt                                                                                                                    <  1:2880433/150‑1 (MQ=255)
                                  ggCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTc                                                                                                                   <  1:317653/150‑1 (MQ=255)
                                    ccGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCTTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTcgc                                                                                                                 <  1:2931977/150‑1 (MQ=255)
                                    ccGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTcgc                                                                                                                 <  1:5023263/150‑1 (MQ=255)
                                       gATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGAt                                                                                                              <  1:373130/150‑1 (MQ=255)
                                          aCGGCGCCGGTGTTGGCTGACCTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATggg                                                                                                           <  2:1514575/150‑1 (MQ=255)
                                            ggCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCt                                                                                                         >  2:1894553/1‑150 (MQ=255)
                                             gcgcCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTg                                                                                                        >  2:3274160/1‑150 (MQ=255)
                                                 cGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCaaa                                                                                                    <  2:1789715/150‑1 (MQ=255)
                                                   gtgtTGGCTGACCTGCGCCGGAACCGCGCGCGGGGCTGCGCTGGCCGGGGCCGCTTGGCGGCCCTACCGCGACGCAAGGCCGTTCCCACAATGTACGTCGTACGCCGAGCCTCTGTGGGCGCGGCCTTGCGTCGCGATGGGCTGCAAAGc                                                                                                  <  2:8355/150‑1 (MQ=37)
                                                   gtgtGGGCTTAACTGCGCCGGCAACTCGCGCGGGGCTGAGCTGGCGGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTCCCCAAAAAGTACGTCTTACGCCGACCCTCGGTGGGAGCGGCCTTGCGCCGCGATGGGCTGCAAAGc                                                                                                  <  1:377190/150‑1 (MQ=255)
                                                    tgtTGGCTTAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCGGGGGCCGATTGGCGGCCCTATCGCGACGCACGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCGGGGGGAGCGGCCTTGCGACGCGATGGGCTGCAAAGCa                                                                                                 <  1:848646/150‑1 (MQ=255)
                                                           tGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGGCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAg                                                                                          >  2:2999862/1‑150 (MQ=255)
                                                            gAACTGCGCCGGCAACTCGCTCGGGGCTGCGCGGGCGGGGGCCGCTTGGCGGCCCTACCGCGACGCAAGGCCGTGCCCACAACGTACGTCGTCCGCCGATCCTCTGTGGGCGCGGCCTTGCGCCGCGATGGGCTGCAAAGCAGACCCAgg                                                                                         <  1:3274160/150‑2 (MQ=255)
                                                             aaCTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGGTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCGTGCGTCGCGATGGGCTGCAAAGCAGACCCAGaa                                                                                        >  1:2172389/1‑150 (MQ=255)
                                                             aaCTGCGCCGGCAACGCGCTCGGGGCTGAGCGGGCCGGGGCCGCTGGGCGGCCCTATCGCGCCGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCACTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGaa                                                                                        <  2:2998748/150‑1 (MQ=38)
                                                               cTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGACCCAGAATc                                                                                      <  1:4525547/150‑1 (MQ=255)
                                                                 gcgcCGGCAACTCGCTGGGGGCTGTGGTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAATGTACGTCGTACGCCGATCCTGTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATccc                                                                                    <  1:4547785/150‑1 (MQ=38)
                                                                 gcgcCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATccc                                                                                    <  1:3677112/150‑1 (MQ=255)
                                                                   gcCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCaa                                                                                  <  2:3716251/150‑1 (MQ=255)
                                                                    ccGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCACGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAAc                                                                                 <  2:923933/150‑1 (MQ=255)
                                                                     cGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAAcc                                                                                >  2:650612/1‑150 (MQ=255)
                                                                     cGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCGTGCGTCGCGATGGGCTGCACAGCAGACCCAGAAACCCAAcc                                                                                >  2:3063838/1‑150 (MQ=255)
                                                                           ctcgctcgGGGCTGAGCTGGCGGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTc                                                                          <  2:5203492/150‑1 (MQ=255)
                                                                            tcgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTcc                                                                         >  2:3230548/1‑150 (MQ=255)
                                                                             cgctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCt                                                                        <  2:3129407/150‑1 (MQ=255)
                                                                               ctcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCATGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCa                                                                      <  2:1962651/150‑1 (MQ=255)
                                                                               ctccGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCa                                                                      >  2:5462715/1‑150 (MQ=255)
                                                                                tcgGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCaa                                                                     >  2:3690310/1‑150 (MQ=255)
                                                                                    ggCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCg                                                                 <  2:1026048/150‑1 (MQ=255)
                                                                                     gCTGAGCTGGCGGGGGCCGGTTGGCGGCCCTATCGCGACGCAAGGCGGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGGGGCCTTGCGTCGCGATGGGCTGCATAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCgg                                                                <  2:214014/150‑1 (MQ=39)
                                                                                     gCTGAGCGGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCgg                                                                <  2:3934969/150‑1 (MQ=255)
                                                                                          gCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCgg                                                           <  2:453242/150‑1 (MQ=255)
                                                                                           gtGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCACGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCggg                                                          <  1:4007980/149‑1 (MQ=255)
                                                                                           cTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCggg                                                          >  2:5423871/1‑150 (MQ=255)
                                                                                            tGGCGGGGGCCGCTTGGCGGCCCTATCGCGACGAAAGGCCGTTCCCACAAAGTACGCCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGc                                                         <  1:1768374/150‑1 (MQ=255)
                                                                                             ggCGGGGGCCGTTTGGCGGCCCTATCGCGCCGAATGGCCGTTCCCAAAACTGACTCCGTCCGCCTATCCTGTGTGGTAGGGGCCTTGCGCCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca                                                        <  1:101912/150‑1 (MQ=255)
                                                                                             ggCGGGGGCCGCTTGGCGGCCCTATCGGCACTCACGGCGGTCCCCCAAACGTACGCCGAACGCCGACCCGCTGTGGGCGGGGCCTTGCGCCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca                                                        <  1:675711/150‑1 (MQ=255)
                                                                                             cgatagggCCGCCAAGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGca                                                        >  2:2361251/6‑150 (MQ=255)
                                                                                                      cGCTTGGCGGCCCTATCGCGACGCAAGGCCATTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAc                                               <  1:4931942/150‑1 (MQ=255)
                                                                                                      cGCTTGGCGGCCCTATCGCGACGCAAGGCCATTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAc                                               <  1:2100527/150‑1 (MQ=255)
                                                                                                       gCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCAcc                                              >  2:3553350/1‑150 (MQ=255)
                                                                                                            gcggcACTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGACCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCACATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTc                                         <  2:387583/150‑1 (MQ=255)
                                                                                                              ggcCCTCTGGGGAGGCAGGGCGGTGCCCACAAGTTTCGTGTTCCGCCGAGCGTGTGGGGGCGGGGCCTGGCGACGGGATGGGCTGCAAAGCAGTCCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCa                                       <  1:4798239/150‑1 (MQ=255)
                                                                                                              ggcCCTATGGCGCCGCAAGGCGGTTCCCACAATGTAGGTCGTACGCCGATCCTGTGTGGGAGCGGCCTTGGGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCa                                       <  1:4376518/150‑1 (MQ=38)
                                                                                                              ggcCCTATCGCGACGCAAGGCCTTTCCCAAAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCa                                       <  1:1311020/150‑1 (MQ=255)
                                                                                                              ggcCCTATCGCGACGCAAGGCCGTTCCCACAAATTACGTCGTACGACGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCa                                       <  2:2087535/150‑1 (MQ=255)
                                                                                                               gcCCTTTCGCCACGCATGGCCGTTCCCAAAACGTACGTGGTACGCCGATCCTCTGTGGGTGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAg                                      <  1:3523966/150‑1 (MQ=39)
                                                                                                                ccatagcgcgACGCAAGGCCGTTCCAACAACGTCGGCCGTCCGCCGAGCATCGGTGGGTGGGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGc                                     <  1:3323889/144‑1 (MQ=255)
                                                                                                                  cTATCGCGACGAATGGCGGTTCCCAAAAAGTACGTCTTACGCGTATCCTCTGTGGGTGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCtt                                   <  1:3602571/150‑1 (MQ=38)
                                                                                                                   taccgctACGCAAGGCCGTTCCCACAACCTACGTCGTACGCCTAACCTCGGGGGGTGGGGCCTTGCGCCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTc                                  <  1:125787/147‑1 (MQ=255)
                                                                                                                   tATCGCGACGCACGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTACAGCTTc                                  <  2:5269364/150‑1 (MQ=255)
                                                                                                                    aTCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATACCAACCACATTCCACAAACCGGATCGGGCACAACGCCACCCCATCCAGCATCt                                 >  1:2661264/1‑150 (MQ=39)
                                                                                                                      cgcgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCACAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGc                               >  2:33377/1‑150 (MQ=255)
                                                                                                                       gcgACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCa                              >  2:676361/1‑150 (MQ=255)
                                                                                                                         gACGCAAGGCCGCTCCCACAGAGGATCGGCGTACGACGTACTTTGTGGGAACGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAAc                            >  2:1704739/1‑150 (MQ=21)
                                                                                                                          aCGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAAct                           >  1:3465979/1‑150 (MQ=255)
                                                                                                                            gCACGGCCGTGCCCACAACTCCCGTCGTACGCCCACCCTCTGTGGGAGCGGCCTCCCGACGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAACCCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcc                         <  1:2215212/150‑1 (MQ=255)
                                                                                                                            gCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAActcc                         >  2:4234043/1‑150 (MQ=255)
                                                                                                                                ggCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGa                     <  2:1816866/150‑1 (MQ=255)
                                                                                                                                ggCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGAATTGCGTCGCGATTGGCTGCAAAGAAGACCCAGAATACCAACCTCATTCCTCAAACAGGATCGGGCACAACTCCGCACCTTCAAGCTTCTGCAAATCCACGa                     >  2:1861122/1‑150 (MQ=255)
                                                                                                                                 gCCGTTCCCACAAAGTACGTCGTACGCCGATCCGCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCAGTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGag                    >  2:3600431/1‑149 (MQ=255)
                                                                                                                                  ccGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATc                   <  2:3352653/150‑1 (MQ=255)
                                                                                                                                  ccGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATc                   >  1:809378/1‑150 (MQ=255)
                                                                                                                                  ccGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATc                   >  1:406986/1‑150 (MQ=255)
                                                                                                                                     ttCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCAcc                <  1:1827381/150‑1 (MQ=255)
                                                                                                                                        ccACAAAGTACTTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCAGTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGc             <  2:4838594/150‑1 (MQ=255)
                                                                                                                                        ccACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCTACTCCTCGATCACCTGc             >  2:4659356/1‑150 (MQ=255)
                                                                                                                                          acaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGTAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCACCTCCTCGATCACCTGCgg           >  2:4495600/1‑150 (MQ=255)
                                                                                                                                          acaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCgg           >  1:3969095/1‑150 (MQ=255)
                                                                                                                                          acaAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCGCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCgg           >  1:5537087/1‑150 (MQ=255)
                                                                                                                                           caAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGc          <  2:2186462/150‑1 (MQ=255)
                                                                                                                                             aaGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCt        <  2:5620987/150‑1 (MQ=255)
                                                                                                                                                   cgtcgtACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATACCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGccc   >  1:662858/1‑150 (MQ=255)
                                                                                                                                                    gtcgtACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGCATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCg  >  1:4238210/1‑150 (MQ=255)
                                                                                                                                                    gtcgtACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCg  >  2:2505916/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAGCTGAAGCCTTCTTTATCTGGCGTGGCGTGAGGCCGGATACGGCGCCGGTGTTGGCTGAACTGCGCCGGCAACTCGCTCGGGGCTGAGCTGGCCGGGGCCGCTTGGCGGCCCTATCGCGACGCAAGGCCGTTCCCACAAAGTACGTCGTACGCCGATCCTCTGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCAAAGCAGACCCAGAATCCCAACCTCATTCCTCAAACCGGATCGGGCACAACTCCACCCCTTCCAGCTTCTGCAACTCCTCGATCACCTGCGGCCTGGCCCG  >  NC_002947/84933‑85229

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: