Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 52,968 A→G 5.0% V150A (GTG→GCG)  pdxA ← 4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091352,9680AG 5.0% 237.8 / 2.5 80V150A (GTG→GCG) pdxA4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base A (38/37);  new base G (2/2);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCACGCTCGCGCCAG  >  NC_000913/52820‑53116
                                                                                                                                                    |                                                                                                                                                    
ggCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCAtcatca                                                                                                                                                      <  2:791933/149‑1 (MQ=255)
 gCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATcac                                                                                                                                                     >  2:534649/1‑149 (MQ=255)
  cAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGCAGTTCTTCGGTCGCCAGCATCATcacc                                                                                                                                                    >  1:259152/1‑149 (MQ=255)
        cAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCttt                                                                                                                                              >  1:1536050/1‑149 (MQ=255)
         aaaTTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCtatt                                                                                                                                             >  1:1149249/1‑149 (MQ=255)
          aaTTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCttttt                                                                                                                                            >  1:1274844/1‑149 (MQ=255)
                 tCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGc                                                                                                                                     <  2:859655/149‑1 (MQ=255)
                 tCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGc                                                                                                                                     >  1:1097060/1‑149 (MQ=255)
                      aaaTCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACg                                                                                                                                <  1:534649/149‑1 (MQ=255)
                        aTCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGct                                                                                                                              >  1:1575136/1‑149 (MQ=255)
                         tCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGctc                                                                                                                             <  2:316696/149‑1 (MQ=255)
                           gTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCAGGCGaaaaaggtg                                                                                                                           <  1:554500/149‑8 (MQ=255)
                           gTGATGCAAAATAGCAATCAATTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGTCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCtt                                                                                                                           <  2:267291/149‑1 (MQ=255)
                            tGATGCAAAATAGCAATCACTTCGAGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTc                                                                                                                          <  2:825531/149‑1 (MQ=255)
                               tGCAAAATAGCAATCACATCGTGCAAAAGTGCAGCGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTTCCAGCGCCACGCGAAGTTCTTCGGTCGCCATCATCATCACCACCTTTTTCTCCTGCGAACGCTCTTCGaa                                                                                                                       >  1:242958/1‑149 (MQ=255)
                                 cAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGaaaa                                                                                                                     >  2:83718/1‑149 (MQ=255)
                                  aaaaTAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGATCGCTCTTCGaaaaa                                                                                                                    >  2:303083/1‑149 (MQ=255)
                                     aTAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTc                                                                                                                 <  2:1097060/149‑1 (MQ=255)
                                      tAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCg                                                                                                                <  2:585954/149‑1 (MQ=255)
                                      tAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCg                                                                                                                >  2:615244/1‑149 (MQ=255)
                                       aGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCgg                                                                                                               <  1:964901/149‑1 (MQ=255)
                                         cAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTa                                                                                                             <  1:144225/149‑1 (MQ=255)
                                         cAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACCCGGTa                                                                                                             <  1:453517/149‑1 (MQ=255)
                                          aaTCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCGCCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTAt                                                                                                            >  1:336114/1‑149 (MQ=255)
                                          aaTCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCGCCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTAt                                                                                                            >  1:341854/1‑149 (MQ=255)
                                          aaTCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGGAAATGCGTCGTTGCCAGCGCCACTAGAAGTTCTTCGGTCGCCAGCATCTTAACCACCTTTTTCGCCTGAGAACGCTCTTCGAAAAACTCGGTAt                                                                                                            >  2:200154/1‑149 (MQ=255)
                                            tCACTTCGTGCAAAAGTGCAGGGCTGATAGCGACAGCGATATCGCGCATCGGTAAATGCGTCCTTGCCAGCTCCACGCGAAGTTCTACGGTCGCCAGCATCATCACCACCTATTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGa                                                                                                          >  1:36188/1‑149 (MQ=255)
                                               cTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCg                                                                                                       <  2:1575136/149‑1 (MQ=255)
                                               cTTCGTGCAAAAGTGCAGGGGTGATACCGTCTGCGATATCGCGCAGCGGTAAATGCTTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCCCCACCTTTTTCGCCTGCGACCGCTCGTCGAAAACCTCCGTATGCCCg                                                                                                       >  1:1047279/1‑149 (MQ=255)
                                                ttCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCgg                                                                                                      >  1:7425/1‑149 (MQ=255)
                                                    tGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTaaa                                                                                                  <  2:1583085/149‑1 (MQ=255)
                                                     gCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTaaaa                                                                                                 <  2:900444/149‑1 (MQ=255)
                                                     gCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTaaaa                                                                                                 <  2:92752/149‑1 (MQ=255)
                                                      cAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAg                                                                                                >  1:1399069/1‑149 (MQ=255)
                                                         aaGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTACGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGaa                                                                                             <  1:860531/149‑1 (MQ=255)
                                                              cAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCa                                                                                        <  1:714216/149‑1 (MQ=255)
                                                              cAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCa                                                                                        <  2:717430/149‑1 (MQ=255)
                                                               aGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAg                                                                                       >  1:1134376/1‑149 (MQ=255)
                                                                 gggTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGcg                                                                                     >  2:727077/1‑149 (MQ=255)
                                                                    tGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGcgtcg                                                                                  >  1:47613/1‑149 (MQ=255)
                                                                     gATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGcgtcgt                                                                                 <  2:223024/149‑1 (MQ=255)
                                                                     gATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGcgtcgt                                                                                 <  2:182193/149‑1 (MQ=255)
                                                                       tAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCCCCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGGATGACCGGTAAAAGGAATGCCAGCGTCGTta                                                                               >  1:135218/1‑149 (MQ=255)
                                                                        aGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTtaa                                                                              <  2:463843/149‑1 (MQ=255)
                                                                        aGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTtaa                                                                              >  1:1031176/1‑149 (MQ=255)
                                                                         gCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCTGCGTCGTtaat                                                                             >  2:566851/1‑149 (MQ=255)
                                                                          cGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGATCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTtaata                                                                            >  1:106058/1‑149 (MQ=255)
                                                                                 gATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCtt                                                                     <  2:1031176/149‑1 (MQ=255)
                                                                                   tatCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTa                                                                   <  2:761229/149‑1 (MQ=255)
                                                                                      cgcgCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGc                                                                <  1:784801/149‑1 (MQ=255)
                                                                                       gcgcAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCa                                                               >  1:1093857/1‑149 (MQ=255)
                                                                                       gcgcAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCa                                                               >  1:681900/1‑149 (MQ=255)
                                                                                        cgcAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCAc                                                              <  2:276503/149‑1 (MQ=255)
                                                                                             cGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGAc                                                         >  1:1146177/1‑149 (MQ=255)
                                                                                             cGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGAc                                                         >  1:1105308/1‑149 (MQ=255)
                                                                                             cGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGAc                                                         >  2:65649/1‑149 (MQ=255)
                                                                                             cGGTAAATGCGACGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGAc                                                         <  2:241614/149‑1 (MQ=255)
                                                                                               gTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCt                                                       >  1:1295066/1‑149 (MQ=255)
                                                                                               gTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTGAATAACGCCTTTATGCACCGGCCCt                                                       >  2:783119/1‑149 (MQ=255)
                                                                                                tAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCtg                                                      <  2:1274844/149‑1 (MQ=255)
                                                                                                  aaTGCGTCGTTGCCAGCGCCTAGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTAGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCtgtg                                                    <  2:293153/149‑1 (MQ=255)
                                                                                                  aaTGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTAGAAAAACTCGGTATGACAGGTAAAAGGAATGCCAGCGTCGTTAATAACCCCTTTATGCACCGGACCtgtg                                                    >  2:920860/1‑149 (MQ=255)
                                                                                                   aTGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATAACAACCTTTTTCGCCTGCGAACGCTCTTCCAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATTCACCGGACCTGTGa                                                   >  2:42744/1‑149 (MQ=255)
                                                                                                        tcgtTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAgc                                              <  2:1307428/149‑1 (MQ=255)
                                                                                                                  cgcCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAAtt                                    <  2:106058/149‑1 (MQ=255)
                                                                                                                     cACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGc                                 <  2:1536050/149‑1 (MQ=255)
                                                                                                                     cACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGc                                 >  1:844418/1‑149 (MQ=255)
                                                                                                                     cACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGc                                 >  2:483294/1‑149 (MQ=255)
                                                                                                                      aCGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGcc                                <  2:1203518/149‑1 (MQ=255)
                                                                                                                          gAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTc                            >  2:387/1‑149 (MQ=255)
                                                                                                                           aaGTTCTTCGGTCGCCAGCATCATCGCCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCa                           <  2:341854/149‑1 (MQ=255)
                                                                                                                           aaGTTCTTCGGTCGCCAGCATCATCGCCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCa                           <  2:336114/149‑1 (MQ=255)
                                                                                                                            aGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAg                          >  1:142887/1‑149 (MQ=255)
                                                                                                                               tcttcGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACa                       >  1:1360887/1‑149 (MQ=255)
                                                                                                                                cttcGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACaa                      >  2:512463/1‑149 (MQ=255)
                                                                                                                                   cGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCa                   <  2:1149249/149‑1 (MQ=255)
                                                                                                                                     gTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATc                 <  2:663188/149‑1 (MQ=255)
                                                                                                                                      tCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCg                >  2:1420015/1‑149 (MQ=255)
                                                                                                                                         ccAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGACAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCAc             >  2:1398790/1‑149 (MQ=255)
                                                                                                                                          cAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCACg            >  2:1045588/1‑149 (MQ=255)
                                                                                                                                           aGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCACGc           <  2:1105308/149‑1 (MQ=255)
                                                                                                                                           aGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCACGc           <  2:1146177/149‑1 (MQ=255)
                                                                                                                                             catcatCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCACGCTc         <  1:615244/149‑1 (MQ=255)
                                                                                                                                                  tcaCCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCCCGCTCGCGcc    >  1:1252589/1‑149 (MQ=255)
                                                                                                                                                  tcaCCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCACGCTCGCGcc    <  1:83718/149‑1 (MQ=255)
                                                                                                                                                  tcaCCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCACGCTCGCGcc    >  1:1355979/1‑149 (MQ=255)
                                                                                                                                                  tcaCCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGGTCAGACAACCATCGCACGCTCGCGcc    >  1:1350042/1‑149 (MQ=255)
                                                                                                                                                    accaccTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCACGCTCGCGCCAg  >  2:662012/1‑149 (MQ=255)
                                                                                                                                                    accaccTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCACGCTCGCGCCAg  <  2:1134376/149‑1 (MQ=255)
                                                                                                                                                    accaccTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCACGCTCGCGCCAg  >  2:1257252/1‑149 (MQ=255)
                                                                                                                                                    accaccTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCACGCTCGCGCCAg  >  2:1488174/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCAATACCAAATTTGGTCCGCAAATCGTGATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGCGCCACGCGAAGTTCTTCGGTCGCCAGCATCATCACCACCTTTTTCGCCTGCGAACGCTCTTCGAAAAACTCGGTATGACCGGTAAAAGGAATGCCAGCGTCGTTAATAACGCCTTTATGCACCGGACCTGTGATCAGCGCGGCAAATTCGCCGTTCAGACAACCATCGCACGCTCGCGCCAG  >  NC_000913/52820‑53116

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑