Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,325,821 G→T 5.7% intergenic (‑79/+579) yjdN ← / ← yjdM metalloprotein superfamily protein/zinc‑ribbon family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,325,8210GT 5.7% 253.1 / 3.3 87intergenic (‑79/+579)yjdN/yjdMmetalloprotein superfamily protein/zinc‑ribbon family protein
Reads supporting (aligned to +/- strand):  ref base G (51/31);  new base T (2/3);  total (53/34)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

TTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAAC  >  NC_000913/4325680‑4325969
                                                                                                                                             |                                                                                                                                                    
ttGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAGAGAGAGGTAGGGACTTAACGGCATGATGGTGACCACAGTTAAGAGAAGCCAGTTAAGTGGAGTCCGGAATTATTACAGATGAGAGGCGGCCTGATTAGAGGCGCGAGGGt                                                                                                                                               >  2:356456/1‑149 (MQ=255)
ttGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGt                                                                                                                                               >  2:229807/1‑149 (MQ=255)
ttGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGt                                                                                                                                               >  2:584376/1‑149 (MQ=255)
   ataataGGCAATCGCGTCGGAACAGTTCCCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGTCCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGACGGCCTGATAAGACGCTCGAGCGTCGc                                                                                                                                            >  1:345507/1‑149 (MQ=255)
    taataGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCTCGCGCGTctta                                                                                                                                           <  1:770553/149‑4 (MQ=255)
      ataGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCTCGCGCGTcttatc                                                                                                                                         <  1:615601/149‑6 (MQ=255)
          gCAAGCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGCCTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCTCGCGCGTCTTATCAGGc                                                                                                                                     <  2:147569/149‑1 (MQ=255)
                gcgTCGGAACAGTTACCGGCAAAAGAGGGGTAGGGACTGAACGGCGTGCTGGGGACCTCAGTTAAGAGAAGCCGGTTTAGGGGAGTTCGCAATTTATACAGATGAGAGTCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCgg                                                                                                                               <  2:302208/149‑1 (MQ=38)
                 cgTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGGCTGATAAGACGCGCGAGGGTCGCATCAGGCAGTGGGc                                                                                                                              >  2:148940/1‑149 (MQ=255)
                    cGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGGCAGGCGGCCTGATAAGACGCGGGAGCGTCGCATCAGCCAGTACGCAcc                                                                                                                           >  2:307006/1‑148 (MQ=255)
                      gAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGt                                                                                                                         >  2:250893/1‑149 (MQ=255)
                          aGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGcc                                                                                                                     <  1:626034/149‑1 (MQ=255)
                           gTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGCTAAGAGAAGCCAGTTAAGTGTAGTGCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCg                                                                                                                    <  2:639336/149‑1 (MQ=255)
                                  gCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGGCCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATgcgg                                                                                                             >  2:725143/1‑149 (MQ=255)
                                   cgaaaGAGGGGTAGGGACTGAACGGCATGATGGGGCCCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTAGCCCGGTGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATgcggc                                                                                                            <  2:435312/147‑1 (MQ=37)
                                      aaGAGAGGTAGGGACTTAACGGCATGATGGTGCCCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTa                                                                                                         <  1:210012/149‑1 (MQ=255)
                                      aaGAGAGGTAGGGACTTAACGGCATGATGGTGACATCAGTTAAGAGAAGCACGCTACGCGCAGCTTGACCTTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCAATGTTGGCGGGTGCGGCGTa                                                                                                         >  2:89051/1‑149 (MQ=25)
                                        gagagGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTaaa                                                                                                       <  2:464026/149‑1 (MQ=255)
                                         agagGTAGGGACTTAACCGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGGTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAAc                                                                                                      <  1:47391/149‑1 (MQ=255)
                                           agGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGc                                                                                                    <  1:639529/149‑1 (MQ=255)
                                               aGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGGAAACGCCTTa                                                                                                >  2:411667/1‑149 (MQ=255)
                                                  gACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                             >  1:700371/1‑149 (MQ=255)
                                                  ctcTTCACGGCATTACGGTGACCTCAGTTAAGGCAAGCCAGTTAATTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                             <  2:662086/147‑1 (MQ=25)
                                                   aCTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                            >  1:702239/1‑149 (MQ=255)
                                                   aCTGAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                            <  2:147815/149‑1 (MQ=255)
                                                          ggCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACg                                                                                     >  1:610252/1‑149 (MQ=255)
                                                           gCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGGGGAGTAAACGCCTTATCGGGCCTACgg                                                                                    >  2:92403/1‑149 (MQ=255)
                                                             atgatgGCGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATCCTGATGAGAGGCGGCCTGTTAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGGCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAg                                                                                  <  2:7024/149‑1 (MQ=25)
                                                              tgatgGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGgggg                                                                                 >  2:330195/1‑149 (MQ=255)
                                                              tgatgGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAgg                                                                                 >  1:177632/1‑149 (MQ=255)
                                                                atgGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCGCTGTTGCCGGATGCGGCGTCAACGCCTTATCCGGCCAGTGGAgggc                                                                               >  2:341614/1‑148 (MQ=25)
                                                                atgGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGt                                                                               <  2:618778/149‑1 (MQ=255)
                                                                   gTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTGCGGAGGGTGCg                                                                            >  1:578043/1‑149 (MQ=255)
                                                                   gTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                            >  1:578103/1‑149 (MQ=255)
                                                                   gTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTGGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCg                                                                            >  1:104945/1‑149 (MQ=255)
                                                                    tGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATGTATACAGTTCAGGTGCGGCTTGATAATACTCGCGAGCGTCTCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCgg                                                                           <  2:44641/149‑1 (MQ=21)
                                                                        cTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGGAt                                                                       >  1:6902/1‑149 (MQ=255)
                                                                        cTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAt                                                                       >  1:168705/1‑149 (MQ=255)
                                                                        cTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAt                                                                       >  2:258443/1‑149 (MQ=255)
                                                                        cTCAGTTAAGAGAAGCCAGTTAAGTGTAGGTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAt                                                                       >  2:460518/1‑149 (MQ=255)
                                                                          cAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAttt                                                                     <  1:507643/149‑1 (MQ=255)
                                                                           aGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTg                                                                    >  1:260033/1‑149 (MQ=255)
                                                                           aGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGCGCGGGaagatg                                                                    >  2:470001/1‑145 (MQ=255)
                                                                           aGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAAGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGATAGTCAACACAGTTGACGGATGCGGCGTCAACGCCTTATCCGGCCTACGGAGGGTGACGGaaatgg                                                                    >  2:324035/1‑145 (MQ=21)
                                                                           aGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTGATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGCAGTTg                                                                    >  2:469998/1‑149 (MQ=255)
                                                                           aGTTAAGAGAAGCCAGTTAAGTGTAGTTCAAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTg                                                                    <  2:155357/149‑1 (MQ=255)
                                                                              tAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAg                                                                 <  1:309692/149‑1 (MQ=255)
                                                                                   gaAGCCAGTTTAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTGGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTg                                                            <  2:674795/149‑1 (MQ=255)
                                                                                   gaAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATTAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGACGTCGGCACTGTTGACGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTg                                                            <  2:353086/149‑1 (MQ=255)
                                                                                    aaGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGa                                                           >  2:436826/1‑149 (MQ=255)
                                                                                    aaGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGa                                                           <  2:76358/149‑1 (MQ=255)
                                                                                     aGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGAt                                                          <  2:581055/149‑1 (MQ=255)
                                                                                              aGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACgggc                                                 >  1:32038/1‑149 (MQ=255)
                                                                                              aGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATGTGTAGGCCTGATAAGACgcgc                                                 >  1:42368/1‑149 (MQ=255)
                                                                                              aGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGAGCGGGAATGTTTAGGCCGGATAAGACgcgc                                                 >  1:357931/1‑149 (MQ=255)
                                                                                              aGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCGACGGAGGGAGCGGGAATTTGTAGGCCTGATAAGACgcgc                                                 >  1:644367/1‑149 (MQ=255)
                                                                                              aGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGACCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAATTTGTAGGCCTGATAAGACgcgc                                                 >  1:737269/1‑149 (MQ=255)
                                                                                                tgtAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCaa                                               <  1:201953/149‑1 (MQ=255)
                                                                                                  tAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCTAGCGTCGAAAAAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGc                                             >  1:696245/1‑149 (MQ=255)
                                                                                                  tAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGc                                             >  1:752159/1‑149 (MQ=255)
                                                                                                  tAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCACATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGAAAAGACGCGCAAGc                                             >  2:588024/1‑149 (MQ=255)
                                                                                                     ttCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTc                                          >  2:439826/1‑149 (MQ=255)
                                                                                                     ttCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTc                                          >  2:206420/1‑149 (MQ=255)
                                                                                                     ttCGAAATTTACACAGACGAGCGGCGGCCGGATAAGACGCGCGCGCGTCGCATAAGGCAGCCGGCACTGTGGACTGATACGGCGCCAACGCCTTAGCCGCCCGAAGGAGGTAGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGctgc                                          >  2:227480/1‑146 (MQ=9)
                                                                                                      tCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCg                                         <  2:337895/149‑1 (MQ=255)
                                                                                                       cGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGc                                        >  2:34403/1‑149 (MQ=255)
                                                                                                       cGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGc                                        >  1:75756/1‑149 (MQ=255)
                                                                                                         aaaTTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCGAGAAAAGACGCGACAGCGGCGCAt                                      >  2:197313/1‑149 (MQ=17)
                                                                                                            tttATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                   <  2:747734/149‑1 (MQ=255)
                                                                                                            tttATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                   <  2:747703/149‑1 (MQ=255)
                                                                                                            cttatcaaGATGCGAGGCCGCCTGATAACACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                   <  2:153551/142‑1 (MQ=14)
                                                                                                             ttATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAgg                                  >  2:427784/1‑149 (MQ=255)
                                                                                                               ataCAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCa                                <  2:610252/149‑1 (MQ=255)
                                                                                                                  cAGATGAGAGGCGGCCTGATAAGACGCGCGAGCATCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCATGCGTCGCATCAGGCAGTc                             >  1:500082/1‑149 (MQ=255)
                                                                                                                   aGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCg                            >  1:65316/1‑149 (MQ=255)
                                                                                                                    gATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCgg                           >  1:387201/1‑149 (MQ=255)
                                                                                                                         gagGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCa                      >  1:700982/1‑149 (MQ=40)
                                                                                                                              tgcCTGATGCGACGCTCGCGCGTCTTATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGcc                 >  2:583357/2‑149 (MQ=11)
                                                                                                                                cttGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGAGAATTTGTACGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCgg               <  2:477845/147‑1 (MQ=255)
                                                                                                                                ccTGATAAGACGCGCCAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCgg               <  2:545782/149‑1 (MQ=21)
                                                                                                                                 cTGATAAGACGCGCGAGCGTCGCAGCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGAGGTACATTGGAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGa              <  1:307006/149‑1 (MQ=11)
                                                                                                                                   gATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAg                                                       >  2:353223/1‑106 (MQ=38)
                                                                                                                                   gATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATg            >  1:558458/1‑149 (MQ=38)
                                                                                                                                     tAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGcc                                                                                                   >  2:621410/1‑60 (MQ=255)
                                                                                                                                     gaaGACGCGCGAGCGTCGCACCAGGCAGCCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATgcg          <  1:98125/148‑1 (MQ=255)
                                                                                                                                      aaGACGCGCGATCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATgcgg         >  2:593624/1‑149 (MQ=38)
                                                                                                                                       aGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTCGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATgcggc        <  1:513569/149‑1 (MQ=38)
                                                                                                                                           gcgcgAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcg                                                    >  2:293385/1‑101 (MQ=38)
                                                                                                                                           gcgcgAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcg                                                    >  2:51365/1‑101 (MQ=38)
                                                                                                                                             gcgAGCGTCGCATGAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAAc  <  1:222159/149‑1 (MQ=38)
                                                                                                                                             |                                                                                                                                                    
TTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAAC  >  NC_000913/4325680‑4325969

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