Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,315,964 T→C 6.3% intergenic (‑406/+252) 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
*CP0092734,315,9640TC 6.3% 429.2 / 12.7 155intergenic (‑406/+252)yjdN/yjdMmetalloprotein superfamily protein/zinc‑ribbon family protein
Reads supporting (aligned to +/- strand):  ref base T (88/52);  new base C (5/4);  total (98/57)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

TGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGG  >  CP009273/4315818‑4316113
                                                                                                                                                  |                                                                                                                                                     
         gATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAAATTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGATGGCACCGTTGCCGGATgc                                                                                                                                           >  1:2561413/1‑150 (MQ=12)
          tgaaGACGCGCTAGCGTAGCAACAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTTCGGATTGTGCGGGATTGTGTCGGCCTGATAAGACGCGNTAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATgcg                                                                                                                                          <  1:4789424/148‑1 (MQ=14)
          aTAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATgcg                                                                                                                                          <  2:1277645/150‑1 (MQ=255)
           tAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATgcgg                                                                                                                                         <  1:1074096/150‑1 (MQ=255)
               aCGCGCAAGCGTCGCATCAGGCAGTCGGCACCGGTGCCGGATGCGGCGTAAACGCCTTATCCTGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCCAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTa                                                                                                                                     <  2:704525/150‑1 (MQ=255)
                   gcACGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCGTATCCGGCCTTAGGATGGTGAGGGAATGTGTTGGCGTGATCAGACGCGCAATAGTCGCATCAGGCAGTTGGCTCTGTTGGCGGATGCGGCGGTAACg                                                                                                                                 >  1:2678665/1‑150 (MQ=14)
                     aaGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGcc                                                                                                                               >  2:4753378/1‑150 (MQ=39)
                     aaGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGcc                                                                                                                               >  2:2704908/1‑150 (MQ=39)
                     aaGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACCGTTGCCGGATGCGGCGTAAACGcc                                                                                                                               >  2:4753484/1‑150 (MQ=21)
                      aGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCt                                                                                                                              <  1:4058547/150‑1 (MQ=39)
                       gCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCtt                                                                                                                             <  1:2700206/150‑1 (MQ=39)
                       gCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCtt                                                                                                                             <  1:2700192/150‑1 (MQ=39)
                           cGCATCAGGCAGTCGGCACCTTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                         <  1:1059875/150‑1 (MQ=37)
                            gCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                        >  1:4223836/1‑150 (MQ=21)
                            gCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGAGCGGGAATTTGTAGGCCTGATAAGACGCGCACGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                        <  1:2163800/150‑1 (MQ=255)
                            gCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACGGATGGGGCGGGACTTTGTAGGCCTGATAGGACGAGCAAGCGTCGCATCAGGCAGTTGGCACGGTTGCCGGATGCGGCGTGAGCGCCTGGTCCg                                                                                                                        >  2:348666/1‑150 (MQ=11)
                                 aGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTCTGCGGGCTTCGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                   <  2:4435801/150‑1 (MQ=11)
                                     aGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCATGGGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGa                                                                                                               <  2:3355287/150‑1 (MQ=255)
                                     aGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTCTCCGGCCTACGGa                                                                                                               >  2:1257670/1‑150 (MQ=255)
                                      gTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAg                                                                                                              <  1:3640904/150‑1 (MQ=37)
                                       tCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAgg                                                                                                             <  1:719188/150‑1 (MQ=40)
                                        cGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGGTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAggg                                                                                                            >  2:4546209/1‑150 (MQ=255)
                                            aCCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                                                        >  2:1398163/1‑150 (MQ=39)
                                            aCCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACAGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                                                        >  1:1371226/1‑150 (MQ=37)
                                             ctgTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCCCTGTTGCCGGATgcggc                                                                                                                                        >  2:7570/3‑117 (MQ=255)
                                             ctgTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATgcggc                                                                                                                                        <  1:7600/115‑1 (MQ=37)
                                             ctgTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCgg                                                                                                       <  2:3422892/148‑1 (MQ=25)
                                              cGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCggg                                                                                                      >  1:1771737/1‑150 (MQ=17)
                                               gTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGa                                                                                                     >  1:2957462/1‑150 (MQ=37)
                                                ttGCCGGATGCGGCGTATACGCCTTCTCCCGCCATCGGAGGTTGCTGGATTTTGTAGGCCTGATAAGACGCGCTAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGaa                                                                                                    <  2:2157935/150‑1 (MQ=11)
                                                ttGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGaa                                                                                                    >  1:1652242/1‑150 (MQ=37)
                                                 tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAt                                                                                                   <  2:2569551/150‑1 (MQ=37)
                                                  gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAtt                                                                                                  >  2:1409247/1‑150 (MQ=37)
                                                        tGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAgg                                                                                            <  1:3360461/150‑1 (MQ=37)
                                                         gcggcgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGc                                                                                           <  1:1589629/150‑1 (MQ=17)
                                                          tggcgTAAACGCCTTATCCGGCCTACGGATGGTGCGGGACTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGcc                                                                                          >  1:4369034/2‑150 (MQ=37)
                                                              gTAAACGCCTTATCCGGCCTACGGATGGTGCGGGACTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGATGGCACCGTTGCCGGATGCGGCGGCAACGCCTTATCCGGCCTACGGAGGGTGAGCGAAGTGGTAGGCCAGAt                                                                                      >  2:3937183/1‑150 (MQ=11)
                                                              gTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGAt                                                                                      >  2:3937375/1‑150 (MQ=37)
                                                              gTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGAt                                                                                      >  2:1911325/1‑150 (MQ=37)
                                                                 aaCGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAg                                                                                   <  1:2602015/150‑1 (MQ=37)
                                                                       ttATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACgcgc                                                                             >  1:2861823/1‑150 (MQ=37)
                                                                       ttATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACgcgc                                                                             >  2:4102785/1‑150 (MQ=37)
                                                                       ttATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACgcgc                                                                             >  2:402638/1‑150 (MQ=37)
                                                                        tATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCa                                                                            >  1:2822444/1‑150 (MQ=37)
                                                                         aTCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCaa                                                                           >  2:4464767/1‑150 (MQ=37)
                                                                         aTCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCaa                                                                           >  2:4653520/1‑150 (MQ=37)
                                                                          tCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAg                                                                          <  2:588726/150‑1 (MQ=17)
                                                                             ggCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGt                                                                       >  1:5085585/1‑150 (MQ=37)
                                                                                cTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGATGGCACAGTTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACCGAGGGTGGGGGGAGAGGGGGGGCTGATAAGACGCGCAAGCGTCgg                                                                    >  2:1919124/1‑149 (MQ=11)
                                                                                 tACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCa                                                                   >  2:4051981/1‑150 (MQ=37)
                                                                                   cGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGACTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                                                                 >  2:5093565/1‑150 (MQ=37)
                                                                                    ggATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa                                                                >  1:2645141/1‑150 (MQ=37)
                                                                                    ggATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa                                                                >  1:3375209/1‑150 (MQ=37)
                                                                                     gATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                                               >  2:3863305/1‑150 (MQ=37)
                                                                                      aTGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAgg                                                              >  1:1780455/1‑150 (MQ=35)
                                                                                          tGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGt                                                          >  1:4589274/1‑150 (MQ=35)
                                                                                            cGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCg                                                        >  2:2941500/1‑150 (MQ=37)
                                                                                               gAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGCTGGCa                                                     >  1:568951/1‑150 (MQ=17)
                                                                                                        ggCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCg                                            >  2:2612650/1‑150 (MQ=37)
                                                                                                          ccTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCtt                                                                                                                             >  2:2112908/1‑67 (MQ=40)
                                                                                                          ccTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCtt                                                                                                                             <  1:2112908/67‑1 (MQ=40)
                                                                                                           cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCACGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGAt                                         >  2:4114628/1‑150 (MQ=37)
                                                                                                           cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGAt                                         >  2:4114569/1‑150 (MQ=37)
                                                                                                           cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGAt                                         >  2:115129/1‑150 (MQ=37)
                                                                                                           cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGAt                                         >  2:3663550/1‑150 (MQ=37)
                                                                                                           cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGAt                                         <  1:3051725/150‑1 (MQ=37)
                                                                                                           cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATGTGTAGGACTGATAAGACGCGCAAGAGTCGCATCAGGCAGTCGGCACCGTTGCCGGag                                         >  1:1325107/1‑149 (MQ=255)
                                                                                                           cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAAAGCCTTATACGGCCTACGGAGGGTGCGGGAAAATGTATGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGGCggtt                                         >  1:1492178/1‑148 (MQ=11)
                                                                                                             gATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgc                                       >  1:2789314/1‑150 (MQ=37)
                                                                                                             gATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGAGgc                                       >  1:4849239/1‑150 (MQ=255)
                                                                                                             gATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgc                                       <  2:340913/150‑1 (MQ=17)
                                                                                                             gATAAGACGAGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgc                                       >  1:1423177/1‑150 (MQ=37)
                                                                                                                aaGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgcggc                                    <  1:3726973/150‑1 (MQ=37)
                                                                                                                aaGACGCGCAAGCGTCGCATCAGGCAGTTCACACTGTTGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTGGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGACGGCACCGTTGCCGGAttctgc                                    >  1:2961244/1‑145 (MQ=11)
                                                                                                                      cgcAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCACGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAAc                              >  2:5145068/1‑150 (MQ=37)
                                                                                                                      cgcAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAAc                              <  1:2373553/150‑1 (MQ=37)
                                                                                                                        cAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGAt                                         <  1:2699311/137‑1 (MQ=37)
                                                                                                                         aaGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGcc                           >  2:1055882/1‑150 (MQ=35)
                                                                                                                         aaGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGGAAAcccc                           >  1:2350619/1‑150 (MQ=255)
                                                                                                                          aGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCt                          >  2:4294959/1‑150 (MQ=34)
                                                                                                                             gTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTAt                       >  2:4973938/1‑150 (MQ=34)
                                                                                                                                 cATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                   <  1:3064285/150‑1 (MQ=34)
                                                                                                                                    cAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTCAACGCCTTTTCCGGCCt                >  1:4862944/1‑150 (MQ=21)
                                                                                                                                     tggCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa               >  2:332132/2‑150 (MQ=21)
                                                                                                                                        cAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCTGATGCGGCGTAAACGCCTTATCCGGCCTACgg            <  2:1777296/150‑1 (MQ=21)
                                                                                                                                         aGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGACTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTCCGGa           >  2:3331155/1‑150 (MQ=21)
                                                                                                                                          gTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAg          >  2:4415863/1‑150 (MQ=34)
                                                                                                                                           ttGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCCCCGTTGCCGGATGCGGCGGACCCGCCTTCTCCGGCCTACGGAgg         >  2:4087088/1‑150 (MQ=11)
                                                                                                                                            tGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAggg        <  1:4464767/150‑1 (MQ=32)
                                                                                                                                            tGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAggg        <  1:4653520/150‑1 (MQ=32)
                                                                                                                                            tGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTAGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAggg        >  2:4560199/1‑150 (MQ=17)
                                                                                                                                            tGGCACAGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAggg        <  1:3450158/150‑1 (MQ=17)
                                                                                                                                             ggCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGt       <  1:631889/150‑1 (MQ=17)
                                                                                                                                              gCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTg      <  1:763559/150‑1 (MQ=17)
                                                                                                                                              gCACTGTTGACGGATGCGGCGTCAAAGCCTTATCCGGCCTTCGGAGTGTTCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTg      <  2:675788/150‑1 (MQ=11)
                                                                                                                                                aCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCt                          >  1:394628/1‑128 (MQ=32)
                                                                                                                                                aCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                   <  1:3065218/135‑1 (MQ=32)
                                                                                                                                                 cTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGACTTGGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGGCAACGCCTTATCCGGCCTACGGAGGGTGCgg   >  1:1563251/1‑150 (MQ=255)
                                                                                                                                                 cTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCgg   >  2:4716045/1‑150 (MQ=32)
                                                                                                                                                 cTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCgg   >  1:3667562/1‑150 (MQ=32)
                                                                                                                                                  |                                                                                                                                                     
TGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGG  >  CP009273/4315818‑4316113

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: