Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,684,210 C→G 5.3% intergenic (+154/+44) yphH → / ← glyA putative DNA‑binding transcriptional regulator/serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,684,2100CG 5.3% 339.9 / 8.3 95intergenic (+154/+44)yphH/glyAputative DNA‑binding transcriptional regulator/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (48/42);  new base G (3/2);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

GTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTCACACATCCAG  >  NC_000913/2684069‑2684356
                                                                                                                                             |                                                                                                                                                  
gTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTTCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCACAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtc                                                                                                                                             >  3:125029/1‑148 (MQ=255)
   gAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac                                                                                                                                          <  4:125029/149‑1 (MQ=255)
        gCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGa                                                                                                                                     >  1:1374876/1‑149 (MQ=255)
                                   tAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCg                                                                                                          <  2:1374876/149‑1 (MQ=255)
                                                 ccTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAcc                                                                                                                >  1:668053/1‑129 (MQ=255)
                                                 ccTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAcc                                                                                                                <  2:668053/129‑1 (MQ=255)
                                                             cGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAAt                                                                                                                    >  1:245549/1‑113 (MQ=255)
                                                             cGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAAt                                                                                                                    <  2:245549/113‑1 (MQ=255)
                                                                ttATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGc                                                                                                         <  1:799268/121‑1 (MQ=255)
                                                                ttATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGc                                                                                                         >  2:799268/1‑121 (MQ=255)
                                                                ttATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCtt                                                                                                       <  1:1453803/123‑1 (MQ=255)
                                                                ttATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCtt                                                                                                       >  2:1453803/1‑123 (MQ=255)
                                                                ttATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTg                                                                                  >  1:1334480/1‑144 (MQ=255)
                                                                ttATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTg                                                                                  <  2:1334480/144‑1 (MQ=255)
                                                                   tCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACa                                                                                                                          <  3:144147/101‑1 (MQ=255)
                                                                   tCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACa                                                                                                                          >  4:144147/1‑101 (MQ=255)
                                                                    cAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAcc                                                                                                                >  3:1410992/1‑110 (MQ=255)
                                                                    cAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAcc                                                                                                                <  4:1410992/110‑1 (MQ=255)
                                                                    cAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGt                                                                         >  2:306850/1‑149 (MQ=255)
                                                                        ccTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCgag                                                                     <  1:573476/149‑1 (MQ=255)
                                                                         cTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCgaga                                                                    >  4:972275/1‑149 (MQ=255)
                                                                               ggAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTa                                                              >  4:472049/1‑149 (MQ=255)
                                                                               ggAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTa                                                              >  3:208967/1‑149 (MQ=255)
                                                                                   aTTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAcc                                                                                                                >  3:569972/1‑95 (MQ=255)
                                                                                   aTTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAcc                                                                                                                <  4:569972/95‑1 (MQ=255)
                                                                                     tGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGa                                                                                                                            >  2:773017/1‑81 (MQ=255)
                                                                                     tGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGa                                                                                                                            <  1:773017/81‑1 (MQ=255)
                                                                                      gCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGCTAACGTGCGCAGATGTCGAGAACTTTACCTTTGa                                                       <  3:972275/149‑1 (MQ=255)
                                                                                       cAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac                                                                                                                                       <  1:420312/68‑1 (MQ=255)
                                                                                       cAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac                                                                                                                                       >  2:420312/1‑68 (MQ=255)
                                                                                       cAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGAt                                                      <  1:306850/149‑1 (MQ=255)
                                                                                         atatatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGAt                                                      <  4:89499/147‑1 (MQ=255)
                                                                                         atatatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGAt                                                      >  3:89499/1‑147 (MQ=255)
                                                                                           atatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATgcgc                                                  <  3:472049/149‑1 (MQ=255)
                                                                                              ttGAATTTGCAAAGTTTTGTAGGCCGGATAATGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTa                                                                                                >  4:1413881/1‑100 (MQ=255)
                                                                                              ttGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTa                                                                                                <  3:1413881/100‑1 (MQ=255)
                                                                                              ttGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCga                                                                      >  3:779958/1‑126 (MQ=255)
                                                                                              ttGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCga                                                                      <  4:779958/126‑1 (MQ=255)
                                                                                                gAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGAt                                             >  1:475021/1‑149 (MQ=255)
                                                                                                 aaTTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATa                                            >  1:801282/1‑149 (MQ=255)
                                                                                                      gCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGc                                       >  4:974566/1‑149 (MQ=255)
                                                                                                       cAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCgg                                                                                          >  2:815667/1‑97 (MQ=255)
                                                                                                       cAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCgg                                                                                          <  1:815667/97‑1 (MQ=255)
                                                                                                       cAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGAt                                                      <  2:739031/133‑1 (MQ=255)
                                                                                                       cAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGAt                                                      >  1:739031/1‑133 (MQ=255)
                                                                                                        aaaGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACa                                                                                                                          <  2:562271/64‑1 (MQ=255)
                                                                                                        aaaGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACa                                                                                                                          >  1:562271/1‑64 (MQ=255)
                                                                                                        aaaGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATgcgc                                                  >  3:685923/1‑136 (MQ=255)
                                                                                                        aaaGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATgcgc                                                  <  4:685923/136‑1 (MQ=255)
                                                                                                          aGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCtt                                     >  4:1064712/1‑147 (MQ=255)
                                                                                                          aGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCtt                                     <  3:1064712/147‑1 (MQ=255)
                                                                                                          aGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTtca                                   <  4:208967/149‑1 (MQ=255)
                                                                                                               tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCAtt                              <  3:1114442/149‑1 (MQ=255)
                                                                                                                  aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGAt                           <  2:801282/149‑1 (MQ=255)
                                                                                                                    gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGa                                                                            >  1:1248364/1‑98 (MQ=37)
                                                                                                                    gCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGa                                                                            <  2:1248364/98‑1 (MQ=25)
                                                                                                                     ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGAt                                                                           >  3:708436/1‑98 (MQ=255)
                                                                                                                     ccGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGAt                                                                           <  4:708436/98‑1 (MQ=255)
                                                                                                                      cGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTg                       >  4:1365560/1‑149 (MQ=255)
                                                                                                                          tAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCgg                                                                                          <  1:438848/78‑1 (MQ=255)
                                                                                                                          tAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCgg                                                                                          >  2:438848/1‑78 (MQ=255)
                                                                                                                           aaGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAg                  >  3:636826/1‑149 (MQ=255)
                                                                                                                             ggCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGa                                                       <  3:911314/110‑1 (MQ=255)
                                                                                                                             ggCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGa                                                       >  4:911314/1‑110 (MQ=255)
                                                                                                                              gCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCAc               >  3:1023823/1‑149 (MQ=255)
                                                                                                                                gTTCATGCCGGATGCGGCGTGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGt             >  1:977904/1‑149 (MQ=255)
                                                                                                                                gTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCa                   >  3:859413/1‑143 (MQ=255)
                                                                                                                                gTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCa                   <  4:859413/143‑1 (MQ=255)
                                                                                                                                gTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGt             <  4:1023823/149‑1 (MQ=255)
                                                                                                                                gTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATCGATGCTGTCCAGCACGt             >  3:1002153/1‑149 (MQ=255)
                                                                                                                                 ttCATGCCGGATGCGGCGTGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAAc                                          >  1:285833/1‑119 (MQ=255)
                                                                                                                                 ttCATGCCGGATGCGGCGTGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAAc                                          <  2:285833/119‑1 (MQ=255)
                                                                                                                                 ttCATGCCGGATGCGGCGTGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGc                                       <  1:1019007/122‑1 (MQ=255)
                                                                                                                                 ttCATGCCGGATGCGGCGTGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGc                                       >  2:1019007/1‑122 (MQ=255)
                                                                                                                                 ttCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTtcatc                                 >  3:858249/1‑128 (MQ=255)
                                                                                                                                 ttCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTtcatc                                 <  4:858249/128‑1 (MQ=255)
                                                                                                                                    aCGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCAtt                              >  4:645433/1‑128 (MQ=255)
                                                                                                                                    aCGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCAtt                              <  3:645433/128‑1 (MQ=255)
                                                                                                                                    aCGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTcaca         >  4:284537/1‑149 (MQ=255)
                                                                                                                                     cgccgcATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTcacac        >  4:973707/1‑149 (MQ=255)
                                                                                                                                      gccgcATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCAtt                              <  4:129187/126‑1 (MQ=255)
                                                                                                                                      gccgcATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCAtt                              >  3:129187/1‑126 (MQ=255)
                                                                                                                                      gccgcATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGAt                           <  1:268036/129‑1 (MQ=255)
                                                                                                                                      gccgcATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGAt                           >  2:268036/1‑129 (MQ=255)
                                                                                                                                      gccgcATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGt             <  3:848212/143‑1 (MQ=255)
                                                                                                                                      gccgcATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGt             >  4:848212/1‑143 (MQ=255)
                                                                                                                                       ccgcATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCa                <  1:681980/139‑1 (MQ=255)
                                                                                                                                       ccgcATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCa                >  2:681980/1‑139 (MQ=255)
                                                                                                                                        cgcATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCga                                                                      >  1:832697/1‑84 (MQ=255)
                                                                                                                                        cgcATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCga                                                                      <  2:832697/84‑1 (MQ=255)
                                                                                                                                          cATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTcaca         >  1:211901/1‑143 (MQ=255)
                                                                                                                                          cATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTcaca         <  2:211901/143‑1 (MQ=255)
                                                                                                                                           aTCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGaa                                                                   <  4:1258519/84‑1 (MQ=255)
                                                                                                                                           aTCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGaa                                                                   >  3:1258519/1‑84 (MQ=255)
                                                                                                                                           aTCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCACCACGTCACACATccag  >  3:927938/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTCACACATCCAG  >  NC_000913/2684069‑2684356

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