Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,485,018 C→G 7.4% A363A (GCC→GCG yhdX → amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,485,0180CG 7.4% 345.2 / 15.3 108A363A (GCC→GCGyhdXamino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (52/48);  new base G (5/3);  total (57/51)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 112 total aligned reads displayed.

TGCGGGTGATCATTCCGCCGCTGACCAGCCAGTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTCAAACCCGATATGCCGCCACCGGTGAAAACC  >  NC_002947/1484870‑1485166
                                                                                                                                                    |                                                                                                                                                    
tGCGGGTGATCATTCCGCCGCTGACCAGCCAGTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCa                                                                                                                                                     <  1:831208/150‑1 (MQ=255)
    ggTGATCATTCCGCCGCTGACCAGCCAGTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCAc                                                                                                                                                 >  1:2278845/1‑150 (MQ=255)
     gTGATCATTCCGCCGCTGACCAGCCAGTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCAcc                                                                                                                                                >  2:282365/1‑150 (MQ=255)
     gTGATCATTCCGCCGCTGACCAGCCAGTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCAcc                                                                                                                                                >  2:2210727/1‑150 (MQ=255)
       gATCATTCCGCCGCTGGCCAGCCAGTACCTGAGCCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCAt                                                                                                                                              <  2:1052547/150‑1 (MQ=255)
         tCATTCCGCCGCTGACCAGCCAGTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGa                                                                                                                                            >  1:1361569/1‑150 (MQ=255)
          cATTCCGCCGCTAACCAGCCAGTGCCTGAACCTGGCGAAAAACTCTTCGCTGGCGGCCGGTATCGGCGACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACGAGACCGGCCGGGCCAGCGAGGTGATTGCCATCAACATTAg                                                                                                                                           >  2:2222/1‑150 (MQ=255)
              ccgccgCTGACCAGCCAGTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTc                                                                                                                                       <  2:2029273/150‑1 (MQ=255)
               ccccacggACCAGGCAGTACCTGCGCCTTGCGCGGAGCTCTTCGCTCGCGGACAGTATCGGCTACCCCGAGATAGTCTCGCTTTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCt                                                                                                                                      <  2:1015645/143‑1 (MQ=255)
                gccgcTGACCAGCCAGTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTa                                                                                                                                     >  1:1124444/1‑150 (MQ=255)
                 ccgcTGACCAGCCAGTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTAt                                                                                                                                    <  2:1145951/150‑1 (MQ=255)
                     tGACCAGCCAGTACCTGAACCTGGCTAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCg                                                                                                                                <  1:282365/150‑1 (MQ=255)
                     tGACCAGCCAGTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCg                                                                                                                                <  1:993580/150‑1 (MQ=255)
                      gACCAGCCAGTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGc                                                                                                                               >  1:924187/1‑150 (MQ=255)
                         cagccagTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCa                                                                                                                                               >  1:485151/1‑131 (MQ=255)
                         cagccagTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCa                                                                                                                                               <  2:485151/131‑1 (MQ=255)
                         cagccagTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCcat                                                                                                                            <  2:1198191/150‑1 (MQ=255)
                         cagccagTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCcat                                                                                                                            >  2:1101342/1‑150 (MQ=255)
                         cagccagTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTATTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCa                                                                                                                                                     >  2:1573392/1‑125 (MQ=255)
                         cagccagTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTATTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCa                                                                                                                                                     <  1:1573392/125‑1 (MQ=255)
                            ccagTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCcatcag                                                                                                                         >  2:2125460/1‑150 (MQ=255)
                              agTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCcatcagca                                                                                                                       <  2:127369/150‑1 (MQ=255)
                               gTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCcatcagcat                                                                                                                      <  1:852447/150‑1 (MQ=255)
                                 aCCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCACCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCcatcagcatca                                                                                                                    <  1:643025/150‑1 (MQ=255)
                                        ccTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTc                                                                                                             <  1:1794142/150‑1 (MQ=255)
                                          tGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCgc                                                                                                           >  2:2281447/1‑150 (MQ=255)
                                          cagCGAAGAGTTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCgc                                                                                                           >  2:2368551/3‑150 (MQ=255)
                                          cagCGAAGAGTTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAACATTTCgc                                                                                                           >  2:818807/3‑150 (MQ=255)
                                              gaagaaCTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGCTCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCtctgct                                                                                                       >  2:1388584/1‑150 (MQ=255)
                                               aagaaCTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGctgctg                                                                                                      >  1:1147558/1‑150 (MQ=255)
                                                 gaaCTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTTGCCAGCAGCATCAGCATTTCGCTGCtgat                                                                                                    >  2:1598381/1‑150 (MQ=255)
                                                 gaaCTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCgccgcttgt                                                                                                    <  2:2278845/150‑4 (MQ=255)
                                                   gttctTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCtgatga                                                                                                  >  2:895991/3‑150 (MQ=255)
                                                     tctTCGCTGGCGGCCGGTATCGGCTACCCTTAGATGGTCTCGCTCTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTTTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAAc                                                                                                <  1:375794/150‑1 (MQ=255)
                                                     tctTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAAc                                                                                                <  2:1967475/150‑1 (MQ=255)
                                                       ttCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTg                                                                                              <  1:174088/150‑1 (MQ=255)
                                                       ttCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTg                                                                                              >  1:2123279/1‑150 (MQ=255)
                                                       ttCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACAAGCCCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGCCTACACCAGGATCAGCATCAGCATTTCGCTGCTGATGAACTg                                                                                              <  1:1999467/150‑1 (MQ=255)
                                                           cTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTac                                                                                          <  2:2135953/150‑1 (MQ=255)
                                                             ggcggcCGGTATCGGCTACCCGGAGATGGTCTCGCTAAAATCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTacaa                                                                                        <  1:1833507/150‑1 (MQ=255)
                                                                ggcCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTAcaacaa                                                                                     >  2:294707/1‑150 (MQ=255)
                                                                 gcCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAg                                                                                    <  2:2123279/150‑1 (MQ=255)
                                                                     gTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCCACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCCTGAGTGTCTATCTCGCCATCAGCATCAGCATGTCGCTGCTGATGAACTGGTACAACAAGCGCa                                                                                <  1:49481/150‑1 (MQ=255)
                                                                              aCCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGa                                                                       <  1:895991/150‑1 (MQ=255)
                                                                                ccGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTa                                                                                                                                     <  2:1708430/86‑1 (MQ=255)
                                                                                ccGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTa                                                                                                                                     >  1:1708430/1‑86 (MQ=255)
                                                                                ccGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATc                                                                     >  2:1590394/1‑150 (MQ=255)
                                                                                 cGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCAATGCGCttgttg                                                                    <  1:2125460/150‑6 (MQ=255)
                                                                                  ggAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGa                                                                   >  2:2204243/1‑150 (MQ=255)
                                                                                         gTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGa                                                            <  1:2368551/150‑1 (MQ=255)
                                                                                         gTCTAGCTGTTCGCCGGCACAGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGa                                                            <  2:1124444/150‑1 (MQ=255)
                                                                                          tctcGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAgc                                                                                   <  2:1215850/126‑1 (MQ=255)
                                                                                          tctcGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAgc                                                                                   >  1:1215850/1‑126 (MQ=255)
                                                                                          tctcGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAg                                                           <  2:72155/150‑1 (MQ=255)
                                                                                            tcGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGa                                                         >  1:589851/1‑150 (MQ=255)
                                                                                            tcGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGa                                                         >  1:2226493/1‑150 (MQ=255)
                                                                                                     gCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGc                                                    >  1:1460227/1‑146 (MQ=255)
                                                                                                     gCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGt                                                <  2:210224/150‑1 (MQ=255)
                                                                                                             cGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCcatca                                                                                                                          >  1:1925490/1‑68 (MQ=255)
                                                                                                             cGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCcatca                                                                                                                          <  2:1925490/68‑1 (MQ=255)
                                                                                                             cGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGccc                                        >  2:1499745/1‑150 (MQ=255)
                                                                                                                 cTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCt                                                                                                                                      >  1:196727/1‑52 (MQ=255)
                                                                                                                 cTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCt                                                                                                                                      <  2:196727/52‑1 (MQ=255)
                                                                                                                  tCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGtt                                   >  1:374782/1‑150 (MQ=255)
                                                                                                                  tCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGAGGTGATGATACGCCCGTGAATGCCCATGtt                                   >  1:1812229/1‑150 (MQ=255)
                                                                                                                    aaCCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTAcaacaa                                                                                     <  2:1694639/98‑1 (MQ=255)
                                                                                                                    aaCCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTAcaacaa                                                                                     >  1:1694639/1‑98 (MQ=255)
                                                                                                                    aaCCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGtttt                                 >  2:300596/1‑150 (MQ=255)
                                                                                                                     aCCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGa                                                            <  1:2036310/122‑1 (MQ=255)
                                                                                                                     aCCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGa                                                            >  2:2036310/1‑122 (MQ=255)
                                                                                                                     aCCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTc                                <  2:374782/150‑1 (MQ=255)
                                                                                                                     aCCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTc                                <  2:1822791/150‑1 (MQ=255)
                                                                                                                      ccAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTAt                                                                                                                                    >  1:41608/1‑49 (MQ=255)
                                                                                                                      ccAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTAt                                                                                                                                    <  2:41608/49‑1 (MQ=255)
                                                                                                                        aggccGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCTTTGCGCTGATCGAGCGGTGAGGATACGCGCGTGAATGCCCATGTTTTCaaa                             <  2:1361569/147‑1 (MQ=255)
                                                                                                                        aGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTAcaacaa                                                                                     >  2:6689/1‑94 (MQ=255)
                                                                                                                        aGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTAcaacaa                                                                                     <  1:6689/94‑1 (MQ=255)
                                                                                                                        aGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGt                                                <  1:2370880/131‑1 (MQ=255)
                                                                                                                        aGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGt                                                >  2:2370880/1‑131 (MQ=255)
                                                                                                                         gACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAAATGTTACAACAAGCGCATTTCGCTGATCGAGCGGTGAGGATACGCCAGTGAATGACCATGTTTTCAAAc                            >  1:1978081/1‑150 (MQ=255)
                                                                                                                          aCCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTCAAAc                            <  2:716818/149‑1 (MQ=255)
                                                                                                                          aCCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTCAAAc                            >  1:716818/1‑149 (MQ=255)
                                                                                                                                  ggccaTCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAg                                                           <  2:505929/110‑1 (MQ=255)
                                                                                                                                  ggccaTCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAg                                                           >  1:505929/1‑110 (MQ=255)
                                                                                                                                  ggccaTCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTCAAACCCGATATgc                   >  1:50386/1‑150 (MQ=255)
                                                                                                                                      aTCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTCAAACCCGATATgccgcc               >  1:114639/1‑150 (MQ=255)
                                                                                                                                       tcatggTGATGGCGATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGc                                                                 <  2:1926953/95‑1 (MQ=255)
                                                                                                                                       tcatggTGATGGCGATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTCaa                              >  1:2088442/5‑134 (MQ=255)
                                                                                                                                       tcatggTGATGGCGATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTCaa                              <  1:2038760/130‑1 (MQ=255)
                                                                                                                                       tcatggTGATGGCGATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTCaa                              >  2:2038760/5‑134 (MQ=255)
                                                                                                                                       tCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTCAAACCCGATATgccgccg              >  1:1153213/1‑149 (MQ=255)
                                                                                                                                        catggTGATGGCGATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTCAAACCCGATATGCCGCCAc             >  2:2012010/4‑150 (MQ=255)
                                                                                                                                         atggTGATGGCGATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTCAAACCCGATATGCCGCCAcc            >  2:1382624/3‑150 (MQ=255)
                                                                                                                                           ggTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGccc                                        <  1:1803725/120‑1 (MQ=255)
                                                                                                                                           ggTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGccc                                        >  2:1803725/1‑120 (MQ=255)
                                                                                                                                              gATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTCAAACCCGATATGCCGCCACCGGTGa       >  2:266944/1‑150 (MQ=255)
                                                                                                                                                tCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTCAAACCCGATATGCCGCCACCGGTGaaa     >  2:591777/1‑150 (MQ=255)
                                                                                                                                                   ccATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGtt                                   >  1:909073/1‑117 (MQ=255)
                                                                                                                                                   ccATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTCAAACCCGATATGCCGCCACCGGTGAAAAcc  >  1:1240376/1‑150 (MQ=255)
                                                                                                                                                   ccATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTCAAACCCGATATGCCGCCACCGGTGAAAAcc  >  1:286663/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCGGGTGATCATTCCGCCGCTGACCAGCCAGTACCTGAACCTGGCGAAGAACTCTTCGCTGGCGGCCGGTATCGGCTACCCGGAGATGGTCTCGCTGTTCGCCGGTACCGTGCTCAACCAGACCGGCCAGGCCATCGAGGTGATCGCCATCACCATGAGTGTCTATCTCGCCATCAGCATCAGCATTTCGCTGCTGATGAACTGGTACAACAAGCGCATTGCGCTGATCGAGCGGTGAGGATACGCCCGTGAATGCCCATGTTTTCAAACCCGATATGCCGCCACCGGTGAAAACC  >  NC_002947/1484870‑1485166

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