Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,229,438 Δ1 bp 5.5% 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. 5.5% 389.1 / 15.6 110intergenic (+76/+15)yjbB/pepEputative Na+/Pi‑cotransporter/peptidase E, alpha‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (55/49);  new base . (3/3);  total (58/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

Only 100 of 113 total aligned reads displayed.

TTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCG  >  NC_000913/4229311‑4229586
                                                                                                                                |                                                                                                                                                    
ttACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAAtt                                                                                                                                           >  1:1108631/1‑139 (MQ=255)
ttACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAAtt                                                                                                                                           <  2:1108631/139‑1 (MQ=255)
ttACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTaaaaa                                                                                                                                 >  4:190006/1‑149 (MQ=255)
  aCAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACg                                                                                                                               <  2:1412576/149‑1 (MQ=255)
   cAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCAc                                                                                                                                        >  2:1454471/1‑139 (MQ=255)
    aGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGt                                                                                                                             <  3:190006/149‑1 (MQ=255)
          cTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAg                                                                                                                       >  1:541436/1‑149 (MQ=255)
             gAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCtt                                                                                                                    >  4:878426/1‑149 (MQ=255)
                                      gggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGCatc                                                                                                                                            <  3:915209/100‑3 (MQ=25)
                                      gggacgAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCGGCatc                                                                                                                                            >  4:915209/1‑98 (MQ=25)
                                           gAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGATCAAGATGCCGGAT‑CATGCATCCGGCAATTCACt                                                                                                                                       <  3:1245616/100‑1 (MQ=255)
                                           gAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACt                                                                                                                                       >  4:1245616/1‑100 (MQ=255)
                                                ttAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCAc                                                                                                                                        <  4:714846/94‑1 (MQ=255)
                                                   atatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTc                                                                                                                   >  4:291409/1‑112 (MQ=255)
                                                   atatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTc                                                                                                                   <  3:291409/112‑1 (MQ=255)
                                                    tatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATGAg                                                                             >  1:732821/1‑149 (MQ=255)
                                                    tatTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAg                                                                             >  1:1439999/1‑149 (MQ=255)
                                                        gtgAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggt                                                                         <  2:732821/149‑1 (MQ=255)
                                                        gtgAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggt                                                                         <  2:1439999/149‑1 (MQ=255)
                                                        gtgAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggt                                                                         <  1:1102830/149‑1 (MQ=255)
                                                         tgAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggtg                                                                        >  2:77183/1‑149 (MQ=255)
                                                            aTCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTc                                                                     >  2:800178/1‑149 (MQ=255)
                                                                tGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                                     >  2:382893/1‑129 (MQ=255)
                                                                tGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                                     <  1:382893/129‑1 (MQ=255)
                                                                  tAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACgg                                                                                                                              >  1:824668/1‑86 (MQ=255)
                                                                    ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                                     <  3:459137/125‑1 (MQ=255)
                                                                    ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                                     >  4:459137/1‑125 (MQ=255)
                                                                    ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAAcaca                                                                                 <  3:1145975/129‑1 (MQ=255)
                                                                    ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCAcc                                                             >  3:386122/1‑149 (MQ=255)
                                                                      ccTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACg                                                                                                                               >  2:981193/1‑81 (MQ=255)
                                                                      ccTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACg                                                                                                                               <  1:981193/81‑1 (MQ=255)
                                                                      ccTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCa                                                           >  4:254119/1‑149 (MQ=255)
                                                                       cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                                     >  4:609348/1‑122 (MQ=255)
                                                                       cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                                     >  2:1131782/1‑122 (MQ=255)
                                                                       cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                                     <  1:1131782/122‑1 (MQ=255)
                                                                       cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                                     <  3:609348/122‑1 (MQ=255)
                                                                       cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggt                                                                         <  1:669408/134‑1 (MQ=255)
                                                                       cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtggt                                                                         >  2:669408/1‑134 (MQ=255)
                                                                       cTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCag                                                          >  2:289475/1‑149 (MQ=255)
                                                                        tGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttc                                                                                                    >  1:463968/1‑106 (MQ=255)
                                                                        tGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttc                                                                                                    <  2:463968/106‑1 (MQ=255)
                                                                         gATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcg                                                        >  4:1194746/1‑149 (MQ=255)
                                                                           tAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgta                                                      <  1:77183/149‑1 (MQ=255)
                                                                            aaGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGctt                                                                                                     <  3:316671/101‑1 (MQ=255)
                                                                            aaGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGctt                                                                                                     >  4:316671/1‑101 (MQ=255)
                                                                            aaGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtag                                                     <  4:386122/149‑1 (MQ=255)
                                                                              gACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagcg                                                   >  3:1345889/1‑149 (MQ=255)
                                                                               aCGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                                     >  3:1049191/1‑114 (MQ=255)
                                                                               aCGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                                     <  4:1049191/114‑1 (MQ=255)
                                                                               aCGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtg                                                                           <  4:1321338/124‑1 (MQ=255)
                                                                               aCGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CAAGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAgtg                                                                           >  3:1321338/1‑124 (MQ=255)
                                                                                  cgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAc                                               <  1:289475/149‑1 (MQ=255)
                                                                                  cgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAc                                               <  3:254119/149‑1 (MQ=255)
                                                                                       aGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTa                                          <  4:1345889/149‑1 (MQ=255)
                                                                                        gCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTAc                                         >  4:745074/1‑149 (MQ=255)
                                                                                          gTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTc                                       >  4:301450/1‑149 (MQ=255)
                                                                                           tCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCa                                      <  3:878426/149‑1 (MQ=255)
                                                                                            cGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                                     >  4:302076/1‑149 (MQ=255)
                                                                                              cATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACtg                                   <  1:800178/149‑1 (MQ=255)
                                                                                               aTCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACtgt                                  >  3:1399664/1‑149 (MQ=255)
                                                                                                  aGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGAt                               >  3:1243108/1‑149 (MQ=255)
                                                                                                   ggCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAc                                               <  2:747551/132‑1 (MQ=255)
                                                                                                   ggCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAc                                               >  1:747551/1‑132 (MQ=255)
                                                                                                    gCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGa                                <  1:611851/146‑1 (MQ=255)
                                                                                                    gCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGa                                >  2:611851/1‑146 (MQ=255)
                                                                                                     cACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCa                            <  3:487231/149‑1 (MQ=255)
                                                                                                      aCTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                                     <  1:1180971/91‑1 (MQ=255)
                                                                                                        tATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGtt                         <  4:872550/149‑1 (MQ=255)
                                                                                                        tATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGtt                         <  3:1194746/149‑1 (MQ=255)
                                                                                                        tATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGtt                         >  3:872550/1‑149 (MQ=255)
                                                                                                         aTGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGt                          <  2:1001903/147‑1 (MQ=255)
                                                                                                         aTGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGt                          >  1:1001903/1‑147 (MQ=255)
                                                                                                            aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAc                                               >  2:1143550/1‑123 (MQ=255)
                                                                                                            aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGAc                                               <  1:1143550/123‑1 (MQ=255)
                                                                                                            aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                                     >  4:622161/1‑133 (MQ=255)
                                                                                                            aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                                     <  1:1088448/133‑1 (MQ=255)
                                                                                                            aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                                     >  2:1088448/1‑133 (MQ=255)
                                                                                                            aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                                     <  3:622161/133‑1 (MQ=255)
                                                                                                            aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCCTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagc                                                         >  4:1047183/1‑113 (MQ=255)
                                                                                                            aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCCTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagc                                                         <  3:1047183/113‑1 (MQ=255)
                                                                                                             aCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCtt                    >  4:1320296/1‑149 (MQ=255)
                                                                                                              cAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                                     <  1:422202/83‑1 (MQ=255)
                                                                                                              cAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTaaa                                                                                     >  2:422202/1‑83 (MQ=255)
                                                                                                              cAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTc                   <  4:1243108/149‑1 (MQ=255)
                                                                                                               aaGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTcc                  >  1:151725/1‑149 (MQ=255)
                                                                                                                 attGCCGGATGCATG‑ATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCagcgtagcgt                                                  <  3:952547/113‑1 (MQ=255)
                                                                                                                  ttGCCGGATGCATG‑ATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTAcc                      <  2:951228/141‑1 (MQ=255)
                                                                                                                  aTGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGt               <  3:1466812/149‑1 (MQ=255)
                                                                                                                  aTGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGt               >  2:1204540/1‑149 (MQ=255)
                                                                                                                    gCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                                     <  2:661130/125‑1 (MQ=255)
                                                                                                                    gCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAg             >  1:861453/1‑149 (MQ=255)
                                                                                                                       ggAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                                     <  1:555687/122‑1 (MQ=255)
                                                                                                                       ggAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCAc                                     >  2:555687/1‑122 (MQ=255)
                                                                                                                             aTGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCaataa     >  1:881684/1‑149 (MQ=255)
                                                                                                                              tGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCaataat    >  4:1192312/1‑149 (MQ=255)
                                                                                                                               gCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATc   <  2:881684/149‑1 (MQ=255)
                                                                                                                               gCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATc   <  3:745074/149‑1 (MQ=255)
                                                                                                                                cATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTa                                                                                       >  1:992722/1‑64 (MQ=255)
                                                                                                                                cATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTa                                                                                       <  2:992722/64‑1 (MQ=255)
                                                                                                                                cATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCg  >  4:1400419/1‑149 (MQ=255)
                                                                                                                                |                                                                                                                                                    
TTACAGTGTGCTGGAACAGCCGGATGAAGATGAAGGACGGGACGAGTATTAATATTGTGAATCCTGTAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCG  >  NC_000913/4229311‑4229586

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