Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,229,438 Δ1 bp 6.6% intergenic (+76/+15) 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,4380C. 6.6% 405.5 / 21.6 121intergenic (+76/+15)yjbB/pepEputative Na+/Pi‑cotransporter/peptidase E, alpha‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (46/65);  new base . (2/6);  total (49/72)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

Only 100 of 122 total aligned reads displayed.

CTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCG  >  NC_000913/4229294‑4229586
                                                                                                                                                 |                                                                                                                                                    
cTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATcc                                                                                                                                                  <  2:817932/149‑1 (MQ=255)
 tGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGCAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCg                                                                                                                                                 <  3:727098/149‑1 (MQ=255)
   ttttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGc                                                                                                                                               <  1:727405/149‑1 (MQ=255)
   ttttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGc                                                                                                                                               >  2:770757/1‑149 (MQ=255)
   ttttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTAGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGc                                                                                                                                               >  2:870668/1‑149 (MQ=255)
   ttttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCCTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGc                                                                                                                                               <  2:766808/149‑1 (MQ=255)
    tttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGCa                                                                                                                                              <  4:402201/149‑1 (MQ=255)
     ttGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGCat                                                                                                                                             <  3:103551/149‑2 (MQ=255)
      tGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGCatc                                                                                                                                            <  4:468693/149‑3 (MQ=255)
          cGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCAc                                                                                                                                        <  4:327442/149‑1 (MQ=255)
                 ttACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTaaaaa                                                                                                                                 >  2:149089/1‑149 (MQ=255)
                     aGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGt                                                                                                                             <  2:76465/149‑1 (MQ=255)
                     aGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGt                                                                                                                             <  2:675039/149‑1 (MQ=255)
                     aGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGt                                                                                                                             <  4:389277/149‑1 (MQ=255)
                     aGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGt                                                                                                                             <  2:577296/149‑1 (MQ=255)
                      gtgtGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTg                                                                                                                            <  2:95151/149‑1 (MQ=255)
                              gAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCtt                                                                                                                    >  3:840558/1‑149 (MQ=255)
                                   gCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGc                                                                                                               <  2:636613/149‑1 (MQ=255)
                                        atgaagatgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAAc                                                                                                          <  1:653774/149‑1 (MQ=255)
                                             gatgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGctt                                                                                                     >  4:491915/1‑149 (MQ=255)
                                              atgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACg                                                                                                                               >  2:10185/1‑122 (MQ=255)
                                              atgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACg                                                                                                                               <  1:10185/122‑1 (MQ=255)
                                              atgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttc                                                                                                    <  1:184704/149‑1 (MQ=255)
                                              atgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttc                                                                                                    <  3:521026/149‑1 (MQ=255)
                                              atgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttc                                                                                                    <  1:148785/149‑1 (MQ=255)
                                               tgaagGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttct                                                                                                   >  1:859408/1‑149 (MQ=255)
                                                  agGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCa                                                                                                <  2:222413/149‑1 (MQ=255)
                                                     acgggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCa                                                                                             <  3:742174/149‑1 (MQ=255)
                                                      cgggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAg                                                                                            <  1:149089/149‑1 (MQ=255)
                                                        ggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGcc                                                                                          <  4:322985/149‑1 (MQ=255)
                                                        ggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGcc                                                                                          <  3:102662/149‑1 (MQ=255)
                                                             attATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                                     <  3:188350/147‑1 (MQ=255)
                                                             aGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                                     <  1:686210/149‑1 (MQ=255)
                                                             aGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                                     <  4:517189/149‑1 (MQ=255)
                                                                aTTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGTATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAacac                                                                                  >  4:176857/1‑149 (MQ=255)
                                                                aTTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAacac                                                                                  >  1:61604/1‑149 (MQ=255)
                                                                aTTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAacac                                                                                  >  3:195371/1‑149 (MQ=255)
                                                                 ttAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAAcaca                                                                                 <  2:236466/149‑1 (MQ=255)
                                                                    atatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCa                                                                                             >  4:134640/1‑134 (MQ=255)
                                                                    atatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATaa                                                                              >  3:634734/1‑149 (MQ=255)
                                                                     tatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGCatc                                                                                                                                            >  4:521141/1‑84 (MQ=21)
                                                                     tatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATGATCCGGCatc                                                                                                                                            <  3:521141/86‑3 (MQ=21)
                                                                      atTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgt                                                                            <  3:428438/149‑1 (MQ=255)
                                                                      atTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgt                                                                            <  3:454496/149‑1 (MQ=255)
                                                                       ttGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtg                                                                           <  2:450209/149‑1 (MQ=255)
                                                                       ttGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtg                                                                           >  2:895437/1‑149 (MQ=255)
                                                                       ttGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtg                                                                           >  1:198414/1‑149 (MQ=255)
                                                                          tgAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggtg                                                                        <  1:536680/149‑1 (MQ=255)
                                                                            aaTCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCgg                                                                                                             <  2:242131/110‑1 (MQ=255)
                                                                            aaTCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCgg                                                                                                             >  1:242131/1‑110 (MQ=255)
                                                                              tCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCGGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCg                                                                    <  4:634734/149‑1 (MQ=255)
                                                                              tCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCg                                                                    <  3:464646/149‑1 (MQ=255)
                                                                                cTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCggg                                                                  <  2:603661/149‑1 (MQ=255)
                                                                                 tgtgGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGc                                                                 >  1:771569/1‑149 (MQ=255)
                                                                                 tGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGct                                                                                                      <  1:834876/112‑1 (MQ=255)
                                                                                 tGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGct                                                                                                      >  2:834876/1‑112 (MQ=255)
                                                                                 tGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCt                                                                                         <  1:388459/125‑1 (MQ=255)
                                                                                 tGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCt                                                                                         >  2:388459/1‑125 (MQ=255)
                                                                                      gCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCAccc                                                            <  1:249272/149‑1 (MQ=255)
                                                                                        cTGATAAGACGCGTCAAGCGTCGCATCAGGCAGTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCag                                                          >  4:462866/1‑149 (MQ=255)
                                                                                            tAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgta                                                      >  4:364437/1‑149 (MQ=255)
                                                                                            tAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgta                                                      >  2:239405/1‑149 (MQ=255)
                                                                                             aaGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtag                                                     >  1:590781/1‑149 (MQ=255)
                                                                                              aGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagc                                                    <  2:136434/149‑1 (MQ=255)
                                                                                               gACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagcg                                                   >  2:207822/1‑149 (MQ=255)
                                                                                                aCGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCgg                                                                   >  4:662505/1‑132 (MQ=255)
                                                                                                aCGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCgg                                                                   <  3:662505/132‑1 (MQ=255)
                                                                                                aCGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagc                                                         >  2:35444/1‑142 (MQ=255)
                                                                                                aCGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagc                                                         <  1:35444/142‑1 (MQ=255)
                                                                                                 cgcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTg                                                 >  3:156885/1‑149 (MQ=255)
                                                                                                    gTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAcc                                              >  1:152312/1‑149 (MQ=255)
                                                                                                     tCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCt                                             >  2:291506/1‑149 (MQ=255)
                                                                                                        aGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTa                                          <  4:669065/149‑1 (MQ=255)
                                                                                                        aGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTa                                          <  1:231111/149‑1 (MQ=255)
                                                                                                        aGCGTCGCATCAGGCACTATGAACAAGATGCCGGACCATGCATCCGGCAATTCACTTAAAAACGCTGACCAGCTTCCAGCGCAACCGCTTCGTCACCAGCCTCAAACACATAAGTGGTGCTCGGGCCCCCCAGCGTAGCGTCACCTTTa                                          <  2:701414/149‑1 (MQ=255)
                                                                                                             cGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                                     <  2:185573/149‑1 (MQ=255)
                                                                                                             cGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                                     <  3:835764/149‑1 (MQ=255)
                                                                                                              gCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACt                                    >  4:713344/1‑149 (MQ=255)
                                                                                                                  cAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagc                                                    >  1:183300/1‑129 (MQ=255)
                                                                                                                  cAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagc                                                    <  2:183300/129‑1 (MQ=255)
                                                                                                                  cAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGa                                <  1:239405/149‑1 (MQ=255)
                                                                                                                     gCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGAccc                             >  1:494524/1‑149 (MQ=255)
                                                                                                                     gCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGATGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCATCGTCGCGTGACCTTTACTCACTGTGATcc                             >  2:330358/1‑149 (MQ=255)
                                                                                                                      cACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCa                            <  3:749638/149‑1 (MQ=255)
                                                                                                                      cACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCa                            <  3:491915/149‑1 (MQ=255)
                                                                                                                      cACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCa                            <  2:61604/149‑1 (MQ=255)
                                                                                                                        cTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGt                          <  3:429508/149‑1 (MQ=255)
                                                                                                                             aaCAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATc                              <  4:400639/140‑1 (MQ=255)
                                                                                                                             aaCAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATc                              >  3:400639/1‑140 (MQ=255)
                                                                                                                             aaCAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATc                              >  3:58020/1‑140 (MQ=255)
                                                                                                                                aaGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTcc                  >  2:636515/1‑149 (MQ=255)
                                                                                                                                    tGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTa              <  3:365963/149‑1 (MQ=255)
                                                                                                                                      ccGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGa            >  1:861267/1‑149 (MQ=255)
                                                                                                                                       cGGATCATGTATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGAc           <  3:176857/149‑1 (MQ=255)
                                                                                                                                        ggATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGAcc          >  1:143945/1‑149 (MQ=255)
                                                                                                                                              aTGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCaataa     >  3:641075/1‑149 (MQ=255)
                                                                                                                                              aTGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCaataa     >  1:587276/1‑149 (MQ=255)
                                                                                                                                              aTGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCaataa     <  1:201323/149‑1 (MQ=255)
                                                                                                                                              aTGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGACCCAGTTACCTTCCGGTAGACCaataa     <  2:494524/149‑1 (MQ=255)
                                                                                                                                                gCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATc   <  4:384261/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CTGTTTTGTTCGGTGGCTTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATCATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCG  >  NC_000913/4229294‑4229586

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