Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,285,234 G→C 5.3% E287Q (GAG→CAG)  yejM → essential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,285,2340GC 5.3% 281.4 / 3.4 88E287Q (GAG→CAG) yejMessential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor
Reads supporting (aligned to +/- strand):  ref base G (51/32);  new base C (2/3);  total (53/35)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCGTACGCCTGCGGCATTAATTACTGCGCTTAATC  >  NC_000913/2285087‑2285381
                                                                                                                                                   |                                                                                                                                                   
aaTCCAGAGTCCGTTTCCGTTCTGTATCCGTTACGCGTACTGCGCTATCGCGATATTGGCACCGGGCAGAATGTGCTGTTGATTCCTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAg                                                                                                                                                   <  1:1059603/150‑1 (MQ=255)
aaTCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAg                                                                                                                                                   <  1:898746/150‑1 (MQ=255)
aaTCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAg                                                                                                                                                   >  2:530937/1‑150 (MQ=255)
   ccAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGATCTGCGCTATCGCGATATGGGCACAGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCaa                                                                                                                                                >  2:537759/1‑150 (MQ=255)
        cGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTCAGCaaaacct                                                                                                                                           <  2:471659/150‑4 (MQ=255)
         gCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTCAGCaaaacctg                                                                                                                                          <  2:356347/150‑5 (MQ=255)
         gCCGATTCCGTTCAGTATCCGTTACTCGAACTGCGCTATCGGGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCTATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTCAGCaaaacctg                                                                                                                                          <  2:238110/150‑5 (MQ=255)
            gTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCg                                                                                                                                       <  1:828266/150‑1 (MQ=255)
                ccGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCa                                                                                                                                   <  2:663711/150‑1 (MQ=255)
                ccGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCa                                                                                                                                   >  2:500069/1‑150 (MQ=255)
                ccGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCa                                                                                                                                   <  2:239871/150‑1 (MQ=255)
                  gTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCAcg                                                                                                                                 >  1:1270704/1‑150 (MQ=255)
                  gTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGAACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTAGCTGAGCAAAATATTTCGTTCAcg                                                                                                                                 >  1:1270705/1‑150 (MQ=255)
                         aTCCGTTAAGCGAACTGCGCTATCGCGATATGGGCCCCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTCAGCAAAATATTTCGTTCACGCGCCata                                                                                                                          >  2:232558/1‑150 (MQ=255)
                         aTCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGCTGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCata                                                                                                                          <  1:511091/150‑1 (MQ=255)
                         aTCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCaga                                                                                                                          >  2:275075/1‑148 (MQ=255)
                            cGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGa                                                                                                                       >  2:1088182/1‑150 (MQ=255)
                                aaGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGAGGGCCTGAACTACTCACGGTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTTGTTAACGCGCCTTATGAGCTc                                                                                                                   >  2:1159067/1‑150 (MQ=255)
                                  gCGAACCGCGCTATCGCGTTATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTCCTAACGCTTCGAGAAGCAGATGACTGCGCTGGCAGGTTTTGCTGAGCAAAATATATCGTTCACGCGCCATATCAGCACCg                                                                                                                 >  1:339174/1‑150 (MQ=255)
                                   cGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCgg                                                                                                                >  2:293692/1‑150 (MQ=255)
                                    gAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCATCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGc                                                                                                               >  2:696918/1‑150 (MQ=255)
                                     aaCTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCa                                                                                                              <  1:1169746/150‑1 (MQ=255)
                                            cTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCAGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTATGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTAc                                                                                                       >  2:466086/1‑150 (MQ=255)
                                             tATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACa                                                                                                      <  1:1257289/150‑1 (MQ=255)
                                              aTCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCCCGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCACCACTACAg                                                                                                     >  1:1501132/1‑150 (MQ=255)
                                                ctcgATATGGGCACCGTGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCTTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGAc                                                                                                   <  1:1324347/148‑1 (MQ=255)
                                                cgcgATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCCCGCGCCATATGAGCTCCGGCACCACTACAgcc                                                                                                   >  2:1497941/1‑148 (MQ=255)
                                                  cgATATGGGCACCGGGCAGAATGTGCTGTTGATTAATGTCGATGGCCTGAAATAATCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGTTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGCCAACACTACAGACaa                                                                                                 >  1:1237923/1‑150 (MQ=255)
                                                        gggCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCatat                                                                                                                         >  1:1632425/1‑120 (MQ=255)
                                                        gggCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCatat                                                                                                                         <  2:1632425/120‑1 (MQ=255)
                                                             ccGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTg                                                                                      <  1:1108294/150‑1 (MQ=255)
                                                               gggCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGc                                                                                    >  1:1496281/1‑150 (MQ=255)
                                                                ggCAGATTGTGCTTTTGATTACTGTCGATGGCCTGAAGTCCTCACGCTTCGAGAAGCTGATGCCTGCGCTGGCAGGTTATGCTGAGCAAAATATTTCGATCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGcc                                                                                   <  1:1296302/150‑1 (MQ=255)
                                                                ggCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGcc                                                                                   >  1:906001/1‑150 (MQ=255)
                                                                ggCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGcc                                                                                   <  2:1391280/150‑1 (MQ=255)
                                                                ggCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGcc                                                                                   >  1:1693922/1‑150 (MQ=255)
                                                                     aaTGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTc                                                                              <  2:766662/150‑1 (MQ=255)
                                                                         tgCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATg                                                                          <  1:683950/150‑1 (MQ=255)
                                                                            tGTTGATTACTGTAGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCa                                                                       >  1:871630/1‑150 (MQ=255)
                                                                             gTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCAt                                                                      >  1:894260/1‑150 (MQ=255)
                                                                                    aCTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCg                                                               >  1:733603/1‑150 (MQ=255)
                                                                                    aCTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGAGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTCGCATCTTCTAAGGCAACTCGccc                                                               >  1:160448/1‑149 (MQ=255)
                                                                                      tGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAg                                                             >  1:1730147/1‑150 (MQ=255)
                                                                                            tGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTatat                                                       >  1:1194440/1‑150 (MQ=255)
                                                                                             ggCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATg                                                      >  2:1663875/1‑150 (MQ=255)
                                                                                             ggCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATg                                                      <  2:1283403/150‑1 (MQ=255)
                                                                                             ggCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATg                                                      >  2:1121506/1‑150 (MQ=255)
                                                                                              gCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATgg                                                     <  1:535899/150‑1 (MQ=255)
                                                                                               ccTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTTTTCTATGGCATCTCGCCGAGCTATATGGa                                                    >  2:341405/1‑150 (MQ=255)
                                                                                                cTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGAc                                                   <  1:1669028/150‑1 (MQ=255)
                                                                                                 tGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACTACGGCATCTTTGGCATGTTCTATGGCATCTCGCCGAGCTATATGGACg                                                  >  1:883694/1‑150 (MQ=255)
                                                                                                  gAACTACTCACGCGGGGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTGAGGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGTCCTGTTCTATGGCTTCTCGCCGAGCTATATGGACgg                                                 <  1:798628/150‑1 (MQ=255)
                                                                                                    actactCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCa                                               >  1:1688877/1‑150 (MQ=255)
                                                                                                     ctactCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCAt                                              >  2:774058/1‑150 (MQ=255)
                                                                                                       actCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTc                                            >  1:909999/1‑150 (MQ=255)
                                                                                                       actCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTc                                            <  1:1343932/150‑1 (MQ=255)
                                                                                                       actCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTc                                            <  2:811365/150‑1 (MQ=255)
                                                                                                       actCACGCTTCGAGAAGCAGATGCCTGCGCTGGCACGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTCGCCTGTTCTATGGCATATCGCCGAGCTATATGGACGGCAtta                                            >  2:315066/1‑149 (MQ=255)
                                                                                                       actCACGCTTCGAGAAGCAGATGCCTGCGATGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCGATGGCATCTCGCCGAGCTATATGGACGGCAtta                                            >  1:111575/1‑149 (MQ=255)
                                                                                                        ctCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCt                                           >  1:877344/1‑150 (MQ=255)
                                                                                                          cACGCTTCGAGAAGCAGATGCCTGCGGTGGCACGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCACCACTACAGACAACGGCATCTTTGGCCTGTTATATGGCATCTCGCCGAGCTATATGGAAGTCATTCTGt                                         >  1:1040407/1‑150 (MQ=255)
                                                                                                          cACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGt                                         >  2:1591619/1‑150 (MQ=255)
                                                                                                           aCGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTc                                        <  2:692910/150‑1 (MQ=255)
                                                                                                              ctgcGAGAAGCAGCTTCCTGCGCTGGCAGGTGTTGCTGAGCAAAATATGTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCACCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGAc                                     <  1:702039/147‑1 (MQ=255)
                                                                                                               ttCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGAcc                                    <  1:1636335/150‑1 (MQ=255)
                                                                                                               ttCGAGAAGCAGAGGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCTTCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGAcc                                    <  2:700764/150‑1 (MQ=255)
                                                                                                                tCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGAccc                                   <  2:1184639/150‑1 (MQ=255)
                                                                                                                 cGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCg                                  >  2:976740/1‑150 (MQ=255)
                                                                                                                  gagaAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGTCGAGCTATATGGACGGCATTCTGTCGACCCGt                                 >  1:731011/1‑150 (MQ=255)
                                                                                                                  gagaAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCGt                                 >  2:1564180/1‑150 (MQ=255)
                                                                                                                   agaAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCGTa                                >  2:1087201/1‑150 (MQ=255)
                                                                                                                      aGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCGTACGc                             >  2:1486659/1‑150 (MQ=255)
                                                                                                                       gCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCGTACGcc                            <  1:638091/150‑1 (MQ=255)
                                                                                                                              ccTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCGTACGCCTGCGGCa                     >  1:705058/1‑150 (MQ=255)
                                                                                                                              ccTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCGTACGCCTGCGGCa                     >  1:191743/1‑150 (MQ=255)
                                                                                                                               cTGCGCTGGCAGGTTTTGCTGCGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCTTCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGCCGGCATTCTGTAGACCCGTACGCCTGCGGCat                    >  2:1191327/1‑150 (MQ=255)
                                                                                                                                tGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCGTACGCCTGCGGCatt                   <  2:894260/150‑1 (MQ=255)
                                                                                                                                  cgcTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGCCGGCATTCTGTCGACCCGTACGCCTGCGGCattaa                 >  1:1257058/1‑150 (MQ=255)
                                                                                                                                      ggCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCGTACGCCTGCGGCATTAATTAc             >  2:217215/1‑150 (MQ=255)
                                                                                                                                      gaaatattttGCTCAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCGTACGCCTGCGGCATTAATTAc             >  2:1350243/7‑150 (MQ=255)
                                                                                                                                        cAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCGTACGCCTGCGGCATTAATTACTg           >  2:1362466/1‑150 (MQ=255)
                                                                                                                                        cAGGTTTCGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGCCCACACTACAGACACCGGCATCTATGGCCTGATCTCTGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCGTACGCCTGCGGCATTAATTACTg           <  2:894457/150‑1 (MQ=255)
                                                                                                                                          ggTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCGTACGCCTGCGGCATTAATTACTgcg         <  1:1663875/150‑1 (MQ=255)
                                                                                                                                           gTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGCCGGCATTCTGTCGACCCGTACGCCTGCGGCATTAATTACTgcgc        >  1:1254241/1‑150 (MQ=255)
                                                                                                                                           gTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGCCCCGTACGGCTGCGGCATTATTTACTgcgc        >  2:122273/1‑150 (MQ=255)
                                                                                                                                              ttGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACATCGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATAAGGACGGCATTCTGTCGACCCGTACGCATGCGGCATTAATTACTGCGCTTa     <  1:293692/150‑1 (MQ=255)
                                                                                                                                                 cTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCGTACGCCTGCGGCATTAATTACTGCGCTTAATc  >  1:296333/1‑150 (MQ=255)
                                                                                                                                                 cTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCGTACGCCTGCGGCATTAATTACTGCGCTTAATc  <  2:1194440/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AATCCAGACGCCGTTTCCGTTCAGTATCCGTTAAGCGAACTGCGCTATCGCGATATGGGCACCGGGCAGAATGTGCTGTTGATTACTGTCGATGGCCTGAACTACTCACGCTTCGAGAAGCAGATGCCTGCGCTGGCAGGTTTTGCTGAGCAAAATATTTCGTTCACGCGCCATATGAGCTCCGGCAACACTACAGACAACGGCATCTTTGGCCTGTTCTATGGCATCTCGCCGAGCTATATGGACGGCATTCTGTCGACCCGTACGCCTGCGGCATTAATTACTGCGCTTAATC  >  NC_000913/2285087‑2285381

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