Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,445,740 C→G 5.3% S216R (AGC→AGG thiD → phosphomethylpyrimidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,445,7400CG 5.3% 235.5 / 2.3 73S216R (AGC→AGGthiDphosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base C (35/34);  new base G (2/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CCGTTCTGTAGACACCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTGCTCC  >  NC_002947/5445591‑5445887
                                                                                                                                                     |                                                                                                                                                   
ccGTTCTGTAGACACCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAgc                                                                                                                                                     >  2:1363888/1‑150 (MQ=255)
  gTTCTGTAGACAACTGCTGATTACGGGCAGTCACGGCGACTATGATGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGCCCTCCCAGCGCCTGCTGGGCAACTACCATGCCTCGGGCTGCACCCTGGCCAccgc                                                                                                                                                   >  2:1020246/1‑150 (MQ=255)
   ttCTGTAGACACCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGcc                                                                                                                                                  >  1:34146/1‑150 (MQ=255)
   ttCTGTAGACACCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGcc                                                                                                                                                  >  2:294389/1‑150 (MQ=255)
    tCTGTAGACACCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGccc                                                                                                                                                 <  2:1405101/150‑1 (MQ=255)
        tAGACACCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTggc                                                                                                                                             <  2:957398/150‑1 (MQ=255)
            cacCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCCGGCCAGGGCGCTggccagg                                                                                                                                         <  2:1409785/150‑4 (MQ=255)
             acCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTgg                                                                                                                                              <  2:841740/144‑1 (MQ=255)
             acCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTgg                                                                                                                                              >  1:841740/1‑144 (MQ=255)
                      ttACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGgcgc                                                                                                                               >  1:561327/1‑150 (MQ=255)
                             cggTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGa                                                                                                                        <  2:34146/150‑1 (MQ=255)
                             cggTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGa                                                                                                                        <  1:1184617/150‑1 (MQ=255)
                                      cgacgaAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGTGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGa                                                                                                               <  2:12267/150‑1 (MQ=255)
                                               cgaaATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGt                                                                                                      >  2:1127397/1‑150 (MQ=255)
                                                      cacaACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCAcc                                                                                               >  2:1232682/1‑150 (MQ=255)
                                                      cacaACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTAACCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCAcc                                                                                               >  1:1387648/1‑150 (MQ=255)
                                                           ccGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCt                                                                                          >  2:597316/1‑150 (MQ=255)
                                                            cGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTg                                                                                         <  2:561327/150‑1 (MQ=255)
                                                             gCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTgg                                                                                        <  1:493735/150‑1 (MQ=255)
                                                                 gTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTa                                                                                    <  1:862507/150‑1 (MQ=255)
                                                                   acacCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTaca                                                                                  <  1:823405/150‑1 (MQ=255)
                                                                     acCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACAcc                                                                                >  1:820872/1‑150 (MQ=255)
                                                                       ccAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTg                                                                              >  1:996540/1‑150 (MQ=255)
                                                                       ccAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTg                                                                              <  2:634040/150‑1 (MQ=255)
                                                                                ccAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCt                                                                     >  1:491042/1‑150 (MQ=255)
                                                                                 cAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCAAACGCACCGCCCTGTACTACACCTGGCGCACCCTg                                                                    >  2:1083938/1‑150 (MQ=255)
                                                                                    cacttcacgtCGACTTGCCAGCGCCTGCCGGGCGGCTAGCATGTCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGt                                                                 <  1:40095/150‑1 (MQ=255)
                                                                                       ttcacttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGAc                                                              >  2:1069948/1‑150 (MQ=255)
                                                                                       ttcacttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGAc                                                              >  1:1313297/1‑150 (MQ=255)
                                                                                       ttcacttGGACTTCCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGAc                                                              <  1:294389/150‑1 (MQ=255)
                                                                                        tcacttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACg                                                             >  2:91666/1‑150 (MQ=255)
                                                                                          acttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGGCGGGGTCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGcc                                                           >  1:885558/1‑150 (MQ=255)
                                                                                          acttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGcc                                                           <  1:1363888/150‑1 (MQ=255)
                                                                                          acttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGcc                                                           <  2:1083290/150‑1 (MQ=255)
                                                                                            ttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGa                                                         >  2:482132/1‑150 (MQ=255)
                                                                                            ttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGa                                                         <  1:1069948/150‑1 (MQ=255)
                                                                                                  ttGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACt                                                   >  1:261216/1‑150 (MQ=255)
                                                                                                    gCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTgg                                                 <  1:138080/150‑1 (MQ=255)
                                                                                                    gCCAGCGCCTGCCGGGCAGCTACCACGGCTCCGCCTGCACTGTGGCCAGCGCCCTGGCCGGCCGCCTTGAGCTGGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTgg                                                 <  1:672152/150‑1 (MQ=255)
                                                                                                     ccAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTggg                                                >  2:1163204/1‑150 (MQ=255)
                                                                                                     ccAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTCGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTggg                                                >  2:377082/1‑150 (MQ=255)
                                                                                                     ccAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCCGGCCAGGGCGCTggccagg                                                                                                                                         <  1:950877/61‑4 (MQ=37)
                                                                                                     ccAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCCGGCCAGGGCGCTggccagg                                                                                                                                         >  2:950877/1‑58 (MQ=25)
                                                                                                          gcCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGGGGCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAgg                                           >  2:1113443/1‑150 (MQ=255)
                                                                                                            cTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGgc                                                                            <  1:1126105/115‑1 (MQ=255)
                                                                                                            cTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGgc                                                                            >  2:1126105/1‑115 (MQ=255)
                                                                                                            cTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCACGGGc                                         <  1:482132/150‑1 (MQ=255)
                                                                                                             tGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACt                                                   <  1:919129/139‑1 (MQ=255)
                                                                                                             tGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACt                                                   >  2:919129/1‑139 (MQ=255)
                                                                                                             tGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGcc                                        <  2:885558/150‑1 (MQ=255)
                                                                                                                    cAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGt                                 <  1:1127397/150‑1 (MQ=255)
                                                                                                                        tACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCGCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTgccg                             >  1:11310/1‑150 (MQ=255)
                                                                                                                          ccATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTgccgcg                           <  1:595376/150‑1 (MQ=255)
                                                                                                                            aTGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTgccgcgcc                         >  1:1228207/1‑150 (MQ=255)
                                                                                                                             tGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCAGCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTgccgcgccg                        >  2:953082/1‑150 (MQ=255)
                                                                                                                               gCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGtt                                                                                                                  >  1:677428/1‑58 (MQ=255)
                                                                                                                               gCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGtt                                                                                                                  <  2:677428/58‑1 (MQ=255)
                                                                                                                                cTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCt                     >  2:272866/1‑150 (MQ=255)
                                                                                                                                cTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGACGCCTGTCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCt                     >  1:114699/1‑150 (MQ=255)
                                                                                                                                       tGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTg              <  2:996540/150‑1 (MQ=255)
                                                                                                                                       tGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTCGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGCGCCGCGCCGCCTGCCCCTg              <  2:459308/150‑1 (MQ=255)
                                                                                                                                         cACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGa            <  1:953082/150‑1 (MQ=255)
                                                                                                                                          caggcggCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGAc           >  1:1265170/6‑150 (MQ=255)
                                                                                                                                           cccTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACt          >  1:263315/1‑150 (MQ=255)
                                                                                                                                            ccTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTgg                                                 <  1:382260/110‑1 (MQ=255)
                                                                                                                                            ccTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTgg                                                 >  2:382260/1‑110 (MQ=255)
                                                                                                                                             cTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTc        <  2:820872/150‑1 (MQ=255)
                                                                                                                                              tGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCt       <  2:491042/150‑1 (MQ=255)
                                                                                                                                              tGGCCAGCGCCCTGGCCGGCCACCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGc                                                                                             >  2:690425/1‑64 (MQ=255)
                                                                                                                                              tGGCCAGCGCCCTGGCCGGCCACCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGc                                                                                             <  1:690425/64‑1 (MQ=255)
                                                                                                                                               ggCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTg      <  1:1163204/150‑1 (MQ=255)
                                                                                                                                                gCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTGc     <  2:1265170/150‑1 (MQ=255)
                                                                                                                                                 ccAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGAAGACGGGATAAGGCGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTGCt    >  2:1106801/1‑150 (MQ=255)
                                                                                                                                                  cAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTGCTc   >  2:420886/1‑150 (MQ=255)
                                                                                                                                                   aGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTGCTcc  >  1:470857/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCGTTCTGTAGACACCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTGCTCC  >  NC_002947/5445591‑5445887

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