Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,684,202 C→T 6.3% intergenic (+146/+52) 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,2020CT 6.3% 259.6 / 6.6 79intergenic (+146/+52)yphH/glyAputative DNA‑binding transcriptional regulator/serine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (39/35);  new base T (3/2);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

GTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTCACAC  >  NC_000913/2684069‑2684350
                                                                                                                                     |                                                                                                                                                    
gTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTTCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCACAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtc                                                                                                                                       >  3:125029/1‑148 (MQ=255)
   gAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac                                                                                                                                    <  4:125029/149‑1 (MQ=255)
        gCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAc                                                                                                                                                      <  1:158088/126‑1 (MQ=255)
        gCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCAc                                                                                                                                                      >  2:158088/1‑126 (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                                                                                                 >  2:1453803/1‑123 (MQ=255)
                                                                ttATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCtt                                                                                                 <  1:1453803/123‑1 (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                                                        >  3:208967/1‑149 (MQ=255)
                                                                               ggAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTa                                                        >  4:472049/1‑149 (MQ=255)
                                                                                   aTTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAcc                                                                                                          <  4:569972/95‑1 (MQ=255)
                                                                                   aTTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCAcc                                                                                                          >  3:569972/1‑95 (MQ=255)
                                                                                     tGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGa                                                                                                                      <  1:773017/81‑1 (MQ=255)
                                                                                     tGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGa                                                                                                                      >  2:773017/1‑81 (MQ=255)
                                                                                      gCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGCTAACGTGCGCAGATGTCGAGAACTTTACCTTTGa                                                 <  3:972275/149‑1 (MQ=255)
                                                                                       cAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac                                                                                                                                 >  2:420312/1‑68 (MQ=255)
                                                                                       cAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac                                                                                                                                 <  1:420312/68‑1 (MQ=255)
                                                                                       cAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGAt                                                <  1:306850/149‑1 (MQ=255)
                                                                                         atatatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGAt                                                >  3:89499/1‑147 (MQ=255)
                                                                                         atatatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGAt                                                <  4:89499/147‑1 (MQ=255)
                                                                                           atatTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATgcgc                                            <  3:472049/149‑1 (MQ=255)
                                                                                              ttGAATTTGCAAAGTTTTGTAGGCCGGATAATGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTa                                                                                          >  4:1413881/1‑100 (MQ=255)
                                                                                              ttGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTa                                                                                          <  3:1413881/100‑1 (MQ=255)
                                                                                              ttGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCga                                                                <  4:779958/126‑1 (MQ=255)
                                                                                              ttGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCga                                                                >  3:779958/1‑126 (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                                            <  4:685923/136‑1 (MQ=255)
                                                                                                        aaaGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATgcgc                                            >  3:685923/1‑136 (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                                    <  2:285833/119‑1 (MQ=255)
                                                                                                                                 ttCATGCCGGATGCGGCGTGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAAc                                    >  1:285833/1‑119 (MQ=255)
                                                                                                                                 ttCATGCCGGATGCGGCGTGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGc                                 >  2:1019007/1‑122 (MQ=255)
                                                                                                                                 ttCATGCCGGATGCGGCGTGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGc                                 <  1:1019007/122‑1 (MQ=255)
                                                                                                                                 ttCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTtcatc                           <  4:858249/128‑1 (MQ=255)
                                                                                                                                 ttCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTtcatc                           >  3:858249/1‑128 (MQ=255)
                                                                                                                                    aCGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCAtt                        <  3:645433/128‑1 (MQ=255)
                                                                                                                                    aCGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCAtt                        >  4:645433/1‑128 (MQ=255)
                                                                                                                                    aCGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTcaca   >  4:284537/1‑149 (MQ=255)
                                                                                                                                     cgccgcATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTcacac  >  4:973707/1‑149 (MQ=255)
                                                                                                                                     |                                                                                                                                                    
GTAGAAAGGCTTCCCGAAGGAAGCCTTGATGATCATAAACGAAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACACGGAGATTGCAATATATTGAATTTGCAAAGTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGAGCACATTGACAGCAAATCACCGTTTCGCTTATGCGTAAACCGGGTAACGTGCGCAGATGTCGAGAACTTTACCTTTGATGCGCTCGATAACGGCTTCATCATTGATGCTGTCCAGCACGTCACAC  >  NC_000913/2684069‑2684350

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