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

GCGCTCCACCAGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAACCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGT  >  NC_000913/4144924‑4145211
                                                                                                                                                     |                                                                                                                                           
gcgcTCCACCAGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCg                                                                                                                                           <  2:501054/151‑1 (MQ=255)
  gcTCCACCAGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGaa                                                                                                                                         <  1:591613/151‑1 (MQ=255)
     ccaccaGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCt                                                                                                                                      <  1:614040/151‑1 (MQ=255)
      caccaGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTc                                                                                                                                     >  2:119965/1‑151 (MQ=255)
       accaGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCg                                                                                                                                    <  2:932091/151‑1 (MQ=255)
        ccaGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCgg                                                                                                                                   >  2:868272/1‑151 (MQ=255)
         caGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGc                                                                                                                                  <  2:803123/151‑1 (MQ=255)
             aGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTgc                                                                                                                              <  1:1007578/151‑1 (MQ=255)
              gcgcgcGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTgcg                                                                                                                             <  2:1098936/151‑1 (MQ=255)
                   cgTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCAc                                                                                                                          >  1:549762/1‑149 (MQ=255)
                    gTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACt                                                                                                                         >  1:326928/1‑149 (MQ=255)
                      tATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAAc                                                                                                                      <  1:811558/150‑1 (MQ=255)
                       aTTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACg                                                                                                                     <  2:211160/150‑1 (MQ=255)
                         ttCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGcc                                                                                                                   >  1:222957/1‑150 (MQ=255)
                          tCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCtt                                                                                                                 >  1:1068509/1‑151 (MQ=255)
                           cGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTg                                                                                                                <  1:1077887/151‑1 (MQ=255)
                               ggTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCg                                                                                                            <  2:549083/151‑1 (MQ=255)
                                    cccGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACg                                                                                                       <  1:255560/151‑1 (MQ=255)
                                       gACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGc                                                                                                      >  1:364606/1‑149 (MQ=255)
                                        aCCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGt                                                                                                   <  1:855572/151‑1 (MQ=255)
                                            gCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCt                                                                                               >  2:829624/1‑151 (MQ=255)
                                            gCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCt                                                                                               <  1:35434/151‑1 (MQ=255)
                                              tataCCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTcc                                                                                                >  2:666903/1‑148 (MQ=255)
                                               ataCCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCtgat                                                                                            <  2:132279/151‑1 (MQ=255)
                                                     cACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAAc                                                                                      >  2:85500/1‑151 (MQ=255)
                                                         gCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAAccgcc                                                                                  >  1:25672/1‑151 (MQ=255)
                                                         gCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAAccg                                                                                    <  1:1050235/149‑1 (MQ=255)
                                                           tCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAAccgcc                                                                                  >  1:210272/1‑149 (MQ=255)
                                                            cggcggGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGa                                                                                <  2:392038/150‑1 (MQ=255)
                                                             ggcggGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCATATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGaa                                                                               >  2:757104/1‑150 (MQ=255)
                                                              gcggGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTCCCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAAc                                                                             <  2:855867/151‑1 (MQ=255)
                                                               cggGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAAcc                                                                            >  2:744782/1‑151 (MQ=255)
                                                                gggTTACCCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCTACACGCCGTTCCTGATGAAAACCGCCGAAACCa                                                                           <  2:494658/151‑1 (MQ=255)
                                                                   ttAACGGAGAACGGCCGTAGCGCGGTGAAGGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCATTTCCCGCAATCAAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCAt                                                                          <  1:1053549/149‑1 (MQ=255)
                                                                       ccGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCTTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGAAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCtt                                                                    <  1:301887/151‑1 (MQ=255)
                                                                       ccGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCAACAAAAGCGTGCAATTCCCGCATTCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCtt                                                                    <  1:515682/151‑1 (MQ=255)
                                                                       ccGAGAACGGACGGAGCGCGGTGAACGTGTTTTCGTTCTTTTGCCTTTACGCCTATCAAAGCGTGCAATTCCCGCAATCGATCCTCGGCGGGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCtt                                                                    <  2:1055761/151‑1 (MQ=255)
                                                                          agaACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGCGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCgg                                                                 >  1:844450/1‑151 (MQ=255)
                                                                            aaCGGACGTAGCGCGGTGAACGTGCTTTCGTTCTTTTGCCTTCACGCCTATCAAGGCGTGCAATTACCGCAATCGACCCGCGGCGTGCGCCCTAACGCCTTGGGCGAAACGCGGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGt                                                                <  1:1021785/150‑1 (MQ=255)
                                                                            aaCGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTa                                                               <  1:626562/151‑1 (MQ=255)
                                                                             agggACTTAGCGGGTTGAACGTTCTGTGTTCCTTTTCCCTGGTCGCCAATCAAAGCGTGGAATTCCCGCAATCAACCCTCGGCGTGCGCACTACCGCCTTGATGGCCACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTa                                                               <  1:920190/148‑1 (MQ=255)
                                                                               ggccGTGGCGCGGTGAACGTGCTTTCGTTCTTTTGCTTTGCCGCCTACAAAAGCGTGAAATTCCCGCAATCGAACCGCGCCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGt                                                                <  1:805083/144‑1 (MQ=255)
                                                                                gACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCgg                                                           >  2:959001/1‑151 (MQ=255)
                                                                                    tAGCGCGGTGAACGTGCTTTCTTTCCTTTGCCTTGACGCCTATCAAAGCTTGCAATTACCGCAATCGACCCTCGGCGTGCGCACTAACGCTTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTAt                                                        <  1:197553/150‑1 (MQ=255)
                                                                                       cgcgGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTcc                                                     <  2:284215/150‑1 (MQ=255)
                                                                                        gcgGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACAAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTcc                                                     <  1:109932/149‑1 (MQ=255)
                                                                                        gcgGTGAACGTGCTTTCGTTCCTTTGCCTTGACCCCTATCAAAGCGTGCAATTACCCCAATCGACCCTCGGCGTGCGCCCTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGc                                                   <  2:132804/151‑1 (MQ=255)
                                                                                             gAACGTGCTTTCGTTCCTTTGCCTGGACGCCTATCAAAGCGTGCAATTCCCGCAATCGACCCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCa                                                  <  1:815007/147‑1 (MQ=255)
                                                                                              aaCGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTCCCGCAATCGACCCTCGGCGTGCGCCCTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCt                                             <  1:194223/151‑1 (MQ=255)
                                                                                                cGTGCTTTCGTTCTTTTGCCTTGACGCCTATCAAAGCGTGCAATTCCCGCAACCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCttt                                           <  2:558493/151‑1 (MQ=255)
                                                                                                   gCTTTCGTTCCTTTGCCTGGCCGCCTTTCAAAGCGTGCAATTCCCGCAATCGACCCTCGGCGTGCCCACTACCTTTTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAAc                                        <  1:1069357/151‑1 (MQ=255)
                                                                                                    cTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACg                                       >  1:10481/1‑151 (MQ=255)
                                                                                                      ttCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTCCCGCAATCGAACCTCGGCGTGCCCACTAACCCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGa                                      <  1:237708/150‑1 (MQ=255)
                                                                                                        cGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATg                                    >  1:1088699/1‑150 (MQ=255)
                                                                                                           tCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAg                                 <  2:326931/150‑1 (MQ=255)
                                                                                                            ccTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAgtg                               >  2:658300/1‑151 (MQ=255)
                                                                                                             cTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAATCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAgtgg                              >  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)
                                                                                                                                                     |                                                                                                                                           
GCGCTCCACCAGTAGCGCGCGTTATTTCGCCGGTTTCCCGACCGGCTATACCGCACTGCTCGGCGGGTTAACCGAGAACGGACGTAGCGCGGTGAACGTGCTTTCGTTCCTTTGCCTTGACGCCTATCAAAGCGTGCAATTACCGCAACCGAACCTCGGCGTGCGCACTAACGCCTTGATCGACACGCCGTTCCTGATGAAAACCGCCGAAACCATTCGCTTCGGTACCGGTATTCCGCAAATCTTTAACGATGAAGTGGTGGTGCCAGCGTTCCTCAACCGTGGCGT  >  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: