Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,015,333 A→T 5.5% intergenic (+19/+21) metE → / ← ysgA 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase/putative carboxymethylenebutenolidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,015,3330AT 5.5% 310.9 / 7.8 91intergenic (+19/+21)metE/ysgA5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase/putative carboxymethylenebutenolidase
Reads supporting (aligned to +/- strand):  ref base A (37/49);  new base T (3/2);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

GGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTAC  >  NC_000913/4015185‑4015480
                                                                                                                                                    |                                                                                                                                                   
ggCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGGACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTaa                                                                                                                                                     <  3:167632/149‑1 (MQ=255)
ggCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTaa                                                                                                                                                     <  1:497900/149‑1 (MQ=255)
ggCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTaa                                                                                                                                                     <  3:136039/149‑1 (MQ=255)
  cAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATa                                                                                                                                                   >  3:380085/1‑149 (MQ=255)
   aaaaCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATac                                                                                                                                                  >  2:189957/1‑149 (MQ=255)
   aaaaCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATac                                                                                                                                                  >  2:419968/1‑149 (MQ=255)
   aaaaCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATac                                                                                                                                                  >  1:417542/1‑149 (MQ=255)
        gCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCAccc                                                                                                                                             >  1:644932/1‑149 (MQ=255)
         cATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCg                                                                                                                                            <  2:803063/149‑1 (MQ=255)
           ttCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCGCCCGGt                                                                                                                                          <  1:691201/149‑1 (MQ=255)
             ccGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGttt                                                                                                                                        <  1:575054/149‑3 (MQ=255)
                gCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATa                                                                                                                                                   <  1:777450/135‑1 (MQ=255)
                gCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATa                                                                                                                                                   >  2:777450/1‑135 (MQ=255)
                   gagCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtctc                                                                                                                                  <  3:161306/149‑1 (MQ=255)
                   gagCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtctc                                                                                                                                  <  4:225705/149‑1 (MQ=255)
                   gagCGCCTGTGGGTCAACCCGGACTGCGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAATATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAAACCAATCCGGGTGGTAATGCCACCCGGTCTTTtctc                                                                                                                                  <  4:535815/149‑1 (MQ=255)
                    agCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCa                                                                                                                                 <  3:751783/149‑1 (MQ=255)
                         cTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACa                                                                                                                            <  4:446002/149‑1 (MQ=255)
                         cTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACa                                                                                                                            >  2:98502/1‑149 (MQ=255)
                          tgtgGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAg                                                                                                                           >  2:146619/1‑149 (MQ=255)
                            tgGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCg                                                                                                                         <  3:109862/149‑1 (MQ=255)
                             gggTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtctc                                                                                                                                  >  2:815324/1‑139 (MQ=255)
                             gggTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtctc                                                                                                                                  <  1:815324/139‑1 (MQ=255)
                                tCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGActt                                                                                                                     >  2:342806/1‑149 (MQ=255)
                                 cAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGActtc                                                                                                                    <  1:98502/149‑1 (MQ=255)
                                  aaCCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCCCCCGGTCTTTTCTCATTACAGCGActtct                                                                                                                   >  1:867168/1‑149 (MQ=255)
                                   aCCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGActtctt                                                                                                                  <  2:485346/149‑1 (MQ=255)
                                     ccGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTcc                                                                                                                <  1:421774/149‑1 (MQ=255)
                                     ccGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTcc                                                                                                                >  3:82233/1‑149 (MQ=255)
                                            gtgGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATa                                                                                                         <  1:496148/149‑1 (MQ=255)
                                              ggCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACt                                                                                                       <  2:646032/149‑1 (MQ=255)
                                               gCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGTGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTg                                                                                                      >  2:230420/1‑149 (MQ=255)
                                               gCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTg                                                                                                      >  1:436756/1‑149 (MQ=255)
                                                  tGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCtt                                                                                                   >  4:762780/1‑149 (MQ=255)
                                                     aaaCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGTGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaaa                                                                                                <  4:605956/149‑1 (MQ=255)
                                                     aaaCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaaa                                                                                                <  3:180744/149‑1 (MQ=255)
                                                     aaaCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaaa                                                                                                <  2:561968/149‑1 (MQ=255)
                                                     aaaCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaaa                                                                                                <  1:378362/149‑1 (MQ=255)
                                                     aaaCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaaa                                                                                                <  1:142015/149‑1 (MQ=255)
                                                         gcgcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCccac                                                                                                             >  3:357435/1‑132 (MQ=255)
                                                         gcgcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCccac                                                                                                             <  4:357435/132‑1 (MQ=255)
                                                             ggCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCa                                                                                        <  1:386553/149‑1 (MQ=255)
                                                             ggCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCa                                                                                        >  3:843862/1‑149 (MQ=255)
                                                             ggCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCa                                                                                        <  4:208179/149‑1 (MQ=255)
                                                             ggCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCa                                                                                        >  1:454320/1‑149 (MQ=255)
                                                                   ccAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGtttgg                                                                                                                                      <  2:509386/97‑5 (MQ=255)
                                                                   ccAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGtttgg                                                                                                                                      >  1:509386/1‑93 (MQ=255)
                                                                      gAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCt                                                                               <  3:824555/149‑1 (MQ=255)
                                                                           ccGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccag                                                                          <  3:30838/149‑1 (MQ=255)
                                                                             gcgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCGCCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagcc                                                                        <  1:189957/149‑1 (MQ=255)
                                                                             gcgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagcc                                                                        <  4:298359/149‑1 (MQ=255)
                                                                             gcgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagcc                                                                        <  1:266482/149‑1 (MQ=255)
                                                                               gcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCGTTCCAGCATACGCTGCCAGCCAt                                                                      >  1:185737/1‑149 (MQ=255)
                                                                                 ggCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCGTTCCAGCATACGCTGCCAGCCAtct                                                                    <  2:185737/149‑1 (MQ=255)
                                                                                     cTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTg                                                                <  2:529752/149‑1 (MQ=255)
                                                                                     cTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTg                                                                <  1:320512/149‑1 (MQ=255)
                                                                                      tGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGc                                                               >  1:98675/1‑149 (MQ=255)
                                                                                      tGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGc                                                               >  1:3763/1‑149 (MQ=255)
                                                                                           aaCATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATg                                                          <  1:342806/149‑1 (MQ=255)
                                                                                            aCATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGGGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGc                                                         >  1:404339/1‑149 (MQ=255)
                                                                                              aTGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCgg                                                       <  1:146619/149‑1 (MQ=255)
                                                                                                  tGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATg                                                   >  1:312553/1‑149 (MQ=255)
                                                                                                   gCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGa                                                  >  4:320563/1‑149 (MQ=255)
                                                                                                    cAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGAt                                                 <  3:647623/149‑1 (MQ=255)
                                                                                                       gcggcgCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGc                                              <  4:709817/149‑1 (MQ=255)
                                                                                                         ggcgCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCAAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTc                                            >  3:690193/1‑149 (MQ=255)
                                                                                                          gcgcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCg                                           <  1:187092/149‑1 (MQ=255)
                                                                                                          gcgcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATACCTCg                                           <  2:287329/149‑1 (MQ=255)
                                                                                                              aGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCg                                       >  4:411721/1‑149 (MQ=255)
                                                                                                               gAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTcc                                                                                                                >  4:277190/1‑75 (MQ=255)
                                                                                                               gAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTcc                                                                                                                <  3:277190/75‑1 (MQ=255)
                                                                                                               gAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGa                                      <  3:126608/149‑1 (MQ=255)
                                                                                                                aaCTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGAt                                     <  2:98675/149‑1 (MQ=255)
                                                                                                                aaCTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGAt                                     <  4:380085/149‑1 (MQ=255)
                                                                                                                aaCTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGAt                                     <  2:3763/149‑1 (MQ=255)
                                                                                                                  cTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATaa                                   >  1:763444/1‑149 (MQ=255)
                                                                                                                       gTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCg                                                        <  3:190298/123‑1 (MQ=255)
                                                                                                                       gTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCg                                                        >  4:190298/1‑123 (MQ=255)
                                                                                                                          gggggTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGtt                           >  1:423700/1‑149 (MQ=255)
                                                                                                                                aaaaTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatgc                     <  1:419968/149‑1 (MQ=255)
                                                                                                                                aaaaTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatgc                     <  3:411721/149‑1 (MQ=255)
                                                                                                                                aaaaTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatgc                     <  2:644932/149‑1 (MQ=255)
                                                                                                                                   aTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatgcatg                  >  1:641734/1‑149 (MQ=255)
                                                                                                                                      cAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGccc               >  1:398630/1‑149 (MQ=255)
                                                                                                                                      aagaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGccc               >  2:647046/4‑149 (MQ=255)
                                                                                                                                       aaaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGAt                                                 <  4:136660/114‑1 (MQ=255)
                                                                                                                                       aaaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGAt                                                 >  3:136660/1‑114 (MQ=255)
                                                                                                                                       aaaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCg              >  3:360091/1‑149 (MQ=255)
                                                                                                                                        gaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCgg             >  4:429715/2‑149 (MQ=255)
                                                                                                                                        aaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCgg             >  2:617343/1‑149 (MQ=255)
                                                                                                                                              gtggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCg       >  1:100842/1‑149 (MQ=255)
                                                                                                                                              gtggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCg       >  2:439319/1‑149 (MQ=255)
                                                                                                                                               tggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCgg      <  2:398630/149‑1 (MQ=255)
                                                                                                                                                  tAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATCGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTa   <  3:265784/149‑1 (MQ=255)
                                                                                                                                                   aaTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTac  >  4:443781/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTAC  >  NC_000913/4015185‑4015480

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: