Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 566,623 C→T 5.2% intergenic (+32/+41) bfr‑I → / ← uvrA bacterioferritin/uvrABC excinuclease protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947566,6230CT 5.2% 351.6 / 6.4 117intergenic (+32/+41)bfr‑I/uvrAbacterioferritin/uvrABC excinuclease protein A
Reads supporting (aligned to +/- strand):  ref base C (59/51);  new base T (2/4);  total (61/56)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.34e-01

Only 100 of 117 total aligned reads displayed.

ATCCTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCAC  >  NC_002947/566478‑566767
                                                                                                                                                    |                                                                                                                                                 
aTCCTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGTAccg                                                                                                                                               >  1:2064466/1‑149 (MQ=255)
  ccTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATc                                                                                                                                            >  2:2096634/1‑150 (MQ=255)
  ccTGCGTGCGCAGCTGCCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCCGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATc                                                                                                                                            >  2:127965/1‑150 (MQ=255)
    tgcgtgcgCAGCTGGCCGACACCGAAGAAGATCACACATACTGGCCAAACACACGGAGGACGCTGTTC‑‑GCCTGTCAGCCCCGTTAAGACTACCTGCAGTCGCAGACGTAAGTTATCCACAGGCTCAAGTAAAAGCCCCTGCCAGCCTCgg                                                                                                                                          >  1:2677533/1‑150 (MQ=255)
            cAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAggg                                                                                                                                  <  1:1733283/150‑1 (MQ=255)
              gCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGTACCGATGGTGCCAGGGGc                                                                                                                                <  2:2064466/150‑1 (MQ=255)
                  gCCGACACCGAAGAACATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCtttt                                                                                                                            <  2:2586248/150‑1 (MQ=255)
                   ccGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGTACCGATGGTGCCAGGGGCtttta                                                                                                                           <  1:2735271/150‑2 (MQ=39)
                   ccGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCttttt                                                                                                                           <  1:1439818/150‑1 (MQ=255)
                    cGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTc                                                                                                                          <  2:2887807/150‑1 (MQ=255)
                      acacCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCAt                                                                                                                        >  2:639737/1‑150 (MQ=255)
                       cacCGAAGAAGATCACACCTACTGGCTTGACAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCAtt                                                                                                                       >  1:2092987/1‑150 (MQ=255)
                            aagaagATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGTACCGATGGTGCCAggg                                                                                                                                  <  2:3015874/134‑1 (MQ=255)
                            aagaagATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGcc                                                                                                                  >  1:2386344/1‑150 (MQ=255)
                              gaagaTCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCgg                                                                                                                <  1:639737/150‑1 (MQ=255)
                               aagaTCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGc                                                                                                               >  2:2196120/1‑150 (MQ=255)
                                 gaTCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTgcgc                                                                                                                   >  2:109035/1‑144 (MQ=255)
                                 gaTCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTgcgc                                                                                                                   <  1:109035/144‑1 (MQ=255)
                                   tCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATa                                                                                                           >  1:1789183/1‑150 (MQ=255)
                                   tCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATa                                                                                                           >  2:1149985/1‑150 (MQ=255)
                                       acCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AACAATACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCa                                                                                                       >  2:2281550/1‑150 (MQ=255)
                                         cTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAgg                                                                                                     <  2:1849363/150‑1 (MQ=255)
                                          tACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGc                                                                                                    <  1:127965/150‑1 (MQ=255)
                                          tACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGc                                                                                                    <  2:1749171/150‑1 (MQ=255)
                                          tACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGc                                                                                                    <  2:1789183/150‑1 (MQ=255)
                                               gCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCaa                                                                                                             >  2:1243237/1‑136 (MQ=255)
                                               gCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCaa                                                                                                             <  1:1243237/136‑1 (MQ=255)
                                               gCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTAt                                                                                               >  2:1421450/1‑150 (MQ=255)
                                                cTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATc                                                                                              <  1:417388/150‑1 (MQ=255)
                                                 ttGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCgg                                                                                                                                          >  1:1952533/1‑105 (MQ=255)
                                                 ttGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATc                                                                                              <  2:628231/149‑1 (MQ=255)
                                                 ttGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATc                                                                                              >  1:628231/1‑149 (MQ=255)
                                                 ttGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTAATGCCCTAtac                                                                                             >  2:3037/1‑148 (MQ=255)
                                                    agaAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGc                                                                                          <  1:2188218/150‑1 (MQ=255)
                                                      aaGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTc                                                                                        <  2:2092987/150‑1 (MQ=255)
                                                               ggCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATgg                                                                               >  1:1313129/1‑150 (MQ=255)
                                                                 ccTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCa                                                                                                                                     <  2:1666527/94‑1 (MQ=255)
                                                                    tGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATTCACAGACTCAAGTAAAAGCCACTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCAGGCAATACTCAGGCCCTATCCATCTCCAGGAATGGTTTc                                                                           >  2:2835751/1‑150 (MQ=255)
                                                                     gATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCa                                                                          >  2:325103/1‑150 (MQ=255)
                                                                     gATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGAAGATGTAAGTTACCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCa                                                                          <  1:214401/150‑1 (MQ=255)
                                                                      aTCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAg                                                                                                      <  2:2226885/121‑1 (MQ=255)
                                                                      aTCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAg                                                                                                      >  1:2226885/1‑121 (MQ=255)
                                                                      aTCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAg                                                                         <  1:325103/150‑1 (MQ=255)
                                                                         aaGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTa                                                                      >  1:1140611/1‑150 (MQ=255)
                                                                           ggCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTCAt                                                                    >  1:510005/1‑150 (MQ=255)
                                                                              cATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTAt                                                                                               >  1:705104/1‑120 (MQ=255)
                                                                              cATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTAt                                                                                               <  2:705104/120‑1 (MQ=255)
                                                                              cATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGc                                                                 >  2:1537896/1‑150 (MQ=255)
                                                                                    gCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGtat                                                           <  1:1537896/150‑1 (MQ=255)
                                                                                      tGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTatat                                                         >  1:962238/1‑150 (MQ=255)
                                                                                      tGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTatat                                                         <  1:1421450/150‑1 (MQ=255)
                                                                                               cTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCt                                                  <  1:2128265/150‑1 (MQ=255)
                                                                                                  ccTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCa                                               >  2:1546535/1‑150 (MQ=255)
                                                                                                    gccAGTCGCAGATGTAAGTTATCACCCGACCCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCAtc                                             <  1:3037/148‑1 (MQ=255)
                                                                                                       aGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCg                                          >  1:2741026/1‑150 (MQ=255)
                                                                                                         tCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCt                                        >  2:308628/1‑150 (MQ=255)
                                                                                                           gCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTaa                                                                     >  1:1984715/1‑119 (MQ=255)
                                                                                                           gCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCtt                                                 >  1:1248250/1‑139 (MQ=255)
                                                                                                           gCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCtt                                                 <  2:1248250/139‑1 (MQ=255)
                                                                                                           gCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTCTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTaa                                                                     <  2:1984715/119‑1 (MQ=255)
                                                                                                            cAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCg                                                               >  2:1521238/1‑124 (MQ=255)
                                                                                                            cAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCg                                                               <  1:1521238/124‑1 (MQ=255)
                                                                                                             aGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATAGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGc                                    >  1:3028155/1‑150 (MQ=255)
                                                                                                                 gTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGt                                                             >  1:1713333/1‑121 (MQ=255)
                                                                                                                 gTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGt                                                             <  2:1713333/121‑1 (MQ=255)
                                                                                                                 gTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCtctt                                >  2:474543/1‑150 (MQ=255)
                                                                                                                    aGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTg                             <  1:308628/150‑1 (MQ=255)
                                                                                                                    aGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTg                             <  2:1140611/150‑1 (MQ=255)
                                                                                                                    aGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCCTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGCCTGCTTCATCTCGCTCAGCTCTTCTg                             <  2:2677533/150‑1 (MQ=255)
                                                                                                                      ttATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGgt                           <  1:474543/150‑1 (MQ=255)
                                                                                                                           cacaGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCa                                                                                                                         >  2:2433691/1‑51 (MQ=255)
                                                                                                                           cacaGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCa                                                                                                                         <  1:2433691/51‑1 (MQ=255)
                                                                                                                           cacaGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTAcc                      >  1:99090/1‑150 (MQ=255)
                                                                                                                            acaGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGt                                                                              <  1:883127/93‑1 (MQ=255)
                                                                                                                            acaGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCg                     >  1:1044009/1‑150 (MQ=255)
                                                                                                                            acaGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCg                     >  1:1926967/1‑150 (MQ=255)
                                                                                                                             caGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGc                    <  1:2196120/150‑1 (MQ=255)
                                                                                                                                aCTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTgc                                                                                                                     <  2:568833/50‑1 (MQ=255)
                                                                                                                                aCTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTgc                                                                                                                     >  1:568833/1‑50 (MQ=255)
                                                                                                                                aCTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAgg                 >  2:1053482/1‑150 (MQ=255)
                                                                                                                                  tCAAGTAAAAGCCCATGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCa               >  1:541239/1‑150 (MQ=255)
                                                                                                                                   cAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATAGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCaa              <  2:3028155/150‑1 (MQ=255)
                                                                                                                                    aaGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATgg                                                                  <  1:1872879/97‑1 (MQ=255)
                                                                                                                                    aaGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATgg                                                                  >  2:1872879/1‑97 (MQ=255)
                                                                                                                                    aaGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGACCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGgt                           >  2:1580441/1‑136 (MQ=255)
                                                                                                                                    aaGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGACCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGgt                           <  1:1580441/136‑1 (MQ=255)
                                                                                                                                     tgaaaaaGCCCCTGGTACCGATGGTGCCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATg                                                                                >  1:2136888/4‑82 (MQ=18)
                                                                                                                                     aGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATg            >  1:1814347/1‑150 (MQ=255)
                                                                                                                                       tAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGAt          <  2:1313129/150‑1 (MQ=255)
                                                                                                                                        aaaaGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATc         <  2:2273846/150‑1 (MQ=255)
                                                                                                                                         aaaGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATgg                                                                  <  2:2048112/92‑1 (MQ=255)
                                                                                                                                         aaaGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATgg                                                                  >  1:2048112/1‑92 (MQ=255)
                                                                                                                                          aaGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAAt             >  1:315321/1‑144 (MQ=255)
                                                                                                                                          aaGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAAt             <  2:315321/144‑1 (MQ=255)
                                                                                                                                          aaGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTg       <  2:2126058/150‑1 (MQ=255)
                                                                                                                                             ccccTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGccc    >  2:2155260/1‑150 (MQ=255)
                                                                                                                                              cccTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCa   >  2:2496886/1‑150 (MQ=255)
                                                                                                                                               ccTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCAc  >  1:1706937/1‑150 (MQ=255)
                                                                                                                                                cTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTg                             <  1:2339999/122‑1 (MQ=255)
                                                                                                                                                cTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTg                             >  2:2339999/1‑122 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATCCTGCGTGCGCAGCTGGCCGACACCGAAGAAGATCACACCTACTGGCTTGAGAAGCAGCAGGGCCTGATCAAGGCCATCGGGCTGG‑‑AAAACTACCTGCAGTCGCAGATGTAAGTTATCCACAGACTCAAGTAAAAGCCCCTGGCACCATCGGTACCAGGGGCTTTTTCATTGCGCCGGCAATACTCAGGCCCTATCCCGCTCCAGCAATGGTTTCAGGTAATGGCCGGTATATGACTGCTTCATCTCGCTCAGCTCTTCTGGTGTACCGCAGGCAATGATCTGCCCAC  >  NC_002947/566478‑566767

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: