Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,229,433:1 +G 7.0% intergenic (+71/+20) yjbB → / ← pepE putative Na+/Pi‑cotransporter/peptidase E, alpha‑aspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,229,4331.G 7.0% 574.6 / 37.0 143intergenic (+71/+20)yjbB/pepEputative Na+/Pi‑cotransporter/peptidase E, alpha‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base . (66/67);  new base G (4/6);  total (70/73)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

Only 100 of 145 total aligned reads displayed.

CCTGAACTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCA  >  NC_000913/4229284‑4229579
                                                                                                                                                      |                                                                                                                                                  
ccTGAACTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAt                                                                                                                                                     >  1:583879/1‑150 (MQ=255)
ccTGAACTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAt                                                                                                                                                     >  2:693924/1‑150 (MQ=255)
  tGAACTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCa                                                                                                                                                  <  2:248777/150‑1 (MQ=255)
       tCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCAt                                                                                                                                             <  2:1134971/150‑1 (MQ=255)
       tCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCAt                                                                                                                                             <  2:1730452/150‑1 (MQ=255)
          cTGTTTTGTACGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCTTCGCATCAGGCACTATGAAAAAGATGCCGGATCAATCATCCg                                                                                                                                          >  2:594676/1‑150 (MQ=39)
            gTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGc                                                                                                                                        <  1:1176351/150‑1 (MQ=255)
            gTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGc                                                                                                                                        >  2:1504143/1‑150 (MQ=255)
            gTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGc                                                                                                                                        <  1:1411967/150‑1 (MQ=255)
            gTTTTGTTCGGTGGCTTACAGGGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGc                                                                                                                                        <  2:91779/150‑1 (MQ=255)
              tttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGCat                                                                                                                                      <  1:1582794/150‑2 (MQ=255)
              tttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGCat                                                                                                                                      <  2:1255211/150‑2 (MQ=255)
              tttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCaa                                                                                                                                      >  2:131809/1‑150 (MQ=255)
               ttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAAt                                                                                                                                     <  1:1398145/150‑1 (MQ=255)
                tGTTCGGTGGCTTACAGTGTTCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTACGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGCatct                                                                                                                                    <  1:1320554/150‑4 (MQ=38)
                tGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGAATCAGGCACAATGAACAAGATGCCGGATCATGCATCCGGCAAtt                                                                                                                                    >  1:759824/1‑150 (MQ=255)
                 gTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTc                                                                                                                                   <  2:1250227/150‑1 (MQ=255)
                   tCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCAc                                                                                                                                 <  1:529270/150‑1 (MQ=255)
                    cGGTGGCTTTCAGTTTGCTGGAACAGCCGGATGAAGATGAAGGACGCGACGAGTATAAATATTGTGAATCCTGTAGGCCTGATACGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACt                                                                                                                                <  1:1074585/150‑1 (MQ=39)
                    cGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACt                                                                                                                                <  2:378505/150‑1 (MQ=255)
                     ggTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACtt                                                                                                                               >  1:1526151/1‑150 (MQ=255)
                     ggTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACtt                                                                                                                               >  2:853723/1‑150 (MQ=255)
                      gTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTa                                                                                                                              >  1:277731/1‑150 (MQ=255)
                      gTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTa                                                                                                                              >  2:1633843/1‑150 (MQ=255)
                       tGGCTTACAGTGTTCTGGAACAGCCGGATGAATATGTAGGACGGGACGAGTATTAATTTTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTaa                                                                                                                             <  2:150557/150‑1 (MQ=255)
                           ttACAGTGTGCTGGAACAGCCGGATGAAGATGAATGACGTGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAAGAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAAc                                                                                                                         <  1:249537/150‑1 (MQ=255)
                               aGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTg                                                                                                                     <  1:1022796/150‑1 (MQ=255)
                                 tgtgCTGGAACAGCCGGATGAAGATGAAGGACGGGACTAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGAc                                                                                                                   <  1:1108737/150‑1 (MQ=255)
                                 tgtgCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGAc                                                                                                                   <  1:1517779/150‑1 (MQ=255)
                                  gtgCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGAcc                                                                                                                  <  1:754815/150‑1 (MQ=255)
                                         aaCAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTcc                                                                                                           <  1:1545593/150‑1 (MQ=255)
                                           cAGCCGGAAGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGCCCAGCTTCCAg                                                                                                         >  1:802534/1‑150 (MQ=255)
                                            aGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGc                                                                                                        <  1:1728300/150‑1 (MQ=255)
                                             gCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCg                                                                                                       >  1:1171469/1‑150 (MQ=255)
                                               cGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGa                                                                                                     <  1:988609/150‑1 (MQ=255)
                                                ggATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGaa                                                                                                    >  1:820212/1‑150 (MQ=255)
                                                   tgaagctgaagGACGGGACGAGTATTAATCTTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGTATCATGCATCCGGCAATTCACTTAAAAACTGTGACCAGCTTCCAGCGGAACCg                                                                                                 <  2:216427/150‑1 (MQ=255)
                                                      agatgaagGACGGGACGAGTATTAATATTGTGATTCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTCACCAGCTTCCAGCGGAACCGctt                                                                                              >  2:1196256/1‑150 (MQ=255)
                                                      agatgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGctt                                                                                              >  1:1035429/1‑150 (MQ=255)
                                                        atgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttct                                                                                            <  2:1410476/150‑1 (MQ=255)
                                                          gaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCttcttc                                                                                          <  2:811509/150‑1 (MQ=255)
                                                            agGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCAGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCAc                                                                                        >  2:921133/1‑150 (MQ=255)
                                                                cgggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAGCGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGc                                                                                    >  1:1252041/1‑150 (MQ=255)
                                                                cgggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGc                                                                                    >  2:1406877/1‑150 (MQ=255)
                                                                  ggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGTCACTATGAACAAGATGCCGGATCATGCATCCGTCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCTCCAGCCt                                                                                  >  2:1370987/1‑150 (MQ=255)
                                                                       aGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAac                                                                             <  2:737551/150‑1 (MQ=255)
                                                                       aGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAac                                                                             <  1:961305/150‑1 (MQ=255)
                                                                         tATTAATATTGTGAATCCTGTATGCCTGATAAGCCGCGGCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAacac                                                                           <  1:166498/150‑1 (MQ=255)
                                                                              atatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAg                                                                      >  1:86919/1‑150 (MQ=255)
                                                                                atTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtg                                                                    >  1:1690878/1‑150 (MQ=255)
                                                                                  tgtgAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggt                                                                  <  1:1594887/150‑1 (MQ=255)
                                                                                    tgAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCAAGCATCCGGCAATTAACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGgtgt                                                                >  1:440442/1‑150 (MQ=255)
                                                                                         ccTGTATGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCggg                                                           <  1:460700/150‑1 (MQ=255)
                                                                                         ccTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCAGTTAAAAACGGTGACCAGCTTCCAGCGGCCCCGCTTCTTCCCCAGCCATAAACACATAAGTGGTGTTCggg                                                           >  1:1091662/1‑150 (MQ=38)
                                                                                         ccTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCggg                                                           <  1:1283009/150‑1 (MQ=255)
                                                                                            gTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCa                                                        >  1:638964/1‑150 (MQ=255)
                                                                                               ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCAccc                                                     <  1:1702826/150‑1 (MQ=255)
                                                                                                gCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGACCCGCTTCTTCACCAGCCTTAAACCCATAAGTGGTGTTCGGTCCACCCa                                                    >  2:1094204/1‑150 (MQ=255)
                                                                                                 ccTGATGCGACGCTTGACGCGTCTTATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCAACAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCag                                                   >  1:877975/1‑150 (MQ=255)
                                                                                                   tGATAAGACGCGTCAAGCGTCGCATCAGGCACTTTGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcg                                                 >  2:911665/1‑150 (MQ=255)
                                                                                                    gATAATACGCGTCAAGAGTCGCATCAGGAACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCAATTAAAAACGGTGACCAGCTTCAATCGGAACCCCTAGATCACCAGCCTTAAACCCATAAGTGGTGTTCGGGCCCCCCcgcgt                                                >  2:157090/1‑150 (MQ=255)
                                                                                                     aTAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgta                                               >  2:544031/1‑150 (MQ=255)
                                                                                                       taGACGCGTCAAGCGTCGCATCAGGCACTAAGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagt                                             <  1:594676/149‑2 (MQ=255)
                                                                                                       aaGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGTATCATGCATACGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagc                                             <  2:1675747/150‑1 (MQ=255)
                                                                                                       aaGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagc                                             >  2:768222/1‑150 (MQ=255)
                                                                                                             cgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAcc                                       <  1:154710/150‑1 (MQ=255)
                                                                                                             cgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCGTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAcc                                       <  2:1051668/150‑1 (MQ=255)
                                                                                                               tCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCtt                                     >  1:695824/1‑150 (MQ=255)
                                                                                                               tCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCtt                                     <  1:500356/150‑1 (MQ=255)
                                                                                                               tCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCtt                                     >  1:382157/1‑150 (MQ=255)
                                                                                                               tCAAGCGTCGCATCAGGCACTATGAACAAGATGCAGGATCATGCATCCGGCAATTCACTTAAAACCTGTGACCAGCTTCCCGCGGAACCGATTCTACAACAGCCTTAAACACATAACAGGTTAACGGGCCACCCAGCGAAGCGTCTCCtt                                     >  1:75558/1‑150 (MQ=255)
                                                                                                                   gCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACt                                 <  2:890704/150‑1 (MQ=255)
                                                                                                                      tCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                              <  2:1429208/150‑1 (MQ=255)
                                                                                                                       cGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAGCGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACt                             <  2:1252041/150‑1 (MQ=255)
                                                                                                                        gCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACtg                            <  1:1159848/150‑1 (MQ=255)
                                                                                                                          aTCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACtgtg                          >  1:1154226/1‑150 (MQ=255)
                                                                                                                             aGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATc                       <  2:438346/150‑1 (MQ=255)
                                                                                                                              ggCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAATGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATcc                      >  2:490706/1‑150 (MQ=255)
                                                                                                                              ggCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATcc                      <  2:312567/150‑1 (MQ=255)
                                                                                                                              ggCACAATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATcc                      <  2:759824/150‑1 (MQ=255)
                                                                                                                                cACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAg                    <  2:583879/150‑1 (MQ=255)
                                                                                                                                cACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACAGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAg                    <  1:131809/150‑1 (MQ=255)
                                                                                                                                 aCTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGt                   >  2:794349/1‑150 (MQ=255)
                                                                                                                                  cTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGtt                  >  2:358307/1‑150 (MQ=255)
                                                                                                                                   tATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATTAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTa                 >  2:193874/1‑150 (MQ=255)
                                                                                                                                     tGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGCCCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTAcc               >  1:1228598/1‑150 (MQ=255)
                                                                                                                                     tGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCCGCGGAACCGCTTCTTCCCCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCCCTGTGATCCAGTTAcc               >  1:1522415/1‑150 (MQ=255)
                                                                                                                                     tGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTAcc               <  1:1504143/150‑1 (MQ=255)
                                                                                                                                      gAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCt              >  1:183984/1‑150 (MQ=255)
                                                                                                                                      gAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCt              <  1:887198/150‑1 (MQ=255)
                                                                                                                                       aaCAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTTGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCtt             <  2:1356370/150‑1 (MQ=255)
                                                                                                                                       aaCAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGGGACCTTTACTCACGGTGATCCAGTTACCtg             >  1:884839/1‑149 (MQ=255)
                                                                                                                                           aGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGATTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAACTGGTGTTCGGGCCACCCAGCGGAGCGTGCCCTGTTCTCCCTGTGATCCAGTTACCTTCCgg         >  2:697467/1‑150 (MQ=255)
                                                                                                                                            gATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGt        <  1:370192/150‑1 (MQ=255)
                                                                                                                                            gATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGt        >  1:711956/1‑150 (MQ=255)
                                                                                                                                              tGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAg      >  1:1612055/1‑150 (MQ=255)
                                                                                                                                                ccGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGAc    <  1:1633843/150‑1 (MQ=255)
                                                                                                                                                 cGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCATAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGAcc   <  2:845671/150‑1 (MQ=255)
                                                                                                                                                  ggATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCa  >  1:150277/1‑150 (MQ=255)
                                                                                                                                                  ggATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCa  >  2:513662/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CCTGAACTCGCTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCA  >  NC_000913/4229284‑4229579

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