Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,265,717 G→T 6.1% R245L (CGC→CTC)  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,7170GT 6.1% 253.6 / 3.5 80R245L (CGC→CTC) folCbifunctional folypolyglutamate synthase/dihydrofolate synthase
Reads supporting (aligned to +/- strand):  ref base G (35/40);  new base T (3/2);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GATCTGGATCCGCCCGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGACCTCCGTAA  >  NC_002947/2265569‑2265861
                                                                                                                                                    |                                                                                                                                                  
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)
                                                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)
                                                       tGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCtt                                                                                            <  1:1167497/150‑1 (MQ=255)
                                                       tGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCtt                                                                                            <  1:2402838/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                                                                    <  1:17472/150‑1 (MQ=255)
                                                                             cGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTg                                                                    <  2:1315807/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                                                       >  2:1635130/1‑130 (MQ=255)
                                                                                                              gCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCaa                                                       <  1:1635130/130‑1 (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                                                       <  2:2323081/106‑1 (MQ=255)
                                                                                                                                   aaggcAGGTCGCGCAGCTCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCaa                                                       >  1:2323081/4‑109 (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                                                                                >  1:343053/1‑81 (MQ=255)
                                                                                                                                      gcAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGAt                                                                                <  2:343053/81‑1 (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)
                                                                                                                                                   cgcgACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATgggg                                                                            >  1:917587/1‑72 (MQ=255)
                                                                                                                                                   cgcgACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATgggg                                                                            <  2:917587/72‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GATCTGGATCCGCCCGGGCCGTTGCTGGACAAGGCCGCCGAACTGGCCGCGCCACTGTTCCTGCGCGGTCGGGAATTCGACCTGGCCAGTGCCGAAGGCTGCTGGAGCTGGCGTGGTACTGCCGCTGATGGTGCGCAGGTCGAGCTGCGCGACCTGCCTTTGCTGGACCTGCCCATGGAAAACGCCACGCTGGCCTTGCAGGCTT‑‑TCCTGCTGATGGGGCTGCCGTGGAACGCGGCGCAAGTCCGCCAGGCGCTGCTGGATACCCGTATCACTGGGCGTCTCGACCTCCGTAA  >  NC_002947/2265569‑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: