Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,145,072 C→T P360S (CCG→TCG)  pflD → putative glycine radical domain‑containing pyruvate formate‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,145,0720CT98.7% 262.4 / NA 76P360S (CCG→TCG) pflDputative glycine radical domain‑containing pyruvate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (34/41);  total (34/42)
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.00e+00

GCGCTCCACCAGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAACCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGT  >  NC_000913/4144924‑4145211
                                                                                                                                                     |                                                                                                                                           
gcgcTCCACCAGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCg                                                                                                                                           <  2:501054/151‑1 (MQ=255)
  gcTCCACCAGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGaa                                                                                                                                         <  1:591613/151‑1 (MQ=255)
     ccaccaGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCt                                                                                                                                      <  1:614040/151‑1 (MQ=255)
      caccaGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTc                                                                                                                                     >  2:119965/1‑151 (MQ=255)
       accaGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCg                                                                                                                                    <  2:932091/151‑1 (MQ=255)
        ccaGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCgg                                                                                                                                   >  2:868272/1‑151 (MQ=255)
         caGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGc                                                                                                                                  <  2:803123/151‑1 (MQ=255)
             aGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTgc                                                                                                                              <  1:1007578/151‑1 (MQ=255)
              gcgcgcGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTgcg                                                                                                                             <  2:1098936/151‑1 (MQ=255)
                   cgTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCAc                                                                                                                          >  1:549762/1‑149 (MQ=255)
                    gTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACt                                                                                                                         >  1:326928/1‑149 (MQ=255)
                      tATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAAc                                                                                                                      <  1:811558/150‑1 (MQ=255)
                       aTTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACg                                                                                                                     <  2:211160/150‑1 (MQ=255)
                         ttCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGcc                                                                                                                   >  1:222957/1‑150 (MQ=255)
                          tCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCtt                                                                                                                 >  1:1068509/1‑151 (MQ=255)
                           cGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTg                                                                                                                <  1:1077887/151‑1 (MQ=255)
                               ggTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCg                                                                                                            <  2:549083/151‑1 (MQ=255)
                                    cccGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACg                                                                                                       <  1:255560/151‑1 (MQ=255)
                                       gACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGc                                                                                                      >  1:364606/1‑149 (MQ=255)
                                        aCCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGt                                                                                                   <  1:855572/151‑1 (MQ=255)
                                            gCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCt                                                                                               >  2:829624/1‑151 (MQ=255)
                                            gCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCt                                                                                               <  1:35434/151‑1 (MQ=255)
                                              tataCCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTcc                                                                                                >  2:666903/1‑148 (MQ=255)
                                               ataCCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCtgat                                                                                            <  2:132279/151‑1 (MQ=255)
                                                     cACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAAc                                                                                      >  2:85500/1‑151 (MQ=255)
                                                         gCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAAccgcc                                                                                  >  1:25672/1‑151 (MQ=255)
                                                         gCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAAccg                                                                                    <  1:1050235/149‑1 (MQ=255)
                                                           tCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAAccgcc                                                                                  >  1:210272/1‑149 (MQ=255)
                                                            cggcggGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGa                                                                                <  2:392038/150‑1 (MQ=255)
                                                             ggcggGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCATATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGaa                                                                               >  2:757104/1‑150 (MQ=255)
                                                              gcggGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTCCCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAAc                                                                             <  2:855867/151‑1 (MQ=255)
                                                               cggGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAAcc                                                                            >  2:744782/1‑151 (MQ=255)
                                                                gggTTACCCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCTACACGCCGTTCCTGATGAAAACCGCCGAAACCa                                                                           <  2:494658/151‑1 (MQ=255)
                                                                   ttAACGGAGAACGGCCGTAGCGCGGTGAAGGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCATTTCCCGCAATCAAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCAt                                                                          <  1:1053549/149‑1 (MQ=255)
                                                                       ccGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCT‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGAAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCtt                                                                    <  1:301887/151‑1 (MQ=255)
                                                                       ccGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCAACAAAAGCGTGCAATTCCCGCATTCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCtt                                                                    <  1:515682/151‑1 (MQ=255)
                                                                       ccGAGAACGGACGGAGCGCGGTGAACGTGTTTTCGTTCT‑TTTGCCTTTACGCCTATCAAAGCGTGCAATTCCCGCAATCGATCCTCGGCGGGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCtt                                                                    <  2:1055761/151‑1 (MQ=255)
                                                                          agaACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGCGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCgg                                                                 >  1:844450/1‑151 (MQ=255)
                                                                            aaCGGACGTAGCGCGGTGAACGTGCTTTCGTTCT‑TTTGCCTTCACGCCTATCAAGGCGTGCAATTACCGCAATCGACCCGCGGCGTGCGCCCTAACGCCTTGGGCGAAACGCGGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGt                                                                <  1:1021785/150‑1 (MQ=255)
                                                                            aaCGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTa                                                               <  1:626562/151‑1 (MQ=255)
                                                                             agggACTTAGCGGGTTGAACGTTCTGT‑GTTCCTTTTCCCTGGTCGCCAATCAAAGCGTGGAATTCCCGCAATCAACCCTCGGCGTGCGCACTACCGCCTTGATGGCCACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTa                                                               <  1:920190/148‑1 (MQ=255)
                                                                               ggccGTGGCGCGGTGAACGTGCTTTCGTTCT‑TTTGCTTTGCCGCCTACAAAAGCGTGAAATTCCCGCAATCGAACCGCGCCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGt                                                                <  1:805083/144‑1 (MQ=255)
                                                                                gACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCgg                                                           >  2:959001/1‑151 (MQ=255)
                                                                                    tAGCGCGGTGAACGTGCTTTCTTTCC‑TTTGCCTTGACGCCTATCAAAGCTTGCAATTACCGCAATCGACCCTCGGCGTGCGCACTAACGCTTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTAt                                                        <  1:197553/150‑1 (MQ=255)
                                                                                       cgcgGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTcc                                                     <  2:284215/150‑1 (MQ=255)
                                                                                        gcgGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACAAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTcc                                                     <  1:109932/149‑1 (MQ=255)
                                                                                        gcgGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACCCCTATCAAAGCGTGCAATTACCCCAATCGACCCTCGGCGTGCGCCCTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGc                                                   <  2:132804/151‑1 (MQ=255)
                                                                                             gAACGTGCTTTCGTTCC‑TTTGCCTGGACGCCTATCAAAGCGTGCAATTCCCGCAATCGACCCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCa                                                  <  1:815007/147‑1 (MQ=255)
                                                                                              aaCGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTCCCGCAATCGACCCTCGGCGTGCGCCCTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCt                                             <  1:194223/151‑1 (MQ=255)
                                                                                                cGTGCTTTCGTTCT‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTCCCGCAACCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCttt                                           <  2:558493/151‑1 (MQ=255)
                                                                                                   gCTTTCGTTCC‑TTTGCCTGGCCGCCTTTCAAAGCGTGCAATTCCCGCAATCGACCCTCGGCGTGCCCACTACCTTTTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAAc                                        <  1:1069357/151‑1 (MQ=255)
                                                                                                    cTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACg                                       >  1:10481/1‑151 (MQ=255)
                                                                                                      ttCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTCCCGCAATCGAACCTCGGCGTGCCCACTAACCCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGa                                      <  1:237708/150‑1 (MQ=255)
                                                                                                        cGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATg                                    >  1:1088699/1‑150 (MQ=255)
                                                                                                           tCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAg                                 <  2:326931/150‑1 (MQ=255)
                                                                                                            cc‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAgtg                               >  2:658300/1‑151 (MQ=255)
                                                                                                             c‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAgtgg                              >  1:885221/1‑151 (MQ=255)
                                                                                                                   ccTTGACGCCTATCAAAGCGTGCAATTCCCGCAATCGACCCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAgtggtggtg                         <  1:364548/151‑1 (MQ=255)
                                                                                                                     ttGCCGCCTATCAAAGCGTGCAATTCCCGCAATCGAACCCGGGCGTGCGCCCTACCGCCTGGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAgtggtggtg                         <  2:301007/149‑1 (MQ=255)
                                                                                                                     ttGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGcc                       >  1:1104443/1‑151 (MQ=255)
                                                                                                                     ttGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGcc                       >  1:1199258/1‑151 (MQ=255)
                                                                                                                      tGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGcc                       >  2:832966/1‑150 (MQ=255)
                                                                                                                         cGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGc                    >  1:200939/1‑150 (MQ=255)
                                                                                                                         cGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCg                   >  2:414615/1‑151 (MQ=255)
                                                                                                                           ccTATCAAAGCGTGCAATTCCCGCAATCGACCCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGt                  <  1:866909/150‑1 (MQ=255)
                                                                                                                           ccTATCAAAGCGTGCAATTACCGCAATCGACCCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGtt                 <  2:44491/151‑1 (MQ=255)
                                                                                                                                cAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTc             >  2:754818/1‑150 (MQ=255)
                                                                                                                                 aaaGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTcccc             >  1:51295/1‑147 (MQ=255)
                                                                                                                                 aaaGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTc             >  1:117588/1‑149 (MQ=255)
                                                                                                                                 aaaGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCa            >  2:964420/1‑150 (MQ=255)
                                                                                                                                 aaaGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCa            <  2:25672/150‑1 (MQ=255)
                                                                                                                                  aaGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAAc          >  2:323814/1‑151 (MQ=255)
                                                                                                                                   aGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAAcc         >  1:1012940/1‑151 (MQ=255)
                                                                                                                                     cGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGt       >  2:331975/1‑151 (MQ=255)
                                                                                                                                          aaTTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGt  >  2:851078/1‑151 (MQ=255)
                                                                                                                                          aaTTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGt  >  1:785901/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
GCGCTCCACCAGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCC‑TTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAACCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGT  >  NC_000913/4144924‑4145211

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