Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,200,011 Δ1 bp 11.3% coding (834/1278 nt) PP_1944 → aminomethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,200,0110G.11.3% 288.8 / 32.1 97coding (834/1278 nt)PP_1944aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (35/51);  new base . (5/6);  total (40/57)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 105 total aligned reads displayed.

ACTGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTGTTC  >  NC_002947/2199862‑2200155
                                                                                                                                                      |                                                                                                                                                
aCTGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAg                                                                                                                                                  <  2:96854/150‑1 (MQ=255)
 cTGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAgc                                                                                                                                                 >  1:2007466/1‑150 (MQ=255)
    aCTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGCGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAgctgg                                                                                                                                              >  1:787205/1‑150 (MQ=255)
    aCTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAgctgg                                                                                                                                              >  1:223806/1‑150 (MQ=255)
      tACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGATTCGACTTTCAAGGTTGCGGGGAGATTTGAGTCGAATAAGGTAAGCGCTTGGTATCGCAGCCCAGTCGAgctgggc                                                                                                                                            >  2:465183/1‑150 (MQ=255)
        cATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAgctgggctg                                                                                                                                          >  1:1414841/1‑150 (MQ=255)
         aTGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAgctgggctgg                                                                                                                                         <  2:1187278/150‑1 (MQ=255)
         aTGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAgctgggctgg                                                                                                                                         <  2:1314732/150‑1 (MQ=255)
                   ttttCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCa                                                                                                                               >  2:404667/1‑150 (MQ=255)
                   ttttCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCa                                                                                                                               >  2:404657/1‑150 (MQ=255)
                    tttCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCAt                                                                                                                              <  2:267285/150‑1 (MQ=255)
                     ttCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATc                                                                                                                             <  2:1515330/150‑1 (MQ=255)
                       cGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCaa                                                                                                                           <  2:2007466/150‑1 (MQ=255)
                          cccAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGtt                                                                                                                        <  1:2366529/150‑1 (MQ=255)
                          cccAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAACtt                                                                                                                        <  1:2403141/150‑1 (MQ=255)
                           ccAGAGATGGCGGGTTATCTCAAGGACTTCCGCGCGGGTATGCCAAGTTTTGCCTCGCCTTTCAATGTTGCGGGGAGCTTTGAGTAGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTc                                                                                                                       >  2:260071/1‑150 (MQ=255)
                            cAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATc                                                                                                                             <  1:159942/143‑1 (MQ=255)
                            cAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATc                                                                                                                             >  2:159942/1‑143 (MQ=255)
                              gagaTGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGAt                                                                                                                    <  2:1414841/150‑1 (MQ=255)
                              gagaTGGCAGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGAt                                                                                                                    <  1:880319/150‑1 (MQ=255)
                                     cGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCGGCCCAGTCGAGCTGGGCTGGGCCaaaa                                                                                                                                  >  1:1705863/1‑129 (MQ=255)
                                     cGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCGGCCCAGTCGAGCTGGGCTGGGCCaaaa                                                                                                                                  <  2:1705863/129‑1 (MQ=255)
                                        gTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTa                                                                                                          <  1:404667/150‑1 (MQ=255)
                                        gTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTa                                                                                                          <  1:404657/150‑1 (MQ=255)
                                          tATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATc                                                                                                        <  2:2207564/150‑1 (MQ=255)
                                              tcAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCg                                                                                                       >  2:585951/1‑147 (MQ=255)
                                                aaGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGt                                                                                                  >  2:2003789/1‑150 (MQ=255)
                                                  ggAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAg                                                                                                                          >  2:623737/1‑124 (MQ=255)
                                                  ggAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAg                                                                                                                          <  1:623737/124‑1 (MQ=255)
                                                  ggAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACt                                                                                                             >  2:1622055/1‑137 (MQ=255)
                                                  ggAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGc                                                                                                <  1:1886540/150‑1 (MQ=255)
                                                  ggAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTCAAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACt                                                                                                             <  1:1622055/137‑1 (MQ=255)
                                                          gcgcgGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATTGGTCGTGGCGCGGTgg                                                                                        >  2:1722621/1‑150 (MQ=255)
                                                           cgcgGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGa                                                                                       >  2:855507/1‑150 (MQ=255)
                                                                  aTGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAagag                                                                                <  1:1137246/150‑1 (MQ=255)
                                                                     ccAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGcc                                                                                                                                      <  1:219848/93‑1 (MQ=255)
                                                                     ccAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGcc                                                                                                                                      >  2:219848/1‑93 (MQ=255)
                                                                     ccAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTAt                                                                                                         <  1:616196/122‑1 (MQ=255)
                                                                     ccAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTAt                                                                                                         >  2:616196/1‑122 (MQ=255)
                                                                          ttttGCCTCGACTTTCTATTTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCaga                                                                        <  2:514164/150‑1 (MQ=255)
                                                                              gCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCt                                                                    <  1:1948012/150‑1 (MQ=255)
                                                                              gCCTCGACTCTTAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCt                                                                    <  2:787205/150‑1 (MQ=255)
                                                                                 tCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGa                                                                 >  1:2261187/1‑150 (MQ=255)
                                                                                 tCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGa                                                                 >  2:103361/1‑150 (MQ=255)
                                                                                  cGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGAc                                                                <  1:1479954/150‑1 (MQ=255)
                                                                                    aCTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACgg                                                              <  1:855507/150‑1 (MQ=255)
                                                                                      tttCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGAt                                                                                                                    <  1:847018/94‑1 (MQ=255)
                                                                                      tttCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGAt                                                                                                                    >  2:847018/1‑94 (MQ=255)
                                                                                          aaTGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATg                                                        >  2:821122/1‑150 (MQ=255)
                                                                                           aTGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATgc                                                       <  2:223806/150‑1 (MQ=255)
                                                                                              ttGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGtt                                                                             <  2:2348573/125‑1 (MQ=255)
                                                                                              ttGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGtt                                                                             >  1:2348573/1‑125 (MQ=255)
                                                                                                   gggAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTg                                               <  2:1473133/150‑1 (MQ=255)
                                                                                                    ggAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGt                                                                                                                         >  1:309433/1‑75 (MQ=255)
                                                                                                        cTTTGCGTCGCATACGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTg                                          >  2:2154019/1‑150 (MQ=255)
                                                                                                        cTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTg                                          <  2:2261187/150‑1 (MQ=255)
                                                                                                         tttGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTgg                                         >  1:1919605/1‑150 (MQ=255)
                                                                                                           tGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGaa                                       <  1:2185092/150‑1 (MQ=255)
                                                                                                               tCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGTTGGGCCAAAACCATCAAGTTCTATCACCACTTTATCGGTCGTGCCGCGTTGGAAAAATAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGt                                   <  1:465183/150‑1 (MQ=255)
                                                                                                               tCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACgg                                                              >  1:216341/1‑123 (MQ=255)
                                                                                                               tCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACgg                                                              <  2:216341/123‑1 (MQ=255)
                                                                                                               tCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGt                                   >  2:1262771/1‑150 (MQ=255)
                                                                                                                cGAATAAGGTAAGCCATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtg                                  <  1:415347/150‑1 (MQ=255)
                                                                                                                  aaTAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtgat                                <  1:260071/150‑1 (MQ=255)
                                                                                                                    tAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtgatga                              >  1:253057/1‑150 (MQ=255)
                                                                                                                         tAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTg                         >  1:1671591/1‑150 (MQ=255)
                                                                                                                           aGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCa                       <  1:976878/150‑1 (MQ=255)
                                                                                                                            gCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAg                      <  2:1919605/150‑1 (MQ=255)
                                                                                                                               aTTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAagag                                                                                >  1:922189/1‑89 (MQ=255)
                                                                                                                               aTTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAagag                                                                                <  2:922189/89‑1 (MQ=255)
                                                                                                                               aTTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGaa                                       <  1:1525521/130‑1 (MQ=255)
                                                                                                                                   gTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGt                                                            >  1:708412/1‑105 (MQ=255)
                                                                                                                                   gTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGt                                                            <  2:708412/105‑1 (MQ=255)
                                                                                                                                   gTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGCTGACGTTGCAGATGTTCa               >  2:917558/1‑150 (MQ=255)
                                                                                                                                   gTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCa               <  1:1391747/150‑1 (MQ=255)
                                                                                                                                      tCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGAt                    >  2:2018565/1‑142 (MQ=255)
                                                                                                                                      tCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGAt                    <  1:2018565/142‑1 (MQ=255)
                                                                                                                                      tCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGc            <  1:2154019/150‑1 (MQ=255)
                                                                                                                                      gcccAGCCCAGCTCGACTGGGCTGCGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCAt              >  2:9000/4‑148 (MQ=255)
                                                                                                                                      gcccAGCCCAGCTCGACTGGGCTGCGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCAt              <  1:9000/145‑1 (MQ=255)
                                                                                                                                        gCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTc                                                                                                    >  1:1170670/1‑60 (MQ=255)
                                                                                                                                        gCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTc                                                                                                    <  2:1170670/60‑1 (MQ=255)
                                                                                                                                        ccAGCCCAGCTCGACTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAAt                                      <  2:1561041/121‑1 (MQ=255)
                                                                                                                                        ccAGCCCAGCTCGACTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAAt                                      >  1:1561041/2‑122 (MQ=255)
                                                                                                                                        ccAGCCCAGCTCGACTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCa               >  2:623881/2‑145 (MQ=255)
                                                                                                                                        ccAGCCCAGCTCGACTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCa               <  1:623881/144‑1 (MQ=255)
                                                                                                                                        ccAGCCCAGCTCGACTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCAt              >  2:554780/2‑146 (MQ=255)
                                                                                                                                        ccAGCCCAGCTCGACTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCAt              <  1:554780/145‑1 (MQ=255)
                                                                                                                                         cAGCCCAGCTCGACTGGGCTGGGCCAAAACCATCAAGtt                                                                                                                        <  1:321769/39‑1 (MQ=17)
                                                                                                                                         cAGCCCAGCTCGACTGGGCTGGGCCAAAACCATCAAGtt                                                                                                                        >  2:321769/1‑39 (MQ=17)
                                                                                                                                         cAGCCCAGCTCGACTGGGCTGCGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTc         <  2:784987/150‑1 (MQ=255)
                                                                                                                                         cAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTc         <  1:2003789/150‑1 (MQ=255)
                                                                                                                                           gCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGc       <  1:1722621/150‑1 (MQ=255)
                                                                                                                                            cccAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATg                                     <  1:2237694/119‑1 (MQ=255)
                                                                                                                                            cccAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATg                                     >  2:2237694/1‑119 (MQ=255)
                                                                                                                                            cccAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCt      >  1:1788860/1‑150 (MQ=255)
                                                                                                                                            cccAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCt      <  1:1262771/150‑1 (MQ=255)
                                                                                                                                               aGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCt          <  1:2114259/143‑1 (MQ=255)
                                                                                                                                               aGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCt          >  2:2114259/1‑143 (MQ=255)
                                                                                                                                                gTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGCTGACGTTGCAGATGTTCATGCCTCGCTGTTc  <  1:917558/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
ACTGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGTCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGTCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTGTTC  >  NC_002947/2199862‑2200155

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: