Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,557,912 C→A 15.0% intergenic (+578/+83) 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,9120CA15.0% 591.5 / 58.2 244intergenic (+578/+83)PP_4882/PP_5717iron ABC transporter permease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (86/112);  new base A (31/9);  total (123/121)
Fisher's exact test for biased strand distribution p-value = 1.00e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

CGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACAGTC  >  NC_002947/5557763‑5558061
                                                                                                                                                     |                                                                                                                                                     
cGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTc                                                                                                                                                       <  1:3382400/150‑1 (MQ=255)
  gTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGc                                                                                                                                                     <  2:4186583/150‑1 (MQ=255)
   tACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCa                                                                                                                                                    >  2:4399150/1‑150 (MQ=255)
     cTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCaaa                                                                                                                                                  <  1:1664362/150‑1 (MQ=255)
       cATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGc                                                                                                                                                <  1:2093867/150‑1 (MQ=255)
       cATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGc                                                                                                                                                <  2:6310005/150‑1 (MQ=255)
        aTTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGGGTTGCGTAGCAAAGcc                                                                                                                                               >  1:1987985/1‑150 (MQ=255)
           ttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCgg                                                                                                                                            <  1:2762396/150‑1 (MQ=255)
            tCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGCCCGGa                                                                                                                                           <  2:7411569/150‑1 (MQ=255)
               cgcTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGGGTTTGTTAGCAAAACCCCGAAAc                                                                                                                                        >  1:3952012/1‑150 (MQ=18)
                  tctcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgg                                                                                                                                     <  1:6622759/150‑1 (MQ=255)
                   ctcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                    <  1:4188154/150‑1 (MQ=255)
                    tcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgggg                                                                                                                                   >  2:2422040/1‑150 (MQ=255)
                    tcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgggg                                                                                                                                   >  2:4866739/1‑150 (MQ=255)
                      aTAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAATCCCGGAAACTGGGGtt                                                                                                                                 >  1:6081891/1‑150 (MQ=255)
                         aGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCa                                                                                                                              <  1:4962004/150‑1 (MQ=255)
                         aGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCa                                                                                                                              <  2:3458116/150‑1 (MQ=255)
                            cACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCCAAGCCCGGAAACTGGGGTTTCAGCt                                                                                                                           >  2:2966496/1‑150 (MQ=255)
                            cACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCt                                                                                                                           <  2:161382/150‑1 (MQ=255)
                             aCTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACGGGGGTTTCAGcgg                                                                                                                          >  2:6349472/1‑148 (MQ=255)
                               tGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCtggt                                                                                                                        <  1:6242828/150‑1 (MQ=255)
                               tGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCtggt                                                                                                                        >  2:1074374/1‑150 (MQ=255)
                               tGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCATAGCAAAGCCCGGCAAAAGGGGTTTCAGCtggt                                                                                                                        >  2:2731816/1‑150 (MQ=17)
                                  gcATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgc                                                                                                                     <  1:5849324/150‑1 (MQ=255)
                                       aGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                                <  1:784525/150‑1 (MQ=255)
                                         ttGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTgt                                                                                                              >  2:1159571/1‑149 (MQ=255)
                                         ttGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGa                                                                                                              <  1:3071078/150‑1 (MQ=255)
                                         ttGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGa                                                                                                              <  2:405005/150‑1 (MQ=255)
                                         ttGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTTGCTGCTTGa                                                                                                              >  2:5191970/1‑150 (MQ=255)
                                            gTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGGGTTGGTTAGCCAAGCCCCGAAACAGGGGGTTCACCAGGTggcggcggtggtg                                                                                                           >  1:5252676/1‑143 (MQ=17)
                                              aCGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGtt                                                                                                         <  1:1169897/150‑1 (MQ=255)
                                                gTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTa                                                                                                       >  2:2846855/1‑150 (MQ=255)
                                                 tCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAg                                                                                                      <  1:2646164/150‑1 (MQ=255)
                                                   tGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTAAGCAGGTGGCTGCTTGAGTGTTTAGaa                                                                                                    >  2:6576636/1‑150 (MQ=255)
                                                    ggTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGaaa                                                                                                   <  1:4375392/150‑1 (MQ=255)
                                                       gtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGGTTTTGTTGGCAAATGCCGGAAAATTGGGTTTTAGGTGGTGGGTGGTTGTTTGTTTTGtaaaaa                                                                                                >  2:229188/1‑148 (MQ=17)
                                                          cAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATg                                                                                             <  2:556453/150‑1 (MQ=255)
                                                           aGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGGTTTGCTTAGCAAAGCCCCGAAAATGGGGGTTCAGCTGGGGGGTGCTTGTGTTTTTTTAAAAAAAttg                                                                                            >  1:1066916/1‑148 (MQ=17)
                                                           aGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATCCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATgg                                                                                            <  1:2618065/150‑1 (MQ=255)
                                                             ccGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGt                                                                                          <  1:2904306/150‑1 (MQ=255)
                                                             ccGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGACGAATTCCCAACCAGAAGCCAAGACCGGCATGACCGGCGTTTGCTTAGCAAAGCCCCGAAACATGTGGTTAAACAGGTGGATGCTGTATAGGTTAGAAAAAAATgggg                                                                                          >  2:903253/1‑149 (MQ=9)
                                                                    aCCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTc                                                                                   <  1:2895374/150‑1 (MQ=255)
                                                                    aCCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGc                                                                                                                                                >  2:3333954/1‑89 (MQ=14)
                                                                    aCCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGc                                                                                                                                                <  1:3333871/89‑1 (MQ=14)
                                                                      cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTg                                                                                 >  2:6202218/1‑150 (MQ=255)
                                                                      cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAAAGGGTTGCATTCTg                                                                                 >  2:584391/1‑150 (MQ=255)
                                                                      cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACGGGGGTTTCAGCTGGTGGCTGCGTGAGTGTTTAGAAAAAAATGGGTTGCATTCTg                                                                                 >  2:5252075/1‑150 (MQ=255)
                                                                      cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTAGCAAAACCCCGAAAACGGGGTTTCACCTTGTGGCTGCTTGAGTGTTTATAAAAAAAAGGGTTGGAttgag                                                                                 >  2:1440638/1‑147 (MQ=14)
                                                                      cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTGCTTCGCCAAGCCCCGAAACAGGGGTTTCAGCGGGGGGCGGCTGGAGTGTTTAGAAAAAAAAGGGGTGCAATCTg                                                                                 >  1:2987029/1‑150 (MQ=17)
                                                                       aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGt                                                                                <  2:6011760/150‑1 (MQ=255)
                                                                       aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTGCGTCGCCAAGCCCCGAAACAGGGGGTTCAGCTGGGGGCGGCGTGAGGGGGTAGAAAACAAATGGTTGCATTCGGt                                                                                >  1:7306329/1‑150 (MQ=14)
                                                                        agaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGtt                                                                               <  2:2481192/150‑1 (MQ=255)
                                                                        agaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGtt                                                                               <  2:7000106/150‑1 (MQ=255)
                                                                          aagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTg                                                                             >  2:3656801/1‑150 (MQ=255)
                                                                           agCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGa                                                                            <  2:3083678/150‑1 (MQ=255)
                                                                               ggCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTGTTAGCAAAGCCCCGCAAAAGGGGTTTAAGCAGGTGGGTGCGTGAGTTTTTTTAAAAAAAAGGGGTGGATTATGTTTTAACAg                                                                        >  1:446288/1‑150 (MQ=12)
                                                                               ggCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTTCTTCGCAAAGCCCCGAAACAGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTATAAAAAAAAGGGTTGGAATCTGTTTTACCAg                                                                        >  2:1883487/1‑150 (MQ=17)
                                                                                gCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAgg                                                                       <  1:5793374/150‑1 (MQ=255)
                                                                                gCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTTCTTCGCAAACCCCCGAAAAGTGGGTTTCAGCTGGTTGCAGCGTGAGTGTTTAGAAAACAAAGGGGTGCATTGTGGTTGAcaaag                                                                       >  2:79037/1‑146 (MQ=14)
                                                                                    aTGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCtt                                                                   <  1:4148529/150‑1 (MQ=255)
                                                                                        tGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTTCGTGGCAAAGCCCCGCAACAGGGGTTTCAGCTGGGGGCGGCGTGAGTGTTTATAAAACAATGGGTTGCATTCTGTTTGTACAGGGCGTACCt                                                               >  2:3304296/1‑150 (MQ=17)
                                                                                        tGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTGCTTCCCAAAGACCCGACACCGGGGGTTCAGCGTGTGGCGGCGTGGGGGTTTATAAAAAAAATGGGTGAAGTCTGTTTTGCCAGGGCGTaacg                                                               >  1:2297581/1‑147 (MQ=12)
                                                                                         gAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTCTTCGCAAAGCCCCGAAACTGGGGTTTCAGCTGGTGGCGGCTGGAGTGTTTAGAAAAAAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTc                                                              >  2:4498295/1‑150 (MQ=18)
                                                                                         gAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTc                                                              <  1:4738452/150‑1 (MQ=255)
                                                                                           aCAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGGATTGTGTTTGTCCAGGGCTTACCTCaa                                                            >  2:6144902/1‑150 (MQ=255)
                                                                                              aGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTa                                                         <  1:3131647/150‑1 (MQ=255)
                                                                                                gTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTAct                                                       <  2:3017130/150‑1 (MQ=255)
                                                                                                 tCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTActc                                                      >  2:6699663/1‑150 (MQ=255)
                                                                                                 tCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTActc                                                      <  2:6661630/150‑1 (MQ=255)
                                                                                                  cTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTActct                                                     <  1:2287071/150‑1 (MQ=255)
                                                                                                   ttGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTTCTTCGCAAAGCCCCGCAACAGGGGTTTAAGCAGGTGGCTGCGTGAGTGTTTAGAAAACAAATGGTTGCCTTATGTTTGAACACGGCTTAACACAAGTaagaga                                                    >  1:6689153/1‑145 (MQ=17)
                                                                                                     gAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATc                                                  <  2:928869/150‑1 (MQ=255)
                                                                                                     gAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATc                                                  <  1:3541894/150‑1 (MQ=255)
                                                                                                      aGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATcc                                                 <  1:6498556/150‑1 (MQ=255)
                                                                                                        gAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCaa                                               <  1:1769730/150‑1 (MQ=255)
                                                                                                         aaTTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGTAAACAAAGGGTTGCATTCTGTTTTGCCAGGGCGTACCACAAGTACTCGATACaaa                                              >  1:6172606/1‑150 (MQ=18)
                                                                                                           ttCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCGTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGa                                            >  2:4448586/1‑150 (MQ=255)
                                                                                                           ttCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACGGGGGTTTCAGCAGGTGGCTGCGTGAGTGTTTAGAAAACAATGGGTTGCAGTCGGTTTGTCCAGGGCGTACCTCAAGTACTCTATCCAAAGa                                            >  2:6959415/1‑150 (MQ=21)
                                                                                                                  ccAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaaa                                     <  1:3775097/150‑1 (MQ=255)
                                                                                                                  ccAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaaa                                     <  2:225259/150‑1 (MQ=255)
                                                                                                                     gAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGaa                                  <  1:3656801/150‑1 (MQ=255)
                                                                                                                           aaGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCa                            <  1:7049634/150‑1 (MQ=255)
                                                                                                                           aaGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAAAAAAGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATACAAAGAAAAAAAAGAACAATCa                            >  2:5463356/1‑150 (MQ=255)
                                                                                                                            aGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAg                           >  2:4957325/1‑150 (MQ=255)
                                                                                                                                    aGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCa                   <  2:5054371/150‑1 (MQ=255)
                                                                                                                                       tCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCa                <  2:2230884/150‑1 (MQ=255)
                                                                                                                                         cGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCaaa              >  1:6583049/1‑150 (MQ=255)
                                                                                                                                         cGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCaaa              >  1:5611319/1‑150 (MQ=255)
                                                                                                                                          gggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGTCCAGGGCTTACCTCAAGTACTCTCTCCAAAGACAAGCAAGAACACTCAGTGCCATCAGCacag             >  2:5325508/1‑147 (MQ=21)
                                                                                                                                          gggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaaa                                     <  1:5499101/126‑1 (MQ=255)
                                                                                                                                          gggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCaaaa             >  2:5680508/1‑150 (MQ=255)
                                                                                                                                          gggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCaaaa             >  2:4700223/1‑150 (MQ=255)
                                                                                                                                           ggCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAAt            >  1:2642569/1‑150 (MQ=255)
                                                                                                                                            gCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCCGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCCGCAAAAtt           >  2:7120457/1‑150 (MQ=255)
                                                                                                                                             cTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAAttt          <  2:2949854/150‑1 (MQ=255)
                                                                                                                                               ttGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGt        <  1:2965667/150‑1 (MQ=255)
                                                                                                                                               ttGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGt        <  2:211972/150‑1 (MQ=255)
                                                                                                                                                tGCTTCGCACAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTa       <  2:4488981/150‑1 (MQ=255)
                                                                                                                                                   ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACAg    >  1:4816438/1‑150 (MQ=255)
                                                                                                                                                     cGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTCCATTCTGATTGACCAGGGCTTACCTCAAGTACTCTATCCAAACAAAAAAAAGAACAATCACTGCAATCAGCAAAATTTGTACAGTc  <  2:2366285/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACAGTC  >  NC_002947/5557763‑5558061

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: