Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,557,910 T→G 11.0% intergenic (+576/+85) PP_4882 → / ← PP_5717 iron ABC transporter permease/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,557,9100TG11.0% 513.6 / 14.4 190intergenic (+576/+85)PP_4882/PP_5717iron ABC transporter permease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (47/116);  new base G (24/2);  total (72/118)
Fisher's exact test for biased strand distribution p-value = 5.99e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 197 total aligned reads displayed.

CGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACAGT  >  NC_002947/5557763‑5558060
                                                                                                                                                   |                                                                                                                                                       
cGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTc                                                                                                                                                       <  1:5199899/150‑1 (MQ=255)
cGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTc                                                                                                                                                       <  2:795374/150‑1 (MQ=255)
cGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTc                                                                                                                                                       <  1:1371379/150‑1 (MQ=255)
cGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTc                                                                                                                                                       <  1:3050343/150‑1 (MQ=255)
cGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTTCtta                                                                                                                                                       >  1:4442698/1‑149 (MQ=255)
  gTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGc                                                                                                                                                     <  2:5278645/150‑1 (MQ=255)
      tCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTGCGTCGCcaac                                                                                                                                                 >  1:3001907/1‑149 (MQ=255)
       cATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACATGGTCTTGAGGAATTCCGAACCAGAAGCCAAGTCCGGCAGGTCCGGGCTTTGCGAAGCAAAGc                                                                                                                                                <  2:3373280/150‑1 (MQ=255)
       cATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGc                                                                                                                                                <  2:4046728/150‑1 (MQ=255)
        aTTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGGTGTGGGTTGCAAAGcc                                                                                                                                               >  1:1582617/1‑150 (MQ=18)
        aTTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGcc                                                                                                                                               <  1:29421/150‑1 (MQ=255)
        aTTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGGTGCGTCGCAAAccc                                                                                                                                               >  2:985134/1‑150 (MQ=255)
         ttttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGGGTTTCTTAGCAAAGccc                                                                                                                                              >  1:804553/1‑150 (MQ=255)
         ttttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGGGTGGCAAAGccc                                                                                                                                              >  1:804541/1‑150 (MQ=255)
          tttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCg                                                                                                                                             <  2:2425691/150‑1 (MQ=255)
          tttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCg                                                                                                                                             <  2:3038223/150‑1 (MQ=255)
           ttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCCGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCgg                                                                                                                                            <  2:3957315/150‑1 (MQ=255)
               cgcTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAAc                                                                                                                                        <  2:3417470/150‑1 (MQ=255)
                gcTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGTCGCAAAGCCCGGAAACt                                                                                                                                       >  1:4624860/1‑150 (MQ=255)
                 ctctCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTg                                                                                                                                      <  1:397465/150‑1 (MQ=255)
                  tctcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgg                                                                                                                                     <  1:2141212/150‑1 (MQ=255)
                  tctcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGTAGCAAAACCCCGaaaaacg                                                                                                                                     >  1:1014586/1‑146 (MQ=17)
                  tctcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTGCTTAGCAAAGCCCCCAAACAgg                                                                                                                                     >  1:3405389/1‑150 (MQ=17)
                   ctcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                    <  2:2456921/150‑1 (MQ=255)
                    tcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgggg                                                                                                                                   >  2:5326333/1‑150 (MQ=255)
                    tcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgggg                                                                                                                                   >  2:3604519/1‑150 (MQ=255)
                    tcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTAGCAAAACCCCGAAAAAgggg                                                                                                                                   >  1:3506867/1‑150 (MQ=17)
                      aTAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTAGCAAAGCCCGGAAACAGGGGtt                                                                                                                                 >  1:3695853/1‑150 (MQ=255)
                       tAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTCTTAGCAAAACCCCGAAACTGGGGttt                                                                                                                                >  1:2628590/1‑150 (MQ=255)
                             aCTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCtg                                                                                                                          <  1:4583159/150‑1 (MQ=255)
                               tGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCtggt                                                                                                                        <  2:60771/150‑1 (MQ=255)
                               tGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCtggt                                                                                                                        >  2:4287532/1‑150 (MQ=255)
                                  gcATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgc                                                                                                                     >  2:5393941/1‑150 (MQ=255)
                                     tCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgctgc                                                                                                                  <  2:4517025/150‑1 (MQ=255)
                                        gTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTGGTTAGCAAATCCCGGAAACGGGGGTTTCAGCTggtgggtgggtg                                                                                                               >  2:3838507/1‑146 (MQ=17)
                                         ttGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGa                                                                                                              <  2:1743600/150‑1 (MQ=255)
                                           ggTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAgt                                                                                                            <  1:106090/150‑1 (MQ=255)
                                                 tCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAg                                                                                                      <  1:3838507/150‑1 (MQ=255)
                                                        tACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACaa                                                                                               <  1:5178765/150‑1 (MQ=255)
                                                         aCAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAAt                                                                                              <  1:4276892/150‑1 (MQ=255)
                                                           aGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATgg                                                                                            <  2:3549142/150‑1 (MQ=255)
                                                              cGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGtt                                                                                         <  1:5116512/150‑1 (MQ=255)
                                                                      cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGGTTTTCTTCGCAAAGCCCGGAAAATGGGGTTTGAGCTGGTGGGTGGGTGTGTGTTTTGAAAACAAAGGGTTGCATTCTg                                                                                >  2:1538698/1‑150 (MQ=17)
                                                                        agaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGtt                                                                              <  2:3474450/150‑1 (MQ=255)
                                                                        agaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGtt                                                                              <  1:5165075/150‑1 (MQ=255)
                                                                         gaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGttt                                                                             <  1:3825965/150‑1 (MQ=255)
                                                                         gaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGttt                                                                             <  2:1007023/150‑1 (MQ=255)
                                                                              aGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGGGTTTGGTGGCAAAACCCCGAAACAGGGGTTTAAGCGGGGGGCTGCGTTTGTTTTTAGAAAAAAATTGGGTTCATTGTGtttttccc                                                                        >  2:485382/1‑145 (MQ=12)
                                                                               ggCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAg                                                                       <  2:1130267/150‑1 (MQ=255)
                                                                                gCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAgg                                                                      <  1:4874408/150‑1 (MQ=255)
                                                                                 cTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAggg                                                                     <  1:3635434/150‑1 (MQ=255)
                                                                                 cTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAggg                                                                     <  2:2061374/150‑1 (MQ=255)
                                                                                 cTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAggg                                                                     <  2:2739297/150‑1 (MQ=255)
                                                                                  ttATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTATAAAACAAAGGGGGGCATTCTGTTTGACCAGGGc                                                                    >  2:2826088/1‑150 (MQ=255)
                                                                                  ttATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGc                                                                    <  1:4342876/150‑1 (MQ=255)
                                                                                   tATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCGGCTTTGGTGTTTATAAAAACATGGGTTGGATTGTGTTTGACCAGGGCt                                                                   >  1:5052674/1‑150 (MQ=17)
                                                                                     tGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGt                                                                               <  1:1297851/136‑1 (MQ=255)
                                                                                      gCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTAGCTTCGCAAATCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTAc                                                                >  2:3274430/1‑150 (MQ=255)
                                                                                        tGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTGCGTCGCAAAGCCCGGAAACATGGGTTTGAGCGGGTGGCGGGTTGGGTGTTTTGTAAACAAAGGGGTTGATTATTTTTGACCAGGGCGTAACt                                                              >  2:1039174/1‑150 (MQ=17)
                                                                                           aCAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCaa                                                           <  2:2893836/150‑1 (MQ=255)
                                                                                               ggTCTTGAGGAATTCCCAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCCGAAAAGGGGGTTTCAGCATGTGGCTGCGTGAGTGTTTAGAAAACAAAGCGTTGCATTCTGTTTGACCAGGGCGAACCACAAGTaa                                                       >  2:4780518/1‑149 (MQ=17)
                                                                                                    tGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTAt                                                  <  2:4583275/150‑1 (MQ=255)
                                                                                                       ggAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTGCTTCGCAAAGCCCCGAAACTGGGGTTTCAGCTGGTGGCTGCGTGAGTGTTTATAAAAAAATGGGTTGGATTCTGTTTTAACAGGGCGTACAAAAAGTACGCTATCCa                                               >  2:5213347/1‑150 (MQ=17)
                                                                                                         aaTTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCaaa                                             <  2:5135890/150‑1 (MQ=255)
                                                                                                         aaTTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCaaa                                             <  1:3297083/150‑1 (MQ=255)
                                                                                                         aaTTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCaaa                                             <  1:934561/150‑1 (MQ=255)
                                                                                                          aTTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAg                                            <  1:11217/150‑1 (MQ=255)
                                                                                                              cGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCGGCTTGAGTGTTTTGAAAACCATGGGTTGGATTCTGTTTGACCAGGGCGTACCACCAGTACTCGATGCCAAGaaaa                                        >  1:4181673/1‑150 (MQ=17)
                                                                                                                aaCCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTGGTAGCAAAGCCCCGCAACGGGGGTTTCAGCGGGTGGGTGGGTGGGTGGTTGTAAAAAACAGGGTTTCATTGTGGTTTTACCGGGCTTTCCTCCACTGCTAGATCCAAAGaaaaaa                                      >  2:210580/1‑150 (MQ=9)
                                                                                                                aaCCAGAAGCCAAGACCGGCAGGTCCGGGCGTTTCGTGGCAAACCCCCGACACGGGGGTGTGATCGTGGGGGGGGGTGAGTGTTTAGAAAAACATGGGTTGGCTTCGGGTTTGCCAGGGCGTACCTCAAGTACTCGATACCAAGaaaaaa                                      >  2:4713487/1‑150 (MQ=9)
                                                                                                                  ccAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaaa                                    <  1:5189847/150‑1 (MQ=255)
                                                                                                                  ccAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaaa                                    <  1:1918038/150‑1 (MQ=255)
                                                                                                                   cAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAAATGGGGTTTCAGCTGGTGGCGGCTTGAGTGTTTAGAAAACAAAGGGTTGCAGTCTGTTTGACCAGGGCTTAACTCAAGTACTCTATACAAAGAAAAAAAAg                                   >  2:1249921/1‑150 (MQ=21)
                                                                                                                     gAAGCCAAGACCGGCAGGTCCGGGCTTTTCATAGCCAAGCCCCGAAACTGGGGTTTAAGCTGGTGGCTGCTTGAGTGTTTAGAAAAAAAAGGGTTGCATTCTGTTTGTCCCGGGCGTAAATCAAGTAAGCAATACAAATACAAAAAAGaa                                 >  2:3456174/1‑150 (MQ=17)
                                                                                                                     gAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACACTGGGTTGCATTCTGTTTGTCCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGaa                                 <  2:3819949/150‑1 (MQ=255)
                                                                                                                     gAAGCCAAGACCGGCAGGTCCGGGCGTTTAGTGGCAAAGCCCCGAAAAGGGGGTTTGAGCGGGTGGGTGCTTGGTTGTTTATAAAACAAAGGGTTGCCTTCTGTTTGTCCAGGGGGTTCCTCAAGTACTAAATCCCAAGAAAAAAAAGaa                                 >  2:1995741/1‑150 (MQ=16)
                                                                                                                        gCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACATTGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAATAAAGAACaa                              >  2:2303601/1‑150 (MQ=255)
                                                                                                                           aaGACCGGCAGGTCCGGGCTTTTCTTGGCAAATCCCGGAAACTGGGGTTTAAGCGGGTGGCTGGGGGTGTGTTTAGAAAACAAGGGGGTGCATTGTGTTTGACCAGGGCTTTACACAAGTACGCGATACCAAGAACAAAAAGaacaaaca                           >  2:4239603/1‑150 (MQ=12)
                                                                                                                            aGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAg                          <  2:5616045/150‑1 (MQ=255)
                                                                                                                              aCCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTg                        <  2:804553/150‑1 (MQ=255)
                                                                                                                              aCCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTg                        <  2:804541/150‑1 (MQ=255)
                                                                                                                                  gCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAAt                    <  2:4913672/150‑1 (MQ=255)
                                                                                                                                  gCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAAt                    <  1:5021346/150‑1 (MQ=255)
                                                                                                                                       tCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTAATCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCa               >  2:2533963/1‑150 (MQ=255)
                                                                                                                                        ccGGGCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCaa              >  1:2849237/1‑150 (MQ=255)
                                                                                                                                             gtttGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTTCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAAttt         <  2:3647294/149‑1 (MQ=255)
                                                                                                                                              tttGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACACTCAGTGCAATCAGCAAAATTTg        <  2:753127/150‑1 (MQ=255)
                                                                                                                                              tttGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTg        <  2:4624860/150‑1 (MQ=255)
                                                                                                                                               ttGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGt       <  2:3752540/150‑1 (MQ=255)
                                                                                                                                               ttGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGt       <  1:2374829/150‑1 (MQ=255)
                                                                                                                                                tGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTa      <  1:5393941/150‑1 (MQ=255)
                                                                                                                                                tGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGAATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTa      <  2:3918712/150‑1 (MQ=255)
                                                                                                                                                 gCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTAc     >  2:3823298/1‑150 (MQ=255)
                                                                                                                                                  cTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACa    <  1:4287532/150‑1 (MQ=255)
                                                                                                                                                  cTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACa    >  2:1909958/1‑150 (MQ=255)
                                                                                                                                                   ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTg                        <  2:1112912/129‑1 (MQ=255)
                                                                                                                                                   ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACAg   >  2:332161/1‑150 (MQ=255)
                                                                                                                                                   ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACAg   >  2:170926/1‑150 (MQ=255)
                                                                                                                                                   ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACAg   >  2:1307516/1‑150 (MQ=255)
                                                                                                                                                   ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCATAAGCAAAATTTGTaaatt  >  1:785392/1‑146 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
CGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACAGT  >  NC_002947/5557763‑5558060

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: