Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,264 A→C 6.4% intergenic (+34/‑183) typA → / → yihL GTP‑binding protein/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,060,2640AC 6.4% 419.7 / 19.3 126intergenic (+34/‑183)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (55/63);  new base C (5/3);  total (60/66)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.24e-01

Only 100 of 128 total aligned reads displayed.

GTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAA  >  NC_000913/4060116‑4060410
                                                                                                                                                    |                                                                                                                                                    
gTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGa                                                                                                                                                      >  4:69801/1‑149 (MQ=255)
gTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGa                                                                                                                                                      >  3:357211/1‑149 (MQ=255)
   cATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                                                                   >  4:223836/1‑149 (MQ=255)
    aTCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                                  >  4:343800/1‑149 (MQ=255)
    aTCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                                  >  3:327708/1‑149 (MQ=255)
    aTCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                                  >  1:677014/1‑149 (MQ=255)
     tCGATGATAACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                                                                                 >  1:108461/1‑149 (MQ=255)
       gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAgg                                                                                                                                               >  2:658911/1‑149 (MQ=255)
        atgatgACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑t                                                                                                                                             <  1:51617/149‑1 (MQ=255)
          gatggCCAACTTGGAGACGTGCTTCGGACCTCTCTCCGTATTGGTAAACGTCACGTGACGGAAAACGATCGTCGCCGCGCCAACCGCCCACCGAAAGACGATTTATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGC‑ttt                                                                                                                                           <  4:436147/149‑1 (MQ=255)
          gatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAAAAAAAACCTGCCATCGCTGGCAGG‑ttt                                                                                                                                           <  1:626389/149‑1 (MQ=255)
             gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tttttt                                                                                                                                        >  4:149511/1‑149 (MQ=255)
                gAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttttt                                                                                                                                     <  2:549913/149‑1 (MQ=255)
                gAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttttttttt                                                                                                                                     >  4:489471/1‑149 (MQ=255)
                 aaCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttatga                                                                                                                                    <  1:522321/149‑5 (MQ=255)
                 aaCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑g                                                                                                                                   >  2:602941/1‑149 (MQ=255)
                  aCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGCCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑Gc                                                                                                                                  <  3:274955/149‑1 (MQ=255)
                    tGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCt                                                                                                                                >  3:641011/1‑149 (MQ=255)
                       tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTc                                                                                                                             <  4:61556/149‑1 (MQ=255)
                       tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTc                                                                                                                             <  3:229484/149‑1 (MQ=255)
                          aaGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTg                                                                                                                          >  2:489383/1‑149 (MQ=255)
                           aGTGGCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGa                                                                                                                         <  4:635427/149‑1 (MQ=255)
                           aGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGa                                                                                                                         <  1:348159/149‑1 (MQ=255)
                               aCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAgcgc                                                                                                                     >  1:376009/1‑149 (MQ=255)
                                cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGca                                                                                                                    <  2:38061/149‑1 (MQ=255)
                                 tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                                  >  1:430295/1‑120 (MQ=255)
                                 tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                                                                  <  2:430295/120‑1 (MQ=255)
                                 tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGcac                                                                                                                   >  4:170303/1‑149 (MQ=255)
                                 tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGcac                                                                                                                   >  4:573558/1‑149 (MQ=255)
                                   cGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCacac                                                                                                                 >  2:563972/1‑149 (MQ=255)
                                    gACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACAcc                                                                                                                >  2:17461/1‑149 (MQ=255)
                                       ctctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑tttt                                                                                                                                          >  4:355462/1‑121 (MQ=255)
                                       ctctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑tttt                                                                                                                                          <  3:355462/121‑1 (MQ=255)
                                        tctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATg                                                                                                            >  1:361568/1‑149 (MQ=255)
                                         ctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGc                                                                                                           <  2:529454/149‑1 (MQ=255)
                                           aTCCGTATTCGTAAACGTCACCTGACGGAAGACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATCAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGa                                                                                                         <  3:644234/149‑1 (MQ=255)
                                           aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGa                                                                                                         <  4:105394/149‑1 (MQ=255)
                                           aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGa                                                                                                         <  4:533845/149‑1 (MQ=255)
                                                  ttCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                                                                  <  3:262824/149‑1 (MQ=255)
                                                   tCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                                                                  <  3:67768/149‑1 (MQ=255)
                                                     gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                                                                               >  3:649324/1‑149 (MQ=255)
                                                      tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                                                                              <  3:115907/149‑1 (MQ=255)
                                                       aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                                                                             <  3:205656/149‑1 (MQ=255)
                                                              aCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCt                                                                                      <  2:376009/149‑1 (MQ=255)
                                                               ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtg                                                                                     <  1:658911/149‑1 (MQ=255)
                                                                cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgt                                                                                    >  4:20300/1‑149 (MQ=255)
                                                                    cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCATACCAATGCGAGTATCTCCCTTATAGGTCTGTGGtt                                                                                <  2:491248/149‑1 (MQ=255)
                                                                     ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                                                                               >  3:88170/1‑149 (MQ=255)
                                                                     ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                                                                               >  3:75214/1‑149 (MQ=255)
                                                                     ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                                                                               >  4:269871/1‑149 (MQ=255)
                                                                      gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                                                                              >  2:277999/1‑149 (MQ=255)
                                                                      gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                                                                              <  2:389278/149‑1 (MQ=255)
                                                                         aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                                                                        <  4:38045/136‑1 (MQ=255)
                                                                         aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                                                                        >  3:38045/1‑136 (MQ=255)
                                                                          aCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                                                                          >  4:587501/1‑149 (MQ=255)
                                                                            gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                                                                        <  2:427350/149‑1 (MQ=255)
                                                                                gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                                                                    >  3:168475/1‑149 (MQ=255)
                                                                                      gcgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt                                                              <  3:69801/149‑1 (MQ=255)
                                                                                       cgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg                                                             <  2:237899/149‑1 (MQ=255)
                                                                                        gcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt                                                            <  4:327708/149‑1 (MQ=255)
                                                                                        gcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt                                                            <  3:20300/149‑1 (MQ=255)
                                                                                         ccAACCGCGCACCGAAAGAGGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tttttttt                                                                                                                                      <  4:511906/75‑1 (MQ=255)
                                                                                         ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCATG‑tgtttttt                                                                                                                                      >  3:511906/1‑75 (MQ=255)
                                                                                            aCCGCGCACCGAAAGAGATTTAATTTCGTTTTTTATTCATAAAAAACCGGCCGGCGATGG‑AGG‑TTTTTTTTTGGCCGGTCGGGGGGCCACACAAATGCAGGTATCTCCC‑TATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa                                                       <  3:33205/149‑1 (MQ=255)
                                                                                             ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa                                                       <  1:424948/149‑1 (MQ=255)
                                                                                             ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa                                                       >  1:65790/1‑149 (MQ=255)
                                                                                               gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGTCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt                                                     >  4:665274/1‑149 (MQ=255)
                                                                                               gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt                                                     <  4:170897/149‑1 (MQ=255)
                                                                                                 gcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTAt                                                   <  3:267935/149‑1 (MQ=255)
                                                                                                  cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATACAGAGCTTTGTCACAATTaac                                                  <  2:613373/149‑3 (MQ=255)
                                                                                                   aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCt                                                 <  3:573558/149‑1 (MQ=255)
                                                                                                          gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGt                                          >  3:565987/1‑149 (MQ=255)
                                                                                                             gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                                                                           >  2:304803/1‑97 (MQ=255)
                                                                                                             gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCAt                                       <  2:332019/149‑1 (MQ=255)
                                                                                                              aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATa                                      >  2:3651/1‑149 (MQ=255)
                                                                                                              aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATa                                      <  1:470698/149‑1 (MQ=255)
                                                                                                              aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATa                                      <  1:268544/149‑1 (MQ=255)
                                                                                                                         tttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAAtt                           <  4:415678/149‑1 (MQ=255)
                                                                                                                             tAGGCATAAAAAACCTGCCAGGGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCAAACCAATGCGTGTATCTCACTTATAAGTCTGTTGTTACGTCAATAGTGAGCTTTGTCACGATTATCTGCAAAGTCATACGCCGTTAATTGctt                       <  4:244425/149‑1 (MQ=255)
                                                                                                                                tCATAAAAAACCTGCCAGCGATGGCAGG‑‑TTTTATTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAGTTGctttctt                   >  4:583972/1‑149 (MQ=255)
                                                                                                                                  aTAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGctttcttt                  >  4:103063/1‑149 (MQ=255)
                                                                                                                                   tAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttt                 >  3:324022/1‑149 (MQ=255)
                                                                                                                                    aaaaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttt                 <  2:340937/148‑1 (MQ=255)
                                                                                                                                    aaaaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttt                 >  1:340937/1‑148 (MQ=255)
                                                                                                                                     aaaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt               <  4:7728/149‑1 (MQ=255)
                                                                                                                                     aaaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt               <  3:200040/149‑1 (MQ=255)
                                                                                                                                     aaaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCtttttt               <  2:677014/149‑1 (MQ=255)
                                                                                                                                      aaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTg              <  4:168475/149‑1 (MQ=255)
                                                                                                                                       taaCCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATAGCTTACTTTTTTgg             >  3:410873/2‑149 (MQ=255)
                                                                                                                                       aaaCCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGt                                <  2:638238/130‑1 (MQ=255)
                                                                                                                                       aaaCCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGt                                >  1:638238/1‑130 (MQ=255)
                                                                                                                                        aaCCTGCCATCGCTGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgc            >  2:694147/1‑149 (MQ=255)
                                                                                                                                        aaCCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                                                                   <  3:117328/94‑1 (MQ=255)
                                                                                                                                         aCCTGCCATCGCTGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcg           >  3:148645/1‑149 (MQ=255)
                                                                                                                                         aCCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcg           >  2:620354/1‑149 (MQ=255)
                                                                                                                                          ccTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgt          <  3:595900/149‑1 (MQ=255)
                                                                                                                                           cTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgta         >  1:188832/1‑149 (MQ=255)
                                                                                                                                                aGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagcgt    >  3:318750/1‑149 (MQ=255)
                                                                                                                                                 gCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagcgta   >  1:139272/1‑149 (MQ=255)
                                                                                                                                                  cGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGgcgtaagcgtaa  <  2:198098/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTT‑GCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAAAGTCATACGCCGTTAATTGCTTTCTTTTTTGGCGTAAGCGTAA  >  NC_000913/4060116‑4060410

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