Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,445,740 C→G 6.4% 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 6.4% 199.3 / 4.7 63S216R (AGC→AGGthiDphosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base C (32/27);  new base G (2/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

CCGTTCTGTAGACACCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTGCTC  >  NC_002947/5445591‑5445886
                                                                                                                                                     |                                                                                                                                                  
ccGTTCTGTAGACACCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAgc                                                                                                                                                    >  1:340975/1‑150 (MQ=255)
 cGTTCTGTAGACACCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAgcg                                                                                                                                                   >  1:241442/1‑150 (MQ=255)
     cTGTAGACACCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCt                                                                                                                                               <  2:340975/150‑1 (MQ=255)
        tAGACGCCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCGCGGGCTGCACCCTGGCCAGCGCCCTggc                                                                                                                                            >  1:322416/1‑150 (MQ=255)
          gACACCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTggccg                                                                                                                                          <  2:109705/150‑1 (MQ=255)
               ctgctgATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGgccgc                                                                                                                                     >  2:44826/1‑150 (MQ=255)
                tgctgATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCTAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGgccgcc                                                                                                                                    >  2:235925/1‑150 (MQ=255)
                      ttACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGgcgc                                                                                                                              <  1:44826/150‑1 (MQ=255)
                          cggcggTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCCGGCCAGGGCGCTggccagggt                                                                                                                                      >  1:203952/1‑133 (MQ=255)
                          cggcggTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCCGGCCAGGGCGCTggccagggt                                                                                                                                      <  2:203952/138‑6 (MQ=255)
                                  aCGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGGGCGCTggccagg                                                                                                                                        <  2:651933/128‑4 (MQ=255)
                                  aCGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGGGCGCTggccagg                                                                                                                                        >  1:651933/1‑125 (MQ=255)
                                   cGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGtt                                                                                                                 <  1:417308/150‑1 (MQ=255)
                                   cGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGtt                                                                                                                 <  2:57202/150‑1 (MQ=255)
                                    ggCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCCTGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTg                                                                                                                <  2:477978/150‑1 (MQ=255)
                                          gaagacgaaATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAgc                                                                                                          >  1:657577/1‑150 (MQ=255)
                                          gaagacgaaATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAgc                                                                                                          >  2:136080/1‑150 (MQ=255)
                                              acgaaATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGTGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCg                                                                                                      >  1:160094/1‑150 (MQ=255)
                                                  aaTCCACAACCGCCTGTACACCCAGGACGGCCAGAACAACATTTGAAATTGTTTCTGCCTGTCTGGTAGGCACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTgcg                                                                                                  >  1:21762/1‑150 (MQ=255)
                                                      cacaACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCAcc                                                                                              >  2:611958/1‑150 (MQ=255)
                                                           ccGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCt                                                                                         <  2:128656/150‑1 (MQ=255)
                                                             gCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTAGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTgg                                                                                       >  1:568971/1‑150 (MQ=255)
                                                              ccTGTACACCCAGGACGGCCAGCACTTCACTTGGAGTTGCCAGCGCCTGCCGGGCACCTACCATGCCTCCGCTTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGGGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGa                                                                                      <  1:132033/150‑1 (MQ=255)
                                                              ccTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGa                                                                                      <  2:644132/150‑1 (MQ=255)
                                                               cTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGAc                                                                                     <  1:611958/150‑1 (MQ=255)
                                                                       ccAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCCGCGCCCTGGCCGGTCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTg                                                                             >  1:91061/1‑150 (MQ=255)
                                                                       ccAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTg                                                                             >  2:992830/1‑150 (MQ=255)
                                                                           gACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGgcgc                                                                         >  1:215885/1‑150 (MQ=255)
                                                                             cGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCAc                                                                       <  1:161698/150‑1 (MQ=255)
                                                                                       ttcacttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGcc                                                                                                       >  2:1043998/1‑108 (MQ=255)
                                                                                       ttcacttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGcc                                                                                                       >  1:547278/1‑108 (MQ=255)
                                                                                       ttcacttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGcc                                                                                                       <  2:547278/108‑1 (MQ=255)
                                                                                       ttcacttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGcc                                                                                                       <  1:1043998/108‑1 (MQ=255)
                                                                                            ttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGGCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGa                                                        <  2:363125/150‑1 (MQ=255)
                                                                                            ttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCCGTTTTAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGa                                                        <  1:136080/150‑1 (MQ=255)
                                                                                            ttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGa                                                        <  2:163419/150‑1 (MQ=255)
                                                                                            ttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGa                                                        >  1:1079361/1‑150 (MQ=255)
                                                                                            ttGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGa                                                        >  1:743492/1‑150 (MQ=255)
                                                                                                    gCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTgg                                                <  2:241442/150‑1 (MQ=255)
                                                                                                     ccAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGcc                                                                                           <  2:751467/106‑1 (MQ=255)
                                                                                                     ccAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGcc                                                                                           >  1:751467/1‑106 (MQ=255)
                                                                                                     ccAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTggg                                               <  1:992830/150‑1 (MQ=255)
                                                                                                      cAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGa                                                        >  2:52854/1‑140 (MQ=255)
                                                                                                      cAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGa                                                        <  1:52854/140‑1 (MQ=255)
                                                                                                      cAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCGCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGc                                              >  2:972560/1‑150 (MQ=255)
                                                                                                          gcCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCTCCCTGCGTGACGCCGAGCAACTGGGCAAgg                                          <  1:362008/150‑1 (MQ=255)
                                                                                                          gcCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAgg                                          <  1:235925/150‑1 (MQ=255)
                                                                                                                  ggCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATAc                                  >  1:128865/1‑150 (MQ=255)
                                                                                                                       cTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCTCGCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCAGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTgcc                             >  1:868963/1‑150 (MQ=255)
                                                                                                                         aCCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCGGGACTACACCTGGCGCACACTGCGCGACGCCAAGCAACTGGGCAAGAGCCAATACGTgacgc                           >  2:540028/1‑150 (MQ=255)
                                                                                                                           cATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTgccgcgc                         >  1:327382/1‑150 (MQ=255)
                                                                                                                                 tCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGAACTCGGGCAATACGTGCCGCGCCGCCTg                   >  1:840148/1‑150 (MQ=255)
                                                                                                                                       tGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTg             <  1:626723/150‑1 (MQ=255)
                                                                                                                                          ccccTGGCCAGCACCCTTGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGAc          <  2:868963/149‑1 (MQ=255)
                                                                                                                                          aCCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGAc          >  2:503606/1‑150 (MQ=255)
                                                                                                                                          aCCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGAGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGAc          >  2:218260/1‑150 (MQ=255)
                                                                                                                                          aCCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCACTATTAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAAGTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGAc          <  1:397067/150‑1 (MQ=255)
                                                                                                                                           cccTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAgg                                          >  1:145955/1‑117 (MQ=255)
                                                                                                                                           cccTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAgg                                          <  2:145955/117‑1 (MQ=255)
                                                                                                                                           cccTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCCGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACt         >  2:1012903/1‑150 (MQ=255)
                                                                                                                                            ccTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACtt        >  1:563542/1‑150 (MQ=255)
                                                                                                                                             cTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGccc                >  1:399585/1‑141 (MQ=255)
                                                                                                                                             cTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGccc                <  2:399585/141‑1 (MQ=255)
                                                                                                                                               ggCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTg     >  2:248388/1‑150 (MQ=255)
                                                                                                                                                  cAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTGCTc  >  2:45255/1‑150 (MQ=255)
                                                                                                                                                    gcgcCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGa                                                                                      <  1:693227/64‑1 (MQ=255)
                                                                                                                                                    gcgcCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGa                                                                                      >  2:693227/1‑64 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCGTTCTGTAGACACCTGCTGATTACCGGCGGTCACGGCGACGAAGACGAAATCCACAACCGCCTGTACACCCAGGACGGCCAGCACTTCACTTGGACTTGCCAGCGCCTGCCGGGCAGCTACCATGGCTCGGGCTGCACCCTGGCCAGCGCCCTGGCCGGCCGCCTGGCGCTCGGCGAGCAGTTGGAAAGCGCCGTGCGCACCGCCCTGGACTACACCTGGCGCACCCTGCGTGACGCCGAGCAACTGGGCAAGGGCCAATACGTGCCGCGCCGCCTGCCCCTGGACTTCTGCTC  >  NC_002947/5445591‑5445886

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: