Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,228,621 T→A 6.2% intergenic (+13/+17) PP_4605 → / ← PP_4606 AraC family transcriptional regulator/ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,228,6210TA 6.2% 639.1 / 22.1 194intergenic (+13/+17)PP_4605/PP_4606AraC family transcriptional regulator/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (83/99);  new base A (5/7);  total (88/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

Only 100 of 197 total aligned reads displayed.

GTTGCGGCTGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTTGTGT  >  NC_002947/5228473‑5228770
                                                                                                                                                    |                                                                                                                                                     
     ggCTGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCgg                                                                                                                                                 >  1:5316431/1‑150 (MQ=255)
     ggCTGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCgg                                                                                                                                                 >  2:1404238/1‑150 (MQ=255)
        tGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGccc                                                                                                                                              >  2:3981810/1‑150 (MQ=255)
        tGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCAGAGCGGccc                                                                                                                                              <  2:3655281/150‑1 (MQ=255)
        tGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTGTCATCTGAAACCGGGGCCGCAGAGCGGccc                                                                                                                                              <  2:3218995/150‑1 (MQ=255)
         gggAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCAGAGCGGcccc                                                                                                                                             <  2:1974587/150‑1 (MQ=255)
         gggAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCAGAGCGGcccc                                                                                                                                             <  2:4027735/150‑1 (MQ=255)
          ggAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCa                                                                                                                                            >  1:2667150/1‑150 (MQ=255)
             gTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCAttt                                                                                                                                         <  1:5402804/150‑1 (MQ=255)
              tAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTg                                                                                                                                        >  2:133610/1‑150 (MQ=255)
                aTGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGcc                                                                                                                                      <  1:637119/150‑1 (MQ=255)
                aTGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGcc                                                                                                                                      >  1:715845/1‑150 (MQ=255)
                aTGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGcc                                                                                                                                      >  2:4735222/1‑150 (MQ=255)
                aTGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCAGTTGcc                                                                                                                                      >  1:3247047/1‑150 (MQ=255)
                aTGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTAGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCAGAGCGGccccggtttca                                                                                                                                      <  2:4190516/150‑8 (MQ=255)
                  gCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTa                                                                                                                                    >  1:222639/1‑150 (MQ=255)
                   ccGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTaa                                                                                                                                   <  2:3345344/150‑1 (MQ=255)
                   ccGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTaa                                                                                                                                   <  1:3493105/150‑1 (MQ=255)
                          gagagaTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGa                                                                                                                            >  2:3443039/1‑150 (MQ=255)
                           agagaTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGaa                                                                                                                           <  2:211370/150‑1 (MQ=255)
                            gagaTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAg                                                                                                                          <  1:361728/150‑1 (MQ=255)
                             agaTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAgg                                                                                                                         >  1:5084780/1‑150 (MQ=255)
                              gaTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATGTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGt                                                                                                                        <  1:2931635/150‑1 (MQ=255)
                               aTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTg                                                                                                                       <  1:4462056/150‑1 (MQ=255)
                               aTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTGTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTg                                                                                                                       <  2:3128694/150‑1 (MQ=255)
                                 tGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCt                                                                                                                     <  1:2980794/150‑1 (MQ=255)
                                  gCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGTGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTc                                                                                                                    <  2:1491606/150‑1 (MQ=255)
                                  gCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTc                                                                                                                    <  2:4961861/150‑1 (MQ=255)
                                  gCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTc                                                                                                                    <  1:2416207/150‑1 (MQ=255)
                                  gCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTc                                                                                                                    <  1:1969684/150‑1 (MQ=255)
                                  gCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTc                                                                                                                    <  1:1413956/150‑1 (MQ=255)
                                    cTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGa                                                                                                                  >  1:2187046/1‑150 (MQ=255)
                                      tATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACa                                                                                                                >  2:5450510/1‑150 (MQ=255)
                                        ttGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACaaa                                                                                                              <  1:4931551/150‑1 (MQ=255)
                                         tgtgCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTGCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTGTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAAc                                                                                                             <  1:3018357/150‑1 (MQ=255)
                                            gCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAATACTCGAAGGTGCTCGACAAACTGa                                                                                                          >  1:4542167/1‑150 (MQ=255)
                                                gTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTg                                                                                                      <  1:953590/150‑1 (MQ=255)
                                                gTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCATGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCACATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTg                                                                                                      <  2:119478/150‑1 (MQ=255)
                                                   tgccagtgccagCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCcg                                                                                                   <  2:2351288/150‑1 (MQ=255)
                                                    gccagtgccagCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCcgc                                                                                                  <  1:1239549/150‑1 (MQ=255)
                                                     ccagtgccagCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGTGGCCCCATTTGCCTAAAACTCGAGGTGCTCGACAAACTGACCTGc                                                                                                     >  2:3232578/1‑145 (MQ=255)
                                                     ccagtgccagCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCcgcg                                                                                                 >  1:5066303/1‑150 (MQ=255)
                                                      cagtgccagCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCAttt                                                                                                                                         <  2:3711653/109‑1 (MQ=255)
                                                      cagtgccagCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCAttt                                                                                                                                         >  1:3711653/1‑109 (MQ=255)
                                                      cagtgccagCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCgcgc                                                                                                >  1:4959667/1‑150 (MQ=255)
                                                      cagtgccagCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCgcgc                                                                                                <  2:746150/150‑1 (MQ=255)
                                                          gccagccaTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcg                                                                                            <  1:1530107/150‑1 (MQ=255)
                                                            cagccaTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTc                                                                                                                    >  2:5392532/1‑124 (MQ=255)
                                                             agccaTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgcca                                                                                         >  1:1003138/1‑150 (MQ=255)
                                                                caTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaatc                                                                                      <  2:2736519/150‑1 (MQ=255)
                                                                         aaCCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAg                                                                             >  1:3786284/1‑150 (MQ=255)
                                                                           ccGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTa                                                                           >  2:4681510/1‑150 (MQ=255)
                                                                              gTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTg                                                                        >  1:1599342/1‑150 (MQ=255)
                                                                                  cGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTc                                                                    <  2:1298318/150‑1 (MQ=255)
                                                                                  cGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACATGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTc                                                                    <  2:1370000/150‑1 (MQ=255)
                                                                                    ccAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCAc                                                                  <  1:1799947/150‑1 (MQ=255)
                                                                                    ccAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCAc                                                                  <  2:2660402/150‑1 (MQ=255)
                                                                                          ccTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCACTCGCCACAAAGTACTGGTTCACCGCCGa                                                            >  1:1016519/1‑150 (MQ=255)
                                                                                          ccTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGa                                                            <  1:4487088/150‑1 (MQ=255)
                                                                                          ccTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGa                                                            <  1:2779529/150‑1 (MQ=255)
                                                                                          ccTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCAACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGa                                                            <  1:399339/150‑1 (MQ=255)
                                                                                            tGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGCCAAACTGACCTGCCGCGCATCGCCACTCGCCACACAGTACGGGTGCACCGCCGAGc                                                          >  2:4072029/1‑150 (MQ=255)
                                                                                            tGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGc                                                          >  1:3326614/1‑150 (MQ=255)
                                                                                             gggggCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCt                                                         <  1:4300632/150‑1 (MQ=255)
                                                                                              ggggCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTg                                                        <  2:3856790/150‑1 (MQ=255)
                                                                                                 gCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTaa                                                     <  2:446897/150‑1 (MQ=255)
                                                                                                 gCCACACCTGGCGAATACCGCAACGCCTATTTCAACTGAAACCGGGGCCGCTCTGCGGCCCCAATTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTCCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGACCTGTaa                                                     <  2:872786/150‑1 (MQ=255)
                                                                                                      acCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaatcgc                                                                                    >  1:3224691/1‑114 (MQ=255)
                                                                                                       ccTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaatcgc                                                                                    >  2:4607550/1‑113 (MQ=255)
                                                                                                          ggCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTg                                            >  2:1242010/1‑150 (MQ=255)
                                                                                                           gCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaatcgc                                                                                    >  2:2957835/1‑109 (MQ=255)
                                                                                                            cGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTc                                          >  1:406092/1‑150 (MQ=255)
                                                                                                             naaTACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCg                                         >  1:2913061/2‑150 (MQ=255)
                                                                                                             gAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCg                                         >  1:399727/1‑150 (MQ=255)
                                                                                                                 aCCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACa                                     >  2:1426127/1‑150 (MQ=255)
                                                                                                                   cGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAgg                                   >  1:2109970/1‑150 (MQ=255)
                                                                                                                   cGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAgg                                   <  1:2674383/150‑1 (MQ=255)
                                                                                                                      aaGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGttc                                <  1:685660/150‑1 (MQ=255)
                                                                                                                         gacTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCcac                                                                                 <  2:3881932/96‑1 (MQ=255)
                                                                                                                         gacTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCcac                                                                                 >  1:3881932/3‑98 (MQ=255)
                                                                                                                         gCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGttcttc                             <  2:2679288/150‑1 (MQ=255)
                                                                                                                           cTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCAc                           <  2:222639/150‑1 (MQ=255)
                                                                                                                                tCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCa                      <  2:5316431/150‑1 (MQ=255)
                                                                                                                                  aTCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGt                    >  1:5141599/1‑150 (MQ=255)
                                                                                                                                  aTCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGt                    <  1:120591/150‑1 (MQ=255)
                                                                                                                                   tCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGtt                   >  2:1770982/1‑150 (MQ=255)
                                                                                                                                   tCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAAACGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGtt                   >  1:25513/1‑150 (MQ=255)
                                                                                                                                    cTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCt                                                         >  2:1164649/1‑111 (MQ=255)
                                                                                                                                     tGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACa                                     <  2:2926737/130‑1 (MQ=255)
                                                                                                                                     tGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAAATGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGc                 >  1:3937933/1‑150 (MQ=255)
                                                                                                                                      gcaaatggggCCGCAGAGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCa                >  2:404268/7‑150 (MQ=255)
                                                                                                                                      gAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCa                <  2:5142969/150‑1 (MQ=255)
                                                                                                                                      gAAACCGGGGCCGCTCAGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCa                <  1:3932992/150‑1 (MQ=255)
                                                                                                                                        aaatggggCCGCAGAGCGGCCCCGTTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTg                                                                        <  1:2021914/88‑1 (MQ=38)
                                                                                                                                          ccGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGgcg            >  2:1984275/1‑150 (MQ=255)
                                                                                                                                             gggCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCAc         <  2:5059153/150‑1 (MQ=255)
                                                                                                                                                  gCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACaaa                                                                              <  2:4263802/76‑1 (MQ=255)
                                                                                                                                                  gCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACaaa                                                                              >  1:4263802/1‑76 (MQ=255)
                                                                                                                                                  gCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTtgt    <  1:1447293/150‑1 (MQ=255)
                                                                                                                                                    tctGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTTgtgt  >  1:4107017/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTGCGGCTGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTTGTGT  >  NC_002947/5228473‑5228770

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: