Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,750,073 G→T 7.3% intergenic (+13/+42) yfjD → / ← grpE UPF0053 family inner membrane protein/heat shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,750,0730GT 7.3% 278.7 / 10.8 82intergenic (+13/+42)yfjD/grpEUPF0053 family inner membrane protein/heat shock protein
Reads supporting (aligned to +/- strand):  ref base G (44/32);  new base T (3/3);  total (47/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

GAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTC  >  NC_000913/2749926‑2750219
                                                                                                                                                   |                                                                                                                                                  
gAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCgg                                                                                                                                                   >  1:803221/1‑149 (MQ=255)
gAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCgg                                                                                                                                                   >  1:364214/1‑149 (MQ=255)
      aTCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTACGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCg                                                                                                                                             >  2:499849/1‑149 (MQ=255)
      aTCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCg                                                                                                                                             >  1:865080/1‑149 (MQ=255)
         ccTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCAt                                                                                                                                          >  4:1118288/1‑149 (MQ=255)
         ccTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCAt                                                                                                                                          >  4:50046/1‑149 (MQ=255)
              cGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAgg                                                                                                                                     <  2:364214/149‑1 (MQ=255)
               gCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCTGCATGCGAATGCCggg                                                                                                                                    <  2:582894/149‑1 (MQ=255)
               gCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAggg                                                                                                                                    <  2:865080/149‑1 (MQ=255)
                 aGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCTGCATGCGAATGCCGGGcc                                                                                                                                  >  4:257184/1‑149 (MQ=255)
                 aGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCTGCATGCGAATGCCGGGcc                                                                                                                                  <  3:257184/149‑1 (MQ=255)
                  ggCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTc                                                                                                                                              <  3:116403/136‑1 (MQ=255)
                  ggCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTc                                                                                                                                              >  4:116403/1‑136 (MQ=255)
                            tgcgtATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAAttatta                                                                                                                       <  3:50046/149‑1 (MQ=255)
                                tATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGt                                                                                                                                <  2:107082/136‑1 (MQ=255)
                                tATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGt                                                                                                                                >  1:107082/1‑136 (MQ=255)
                                tATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGaa                                                                                                                   >  4:256297/1‑149 (MQ=255)
                                tATTGGCGAGTACGATATCGATATTCTCGACGTACACGACAATATGATTAAGCAGGTAAACGCTTTTGCTGCGAAACCGCTGCGCGAGAGTGTGGCGGATTAACGAAAACGGCCCGGCCTTCGCATGCAGGGCCGCGTATTATTACGaa                                                                                                                   >  4:474497/1‑149 (MQ=255)
                                  ttGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAg                                                                                                                 >  2:602614/1‑149 (MQ=255)
                                            cgatatcgatatTCTCGACGTACAGGACAATATGATTAAGCAGTTAAAAGTTTTTCCTGTGAAACCACTGCGCAAGAGTGTGACTGAGTAACGAAAACGGCCCGGCATTCGCAGGCAGGGCCGTGAATTATTACGAAAGCAGAAATTaa                                                                                                       <  2:901529/149‑1 (MQ=255)
                                            cgatatcgatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTaa                                                                                                       <  1:1096754/149‑1 (MQ=255)
                                             gatatcgatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGaaa                                                                                                                  >  4:188974/1‑137 (MQ=255)
                                             gatatcgatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGaaa                                                                                                                  <  3:188974/137‑1 (MQ=255)
                                                atcgatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCtt                                                                                                   >  4:604998/1‑149 (MQ=255)
                                                 tcgatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCttt                                                                                                  <  1:499849/149‑1 (MQ=255)
                                                   gatatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTg                                                                                                >  1:431710/1‑149 (MQ=255)
                                                    atatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCtt                                                                                                   >  4:346739/1‑145 (MQ=255)
                                                    atatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCtt                                                                                                   <  3:346739/145‑1 (MQ=255)
                                                     tatTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCt                                                                                              >  4:631255/1‑149 (MQ=255)
                                                      atTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCtt                                                                                             >  1:327571/1‑149 (MQ=255)
                                                        tctcGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTc                                                                                           <  2:327571/149‑1 (MQ=255)
                                                                  cAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTaa                                                                                 >  4:629315/1‑149 (MQ=255)
                                                                   aGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAAc                                                                                >  3:817406/1‑149 (MQ=255)
                                                                     gACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCa                                                                              >  2:148393/1‑149 (MQ=255)
                                                                     gACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCa                                                                              >  4:344951/1‑149 (MQ=255)
                                                                     gACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCa                                                                              >  1:634723/1‑149 (MQ=255)
                                                                      aCAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCAt                                                                             <  3:631255/149‑1 (MQ=255)
                                                                             gATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCAt                                                                                                                                          >  2:688190/1‑81 (MQ=255)
                                                                             gATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCAt                                                                                                                                          <  1:688190/81‑1 (MQ=255)
                                                                              gtttcacAGGAAAAACTTTTACCTGCGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCAc                                                                     >  3:417764/7‑149 (MQ=255)
                                                                                  aGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAg                                                                                                      <  3:1046683/112‑1 (MQ=255)
                                                                                  aGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAg                                                                                                      >  4:1046683/1‑112 (MQ=255)
                                                                                  aGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAAt                                                                 <  3:256297/149‑1 (MQ=255)
                                                                                   gCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATc                                                                >  1:574426/1‑149 (MQ=255)
                                                                                    cAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCg                                                               >  3:980059/1‑149 (MQ=255)
                                                                                     aGGTAAAAGTTTCTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGGGAATTATTACGAAAGCAGAAATTAAGATTTTGCTTTCGCTACAGTAACTATCGCCGCACGAATCGt                                                              >  4:968609/1‑149 (MQ=255)
                                                                                      ggTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTa                                                             >  3:1115301/1‑149 (MQ=255)
                                                                                      ggTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTa                                                             <  3:719346/149‑1 (MQ=255)
                                                                                          aaaGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGAc                                                         >  4:26451/1‑149 (MQ=255)
                                                                                              tttttCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCAtt                                                     >  3:1091960/1‑149 (MQ=255)
                                                                                                  tCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGc                                                 <  2:574426/149‑1 (MQ=255)
                                                                                                  tCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGACTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGc                                                 >  3:974969/1‑149 (MQ=255)
                                                                                                    cTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGt                                               >  4:312175/1‑149 (MQ=255)
                                                                                                     tgtgAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGta                                              <  3:243593/149‑1 (MQ=255)
                                                                                                     tgtgAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGta                                              <  1:587739/149‑1 (MQ=255)
                                                                                                     tgtgAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGta                                              >  4:454696/1‑149 (MQ=255)
                                                                                                        gAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATaa                                           >  2:43389/1‑149 (MQ=255)
                                                                                                              gCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTc                                     <  4:980059/149‑1 (MQ=255)
                                                                                                              gCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTc                                     >  3:599493/1‑149 (MQ=255)
                                                                                                               cTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCt                                    <  4:417764/149‑1 (MQ=255)
                                                                                                                   gcgAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCAt                                >  4:81531/1‑149 (MQ=255)
                                                                                                                     gagaGTGTGGCGGAGTAACGAAACCGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATaa                              <  3:454696/149‑1 (MQ=255)
                                                                                                                           gtgGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCa                        <  3:26451/149‑1 (MQ=255)
                                                                                                                               cGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACTCTTCTGCATAATGCCCAGTAc                    <  3:629315/149‑1 (MQ=255)
                                                                                                                                 gAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGt                  >  1:239524/1‑149 (MQ=255)
                                                                                                                                  aGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGtt                 <  2:431710/149‑1 (MQ=255)
                                                                                                                                      aCGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCt             <  3:344951/149‑1 (MQ=255)
                                                                                                                                        gAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGcc                          >  4:866050/1‑134 (MQ=255)
                                                                                                                                        gAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGcc                          <  3:866050/134‑1 (MQ=255)
                                                                                                                                        gAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCa                        <  1:786100/136‑1 (MQ=255)
                                                                                                                                        gAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCa                        >  2:786100/1‑136 (MQ=255)
                                                                                                                                         aaaaCGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGgc          <  2:803221/149‑1 (MQ=255)
                                                                                                                                          aaaCGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGgcg         >  4:739709/1‑149 (MQ=255)
                                                                                                                                           aaCGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGgcgc        >  1:307339/1‑149 (MQ=255)
                                                                                                                                              ggCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAAc     <  3:739709/149‑1 (MQ=255)
                                                                                                                                               gCCCTGCATGCGAATGCCGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACg    >  3:630313/1‑149 (MQ=255)
                                                                                                                                               gCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACg    >  2:781202/1‑149 (MQ=255)
                                                                                                                                                cccTGCATGCGAATGCCGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTAc                                                            >  2:1145815/1‑92 (MQ=39)
                                                                                                                                                cccTGCATGCGAATGCCGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTAc                                                            <  1:1145815/92‑1 (MQ=39)
                                                                                                                                                 ccGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGgc          <  1:801958/141‑1 (MQ=255)
                                                                                                                                                 ccGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGgc          >  2:801958/1‑141 (MQ=255)
                                                                                                                                                 ccGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGtc  >  4:3462/1‑149 (MQ=255)
                                                                                                                                                 ccGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGtc  <  2:307339/149‑1 (MQ=255)
                                                                                                                                                 ccGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGtc  >  4:101123/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAAGAGATCCCTGTCGCAGGCACCCGCGTGCGTATTGGCGAGTACGATATCGATATTCTCGACGTACAGGACAATATGATTAAGCAGGTAAAAGTTTTTCCTGTGAAACCGCTGCGCGAGAGTGTGGCGGAGTAACGAAAACGGCCCGGCATTCGCATGCAGGGCCGTGAATTATTACGAAAGCAGAAATTAAGCTTTTGCTTTCGCTACAGTAACCATCGCCGCACGAATCGTACGACCATTCAGCGTATAACCCTTCTGCATAATGCCCAGTACGTTACCTGGCGCAACGTC  >  NC_000913/2749926‑2750219

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