Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,265,714 T→A 5.5% L244Q (CTG→CAG)  folC → bifunctional folypolyglutamate synthase/dihydrofolate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,265,7140TA 5.5% 269.6 / 3.7 87L244Q (CTG→CAG) folCbifunctional folypolyglutamate synthase/dihydrofolate synthase
Reads supporting (aligned to +/- strand):  ref base T (37/45);  new base A (3/2);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTGATCTGGATCCGCCCGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGACCTCCGTAA  >  NC_002947/2265567‑2265861
                                                                                                                                                   |                                                                                                                                                     
gTGATCTGGATCCGCCCGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTgc                                                                                                                                                     <  1:790889/150‑1 (MQ=255)
gTGATCTGGATCCGCCCGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTgc                                                                                                                                                     <  2:40342/150‑1 (MQ=255)
  gagcTGGATCCGCCCGGGCCGTTGCGGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGGCGCTGATGGTGCGCAGGTCGAGCTgcgc                                                                                                                                                   <  2:36604/147‑1 (MQ=255)
   aTCTGGATCCGCCCGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGcgcg                                                                                                                                                  >  1:774362/1‑150 (MQ=255)
   aTCTGGATCCGCCCGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCAAAGAGCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGCTCGACCTgcgca                                                                                                                                                  <  1:635573/150‑2 (MQ=255)
       gtaggctcccGTGCCGTTGTTGGACAAGGCCGCCGAACTGGGCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGTCGCTGATGGGGCGCACGTCGCGCAGCTCGACCt                                                                                                                                              <  2:869956/143‑1 (MQ=255)
       ggATCCGCCCGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCNACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCt                                                                                                                                              >  1:187817/1‑150 (MQ=255)
               ccGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGc                                                                                                                                      <  2:457285/150‑1 (MQ=255)
                cGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCt                                                                                                                                     >  1:2010074/1‑150 (MQ=255)
                      gTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGTGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTgg                                                                                                                                   >  2:1751095/1‑146 (MQ=255)
                      gTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTgg                                                                                                                                   <  1:1751095/146‑1 (MQ=255)
                       ttGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGAc                                                                                                                                 >  2:1264550/1‑147 (MQ=255)
                       ttGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGAc                                                                                                                                 <  1:1264550/147‑1 (MQ=255)
                       ttGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTg                                                                                                                              >  1:1985042/1‑150 (MQ=255)
                            ggACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCAt                                                                                                                         <  1:1036338/150‑1 (MQ=255)
                            ggACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCAt                                                                                                                         <  2:1170638/150‑1 (MQ=255)
                             gACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGTTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATg                                                                                                                        >  2:749931/1‑150 (MQ=255)
                                  ggCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGaaaa                                                                                                                   <  2:2226051/150‑1 (MQ=255)
                                  ggCCGCCGAACTGGCCGCGCAAGTGTTCGTGCGCGGTCGGGAATTCGACCTTGCCAGTGCCGAAGGCTGCTGGGGCTCGCGTGGTACCGCCGCTGGTGGTGCGCAGGTCGAGCTGCGCGCCGTGCCTTTGCTGGACCTGCCCATGGaaaa                                                                                                                   <  2:187817/150‑1 (MQ=255)
                                       ccgAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCa                                                                                                              <  1:1465944/150‑1 (MQ=255)
                                             tGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCGGGCCAGTGCCGAATGCTGCTGGAGCTGGCGTAGTACTGCCGCTGGTGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTgg                                                                                                        >  2:495954/1‑150 (MQ=255)
                                               gCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTgc                                                                                                                                                     <  2:178400/103‑1 (MQ=255)
                                               gCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTgc                                                                                                                                                     >  1:178400/1‑103 (MQ=255)
                                                  gcgcCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTg                                                                                                   <  2:2004676/150‑1 (MQ=255)
                                                     ccACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGGCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCCGGGCGTGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAg                                                                                                >  1:782224/1‑150 (MQ=255)
                                                     ccACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAg                                                                                                <  2:852619/150‑1 (MQ=255)
                                                      cACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAgg                                                                                               >  1:718657/1‑150 (MQ=255)
                                                       aCTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGc                                                                                              >  2:664419/1‑150 (MQ=255)
                                                         tGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCt                                                                                                                                                       >  2:1810483/1‑91 (MQ=255)
                                                         tGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCt                                                                                                                                                       <  1:1810483/91‑1 (MQ=255)
                                                         tGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCtt                                                                                            <  1:2402838/150‑1 (MQ=255)
                                                         tGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCtt                                                                                            <  1:1167497/150‑1 (MQ=255)
                                                           ttCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑Tc                                                                                        >  1:464365/1‑150 (MQ=255)
                                                            tCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑Tcc                                                                                       <  2:774362/150‑1 (MQ=255)
                                                             ccTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCct                                                                                      <  2:2010074/150‑1 (MQ=255)
                                                              cTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTg                                                                                                         <  2:1822230/132‑1 (MQ=255)
                                                              cTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTg                                                                                                         >  1:1822230/1‑132 (MQ=255)
                                                              cTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGc                                                                                                       >  1:1191895/1‑134 (MQ=255)
                                                              cTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGc                                                                                                       <  2:1191895/134‑1 (MQ=255)
                                                              cTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCctg                                                                                     >  2:1959158/1‑150 (MQ=255)
                                                                        gggAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGc                                                                           >  1:1851447/1‑150 (MQ=255)
                                                                         ggAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCt                                                                          >  2:1185669/1‑150 (MQ=255)
                                                                           aaTTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCctgc                                                                                    <  2:2156119/138‑1 (MQ=255)
                                                                           aaTTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCctgc                                                                                    >  1:2156119/1‑138 (MQ=255)
                                                                              tCGACCTGGCCAGTGCCGAAGGCTGTTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGt                                                                     <  1:749931/150‑1 (MQ=255)
                                                                               cGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTg                                                                    <  2:1315807/150‑1 (MQ=255)
                                                                               cGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTg                                                                    <  1:17472/150‑1 (MQ=255)
                                                                                gACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTgg                                                                   >  2:2403513/1‑150 (MQ=255)
                                                                                 aCCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGa                                                                  >  2:835953/1‑150 (MQ=255)
                                                                                  ccTGGTCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGaa                                                                 <  1:25921/150‑1 (MQ=255)
                                                                                   cTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAAc                                                                >  2:979791/1‑150 (MQ=255)
                                                                                      gCCAGTCCCGAAGCCTGCTGGAGCTGCCCTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTTCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAAcgcc                                                             <  1:1189903/150‑2 (MQ=255)
                                                                                        cAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGTGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACgcggc                                                           >  2:2363523/1‑150 (MQ=255)
                                                                                          gTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGgcgc                                                         >  2:1392581/1‑150 (MQ=255)
                                                                                          ctGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGgcgc                                                         >  2:1657119/2‑150 (MQ=255)
                                                                                                      gctgGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAgg                                             <  1:2323977/150‑1 (MQ=255)
                                                                                                       ctgGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGgc                                            >  2:580753/1‑150 (MQ=255)
                                                                                                        tgGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGgcg                                           <  1:2019796/150‑1 (MQ=255)
                                                                                                         ggAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGgcgc                                          >  1:825664/1‑150 (MQ=255)
                                                                                                         ggAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGgcgc                                          <  2:718657/150‑1 (MQ=255)
                                                                                                            gCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCgctgc                                       <  2:2090605/150‑1 (MQ=255)
                                                                                                            gCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCCACCTGCCTTGGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCgctgc                                       <  1:1969726/150‑1 (MQ=255)
                                                                                                               ggCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTgg                                    <  2:1765009/150‑1 (MQ=255)
                                                                                                                gCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCaa                                                       <  1:1635130/130‑1 (MQ=255)
                                                                                                                gCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCaa                                                       >  2:1635130/1‑130 (MQ=255)
                                                                                                                gCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGa                                   <  1:1392581/150‑1 (MQ=255)
                                                                                                                         tGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATc                         >  2:1496942/1‑150 (MQ=255)
                                                                                                                              cTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTg                     <  2:464365/150‑1 (MQ=255)
                                                                                                                              cTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTg                     <  1:1077867/150‑1 (MQ=255)
                                                                                                                              cTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTg                     <  1:664419/150‑1 (MQ=255)
                                                                                                                              cTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCCAACAGGCTTAATACTGGAGAT‑‑GGTTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTg                     >  1:171923/1‑150 (MQ=255)
                                                                                                                               tGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTgg                    <  1:835953/150‑1 (MQ=255)
                                                                                                                                     aaggcAGGTCGCGCAGCTCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCaa                                                       >  1:2323081/4‑109 (MQ=255)
                                                                                                                                     aaggcAGGTCGCGCAGCTCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCaa                                                       <  2:2323081/106‑1 (MQ=255)
                                                                                                                                     aaggcAGGTCGCGCAGCTCGACCTGAAATTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGtct              >  1:1873492/4‑150 (MQ=255)
                                                                                                                                        gcAGGTCGCGCAGCTCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGa           >  2:951684/1‑150 (MQ=255)
                                                                                                                                        gcAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCa                                                                                                 >  2:1282114/1‑66 (MQ=255)
                                                                                                                                        gcAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCa                                                                                                 <  1:1282114/66‑1 (MQ=255)
                                                                                                                                        gcAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑Tc                                                                                        >  2:604094/1‑73 (MQ=255)
                                                                                                                                        gcAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑Tc                                                                                        <  1:604094/73‑1 (MQ=255)
                                                                                                                                        gcAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGAt                                                                                <  2:343053/81‑1 (MQ=255)
                                                                                                                                        gcAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGAt                                                                                >  1:343053/1‑81 (MQ=255)
                                                                                                                                        gcAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGa           <  1:2403513/150‑1 (MQ=255)
                                                                                                                                         cAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGAc          <  2:825664/150‑1 (MQ=255)
                                                                                                                                            gTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGACCTc       >  2:2076205/1‑150 (MQ=255)
                                                                                                                                               gAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCctgct                                                                                   >  2:598573/1‑71 (MQ=255)
                                                                                                                                               gAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCctgct                                                                                   <  1:598573/71‑1 (MQ=255)
                                                                                                                                                 gCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGACCTCCGTaa  >  2:1699076/1‑150 (MQ=255)
                                                                                                                                                 gCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATAACTGGGCGTCTCGACCCCCgaca  >  2:373482/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTGATCTGGATCCGCCCGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGACCTCCGTAA  >  NC_002947/2265567‑2265861

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: