Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 336,866 G→C 5.0% G30A (GGC→GCC)  yahF → putative NAD(P)‑binding succinyl‑CoA synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913336,8660GC 5.0% 415.8 / 11.2 120G30A (GGC→GCC) yahFputative NAD(P)‑binding succinyl‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base G (59/55);  new base C (3/3);  total (62/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-01

Only 100 of 120 total aligned reads displayed.

AACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACATCCCTGGGAACGCAACTCAGGCAATAGAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATT  >  NC_000913/336604‑337134
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                            
aaaaacGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTACGACATCCCTGGGAACGCAACTCAGGCAATAGAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAgg                                                                                                                                                                                                                                                                   <  1:75837/271‑1 (MQ=255)
                                          ttatttatGACATCCCTGGGAACGCAACTCAGGCAATAGAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTgaa                                                                                                                                                                                                                        <  2:47577/275‑1 (MQ=255)
                                                atGACATCCCTGGGAACGCAACTCAGGCAATAGAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCCGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGTGTGCCACCTAAATGCAAAAAGGCGTTCTGAAAAAtt                                                                                                                                                                                                                   >  4:368497/1‑274 (MQ=255)
                                                       cccTGGGAACGCAACTCAGGCAATAGAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGAc                                                                                                                                                                                                                                                >  2:168641/1‑238 (MQ=255)
                                                              aaCGCAACTCAGGCAATAGAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGc                                                                                                                                                                                                    >  2:336019/1‑275 (MQ=255)
                                                                                               atTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTAGGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCtga                                                                                                                                                                   <  2:384430/275‑1 (MQ=255)
                                                                                                                        tGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCg                                                                                                                                                                                                                                                                            <  6:54828/145‑1 (MQ=255)
                                                                                                                        tGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCg                                                                                                                                                                                                                                                                            >  5:54828/1‑145 (MQ=255)
                                                                                                                            gTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCa                                                                                                                                                                                                                                                                     <  7:949937/148‑1 (MQ=255)
                                                                                                                                  catggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCAtt                                                                                                                                                                                                                                                               >  5:130600/1‑148 (MQ=255)
                                                                                                                                  catggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGAACTAATGATTGTCCACAATGGTAAATCGGGTGCGGAAAAc                                                                                                                                >  4:83320/1‑275 (MQ=255)
                                                                                                                                    tggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAACTGAATAAAGGCGTGCTGAAGAATTTAGGCCTGCTGACGCCCGAGCTGGAGCAGGCGAAAAACGGCGCCCTGATGATTGTCATCAACGGTAAATAGGGGGCCGGCAACgg                                                                                                                              >  4:120519/1‑274 (MQ=255)
                                                                                                                                     ggcagagtgTTGTGCGCTGAGCCAAAGCCACGGAGGGAATTATGTCAGTTAAAATATGCATTAAAACGCAGACCTAGTTTGAGTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTAGACGGCGTCGAGCAGGAATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAg                                                                                                                             <  4:257715/267‑1 (MQ=255)
                                                                                                                                        cGATGGTTGTGCGCTGAGCCAAACATACGGAGGGAATTATGTCAGTTAAAATAGTCATTAACCCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAAt                                                                                                                                                     <  1:8550/248‑1 (MQ=255)
                                                                                                                                         gATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGt                                                                                                                         <  4:451930/275‑1 (MQ=255)
                                                                                                                                                gtgCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGtggcgatggcga                                                                                                                                                                                                                                                 >  6:1091730/1‑148 (MQ=255)
                                                                                                                                                   cgcTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGga                                                                                                               <  2:368498/275‑1 (MQ=255)
                                                                                                                                                    gcTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGga                                                                                                               <  2:34279/274‑1 (MQ=255)
                                                                                                                                                        aGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTCATGGTGGAGATTg                                                                                                          >  2:310640/1‑275 (MQ=255)
                                                                                                                                                         gCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGa                                                                                                                                                                                                                                        <  7:1100058/148‑1 (MQ=255)
                                                                                                                                                             aaCCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATaa                                                                                                                                                                                                                                    <  5:897268/148‑1 (MQ=255)
                                                                                                                                                             aaCCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTgaagaa                                                                                                     >  1:42232/1‑275 (MQ=255)
                                                                                                                                                                  aCGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTAGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACCGCGACCTGATGATTGTCATCCATGGTACATCGGGTGCGGACCAAGAAAAGTTACTTGTGGACTTTGAAGGACttt                                                                                                 >  4:330794/1‑272 (MQ=255)
                                                                                                                                                                   cGGAGGGAATTATGTCAGTTAAAATAGTCATTAAAACGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTc                                                                                                                                                                                                                                                                          >  8:963596/1‑104 (MQ=255)
                                                                                                                                                                   cGGAGCGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTc                                                                                                                                                                                                                                                                          <  7:963596/104‑1 (MQ=255)
                                                                                                                                                                     gAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGGCGGCGGCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCaa                                                                                             <  1:94295/275‑1 (MQ=255)
                                                                                                                                                                     gAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGCCGTCGAg                                                                                                                                                                                                                                                                       >  8:911366/1‑105 (MQ=255)
                                                                                                                                                                     gAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGCCGTCGAg                                                                                                                                                                                                                                                                       <  7:911366/105‑1 (MQ=255)
                                                                                                                                                                      atggAATTATTTAACTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCaa                                                                                             <  2:330795/272‑1 (MQ=255)
                                                                                                                                                                        ggAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTga                                                                                                                                                                                                                         <  5:1091730/148‑1 (MQ=255)
                                                                                                                                                                        ggAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCa                                                                                                                           >  4:199931/1‑242 (MQ=255)
                                                                                                                                                                         gAATTATGTCAGTTAAAATAGTCATTAAACCGGATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAg                                                                                                                                                                                                                                                                    >  8:539102/1‑104 (MQ=255)
                                                                                                                                                                         gAATTATGTCAGTTAAAATAGTCATTAAACCGGATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAg                                                                                                                                                                                                                                                                    <  7:539102/104‑1 (MQ=255)
                                                                                                                                                                         gAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCAtt                                                                                                                                                                                                                                                               <  6:1001297/109‑1 (MQ=255)
                                                                                                                                                                          aaTTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTTCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACAcc                                                                                         >  2:171677/1‑274 (MQ=255)
                                                                                                                                                                            ttATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTa                                                                                                                       <  2:199932/242‑1 (MQ=255)
                                                                                                                                                                                gTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTa                                                                                                                                                                                                                 >  6:717075/1‑148 (MQ=255)
                                                                                                                                                                                  cAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATaaa                                                                                                                                                                                                                                   >  8:863905/1‑128 (MQ=255)
                                                                                                                                                                                  cAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATaaa                                                                                                                                                                                                                                   <  7:863905/128‑1 (MQ=255)
                                                                                                                                                                                  cAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAgg                                                                                                                                                                                                               <  6:130600/148‑1 (MQ=255)
                                                                                                                                                                                    gTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGCAAAGGAGCAGTTACGGGTGGCGCTTGAAGAAATGTGAAACACcaaagcgcaaa                                                                               >  4:384491/1‑274 (MQ=255)
                                                                                                                                                                                           tAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGccc                                                                        >  2:75236/1‑272 (MQ=255)
                                                                                                                                                                                            aGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCa                                                                                                                                                                                                                                                                 <  8:918919/88‑1 (MQ=255)
                                                                                                                                                                                              tCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGcc                                                                                                                                                                                                   <  5:717075/148‑1 (MQ=255)
                                                                                                                                                                                                 ttAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGaa                                                                                                                                                                                                                                            <  6:793700/104‑1 (MQ=255)
                                                                                                                                                                                                 ttAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGaa                                                                                                                                                                                                                                            >  5:793700/1‑104 (MQ=255)
                                                                                                                                                                                                 ttAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTgaaga                                                                                                                                                                                                                      <  6:19676/126‑1 (MQ=255)
                                                                                                                                                                                                 ttAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTgaaga                                                                                                                                                                                                                      <  8:362449/126‑1 (MQ=255)
                                                                                                                                                                                                 ttAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTgaaga                                                                                                                                                                                                                      >  5:19676/1‑126 (MQ=255)
                                                                                                                                                                                                 ttAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTgaaga                                                                                                                                                                                                                      >  7:362449/1‑126 (MQ=255)
                                                                                                                                                                                                  tAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTACCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATg                                                                                                                                                                                                                                         <  7:23746/106‑1 (MQ=255)
                                                                                                                                                                                                  tAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTACCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATg                                                                                                                                                                                                                                         >  8:23746/1‑106 (MQ=255)
                                                                                                                                                                                                         aaTACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGtggcgatg                                                                                                                                                                                                                                                     >  5:255156/1‑87 (MQ=255)
                                                                                                                                                                                                         aaTACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGtggcgatg                                                                                                                                                                                                                                                     <  6:255156/87‑1 (MQ=255)
                                                                                                                                                                                                         aaTACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGa                                                                                                                                                                                                       >  4:150260/1‑133 (MQ=255)
                                                                                                                                                                                                         aaTACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGATACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGAGCAAATCGGCTCGCACGAGGCGCGTTa                                                         >  1:54940/1‑275 (MQ=255)
                                                                                                                                                                                                          aTACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGaaa                                                                                                                                                                                                                                           >  5:626292/1‑96 (MQ=255)
                                                                                                                                                                                                             ccTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGcc                                                                                                                                                                                                   <  2:150261/133‑1 (MQ=255)
                                                                                                                                                                                                               tATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGa                                                                                                                                                                                  >  5:870163/1‑148 (MQ=255)
                                                                                                                                                                                                                 ttttGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACAAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACtc                                                 >  4:106941/1‑274 (MQ=255)
                                                                                                                                                                                                                  tttGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGa                                                                                                                                                                                                                                        <  8:1094510/91‑1 (MQ=255)
                                                                                                                                                                                                                  tttGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGa                                                                                                                                                                                                                                        >  7:1094510/1‑91 (MQ=255)
                                                                                                                                                                                                                       ttCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGtggcgatggcg                                                                                                                                                                                                                                                  >  7:80639/1‑76 (MQ=255)
                                                                                                                                                                                                                        tCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGt                                                                                                                                                                                                                                                                           >  5:857792/1‑50 (MQ=255)
                                                                                                                                                                                                                        tCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTa                                                                                                                                                                                                                 <  8:198014/108‑1 (MQ=255)
                                                                                                                                                                                                                        tCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTa                                                                                                                                                                                                                 >  7:198014/1‑108 (MQ=255)
                                                                                                                                                                                                                           gTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAgcg                                        >  4:263139/1‑274 (MQ=255)
                                                                                                                                                                                                                             ctcGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAgg                                                                                                                                                                                                                                 >  7:753972/1‑87 (MQ=255)
                                                                                                                                                                                                                             ctcGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAAttt                                                                                                                                                                                                                  >  8:943230/1‑102 (MQ=255)
                                                                                                                                                                                                                             ctcGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAAttt                                                                                                                                                                                                                  <  7:943230/102‑1 (MQ=255)
                                                                                                                                                                                                                             ctcGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCACTTGTGGCGATGGCGACCGAAATGAATAAAgg                                                                                                                                                                                                                                 <  8:753972/87‑1 (MQ=255)
                                                                                                                                                                                                                                   gATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAgg                                                                                                                                                                                                               <  5:235407/99‑1 (MQ=255)
                                                                                                                                                                                                                                        cTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGAcc                                                                                                                                                                      <  8:559878/135‑1 (MQ=255)
                                                                                                                                                                                                                                        cTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGtcat                                                                                                                                                         <  7:989630/148‑1 (MQ=255)
                                                                                                                                                                                                                                         tATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCa                                                                                                                                                                                       <  7:447079/117‑1 (MQ=255)
                                                                                                                                                                                                                                          aTCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGaaa                                                                                                                                                                                                                                           >  6:680150/1‑64 (MQ=255)
                                                                                                                                                                                                                                             tcCACGCGTGAAAATAAACGCGACAGCGGCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGa                     <  4:171676/275‑1 (MQ=255)
                                                                                                                                                                                                                                              ccACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGaa                                                                                                                                                                                                                                            >  7:589936/1‑59 (MQ=255)
                                                                                                                                                                                                                                              ccACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGaa                                                                                                                                                                                                                                            <  8:589936/59‑1 (MQ=255)
                                                                                                                                                                                                                                               cACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAGGGCg                                                                                                                                                                                                                               >  6:493181/1‑71 (MQ=255)
                                                                                                                                                                                                                                                aCGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTg                                                                                                                                                                                                                          <  6:77311/75‑1 (MQ=255)
                                                                                                                                                                                                                                                aCGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTg                                                                                                                                                                                                                          >  5:77311/1‑75 (MQ=255)
                                                                                                                                                                                                                                                 cgcgTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAgg                                                                                                                                                                                                                                 <  5:565666/67‑1 (MQ=255)
                                                                                                                                                                                                                                                 cgcgTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAgg                                                                                                                                                                                                                                 >  6:565666/1‑67 (MQ=255)
                                                                                                                                                                                                                                                 cgcgTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGt                 >  2:57198/1‑275 (MQ=255)
                                                                                                                                                                                                                                                     tGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCCCGAGGCGCGTTACGCCACTATTGGCATCGCCAAAAAGCATATCCCGGAAAGTAAcc             >  2:153675/1‑275 (MQ=255)
                                                                                                                                                                                                                                                     tGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGATACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCCTGTTACGCCAGCATTCGCCACCACAACAAGCATTTCCCGGAAAGTAAcc             >  4:150861/1‑275 (MQ=255)
                                                                                                                                                                                                                                                     tGCAAAAAAACTCGACGGGGTCGGGCAGGAATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCat                             <  4:212876/259‑1 (MQ=255)
                                                                                                                                                                                                                                                         aaTAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAgg                                                                                                                                                                                                               <  6:884625/77‑1 (MQ=255)
                                                                                                                                                                                                                                                          aTAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTgatgat                                                                                                                                                               >  8:524963/1‑124 (MQ=255)
                                                                                                                                                                                                                                                          aTAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTgatgat                                                                                                                                                               <  7:524963/124‑1 (MQ=255)
                                                                                                                                                                                                                                                          aTAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCg        <  2:263140/275‑1 (MQ=255)
                                                                                                                                                                                                                                                              gcTCGACGCCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAgg                                                                                                                                                                                                               <  8:311699/71‑1 (MQ=255)
                                                                                                                                                                                                                                                               cTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCa                                                                                                                                                                                       >  8:428470/1‑95 (MQ=255)
                                                                                                                                                                                                                                                               cTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCa                                                                                                                                                                                       <  7:428470/95‑1 (MQ=255)
                                                                                                                                                                                                                                                                tCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAAt                                                                                                                                                                                                                    >  5:583611/1‑65 (MQ=255)
                                                                                                                                                                                                                                                                tCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAAt                                                                                                                                                                                                                    <  6:583611/65‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cGACGGAGTGGGGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGAtt  <  4:75235/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aCGCCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCggg                                                                                                                                          <  6:633108/136‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aCGCCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCggg                                                                                                                                          >  5:633108/1‑136 (MQ=255)
                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                            
AACAACGCTGCCGGGAAAATATTCATCGTTTTATTCACAACATTATTTATGACATCCCTGGGAACGCAACTCAGGCAATAGAAAAAATTAAACATATTGGCTCCAGTTCCGGCTGCGACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCAGTTACTGGTGGAGATTGAAGAACTGTTCAACACCAAAGCGCAAAGCGGCTCGCACGAGGCGCGTTACGCCACTATTGGCAGCGCCAAAAAGCATATCCCGGAAAGTAACCTGGCGGTGATT  >  NC_000913/336604‑337134

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