Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,146,023 A→G T677A (ACG→GCG)  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,146,0230AG99.3% 534.4 / NA 142T677A (ACG→GCG) pflDputative glycine radical domain‑containing pyruvate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (74/67);  minor base T (0/1);  total (74/68)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

GGTGGTTGGCGCGACGCCAGACGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTACGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACGTGGTGAACGCCGACACGTTGCGGGAA  >  NC_000913/4145875‑4146166
                                                                                                                                                    |                                                                                                                                               
ggtggtTGGCGCGACGCCAGACGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCg                                                                                                                                               >  3:154078/1‑151 (MQ=255)
ggtggtTGGCGCGACGCCAGACGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCg                                                                                                                                               >  3:592929/1‑151 (MQ=255)
   ggtTGGCGCGACGCCAGACGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGtt                                                                                                                                             <  2:440183/150‑1 (MQ=255)
      tGGCGCGACGCCAGACGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCt                                                                                                                                          <  2:516636/150‑1 (MQ=255)
       ggCGCGACGCCAGACGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTg                                                                                                                                         <  1:1095454/150‑1 (MQ=255)
           tgACGCCAGACGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACg                                                                                                                                     >  1:565127/2‑150 (MQ=255)
             aCGCCAGACGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTg                                                                                                                                   >  1:778926/1‑150 (MQ=255)
             aCGCCAGACGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTg                                                                                                                                   >  2:333093/1‑150 (MQ=255)
              cGCCAGACGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGa                                                                                                                                  >  2:1471365/1‑150 (MQ=255)
                ccAGACGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAAt                                                                                                                               <  2:1219984/151‑1 (MQ=255)
                  aGACGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTc                                                                                                                             <  1:1405389/151‑1 (MQ=255)
                  aGACGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTc                                                                                                                             >  2:999156/1‑151 (MQ=255)
                   gACGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTc                                                                                                                             >  4:208151/1‑150 (MQ=255)
                    aCGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGCAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCAc                                                                                                                           >  1:1487298/1‑151 (MQ=255)
                    aCGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCAc                                                                                                                           >  2:959716/1‑151 (MQ=255)
                          gTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGc                                                                                                                     <  1:931364/151‑1 (MQ=255)
                          gTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGCCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGc                                                                                                                     >  1:1521831/1‑151 (MQ=255)
                           ttttGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCg                                                                                                                    >  2:1067031/1‑151 (MQ=255)
                             ttGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGa                                                                                                                   <  1:806053/150‑1 (MQ=255)
                                ccGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCAGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGAccc                                                                                                                >  1:1737223/1‑150 (MQ=255)
                                 cGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTg                                                                                                              >  3:455771/1‑151 (MQ=255)
                                 cGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTg                                                                                                              <  1:186071/151‑1 (MQ=255)
                                        cAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAgg                                                                                                         >  1:1585229/1‑149 (MQ=255)
                                        cAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTg                                                                                                       >  2:1542142/1‑151 (MQ=255)
                                        cAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTg                                                                                                       <  1:1106573/151‑1 (MQ=255)
                                          gCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGaa                                                                                                     <  3:338603/151‑1 (MQ=255)
                                          gCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGaa                                                                                                     <  3:338621/151‑1 (MQ=255)
                                           cTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGagg                                                                                                          >  4:212184/1‑143 (MQ=255)
                                            tGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGc                                                                                                   >  1:1607120/1‑151 (MQ=255)
                                              gCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAATGGCCAACAGCGGTGCTGAAGTCAGTCAGTACGCTCGAAAACACGCTGCTGTCTAAGGGTTCGTTGCTGAACGTGAAATTCACTCCTGCGACCCTGGAACGTGAAGCa                                                                                                  <  1:592646/150‑1 (MQ=255)
                                              gCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAg                                                                                                 <  2:448461/151‑1 (MQ=255)
                                                   cggcggCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGTTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGAtt                                                                                             <  1:976554/150‑1 (MQ=255)
                                                   cggcggCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGAtt                                                                                             >  1:300340/1‑150 (MQ=255)
                                                     gcggcTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACg                                                                                          >  2:1212131/1‑151 (MQ=255)
                                                     gcggcTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACg                                                                                          <  1:866040/151‑1 (MQ=255)
                                                              cACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGc                                                                                 >  3:592259/1‑151 (MQ=255)
                                                              cACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGc                                                                                 <  2:1451013/151‑1 (MQ=255)
                                                               aCCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGct                                                                                >  3:272466/1‑150 (MQ=255)
                                                               aCCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGcc                                                                                <  3:233292/151‑1 (MQ=255)
                                                                cccATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCg                                                                               >  2:952880/1‑151 (MQ=255)
                                                                  cATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGAc                                                                             >  2:1291184/1‑151 (MQ=255)
                                                                     gCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGActt                                                                           <  1:1277572/150‑1 (MQ=255)
                                                                       tGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGActtctt                                                                        >  2:787845/1‑151 (MQ=255)
                                                                        gggCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGActtctt                                                                        >  2:1309629/1‑150 (MQ=255)
                                                                        gggCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGActtct                                                                         <  1:1189985/149‑1 (MQ=255)
                                                                        gggCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGActtct                                                                         <  1:1257330/149‑1 (MQ=255)
                                                                        gggCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTa                                                                       <  2:796379/151‑1 (MQ=255)
                                                                            cAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACggg                                                                   >  1:351190/1‑151 (MQ=255)
                                                                             aGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGc                                                                  >  2:200902/1‑151 (MQ=255)
                                                                                 cGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGttt                                                              <  1:678000/151‑1 (MQ=255)
                                                                                   cacaAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTAc                                                            >  3:47205/1‑151 (MQ=255)
                                                                                    acaAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTAcc                                                           <  1:607683/151‑1 (MQ=255)
                                                                                    acaAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTAcc                                                           >  4:211890/1‑151 (MQ=255)
                                                                                     caAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTAccc                                                          >  2:705371/1‑151 (MQ=255)
                                                                                      aaGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCa                                                         <  3:73768/151‑1 (MQ=255)
                                                                                      aaGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCa                                                         <  4:117898/151‑1 (MQ=255)
                                                                                        gggCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGc                                                       <  3:278386/151‑1 (MQ=255)
                                                                                          gCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCtt                                                     <  2:78838/151‑1 (MQ=255)
                                                                                          gCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCtt                                                     <  1:828993/151‑1 (MQ=255)
                                                                                           ccAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGGCCCTGGAAGGTGAAGCAAGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCtt                                                     <  1:73315/150‑1 (MQ=255)
                                                                                            cAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTaa                                                   <  4:55006/151‑1 (MQ=255)
                                                                                            cAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTa                                                    >  1:1392181/1‑150 (MQ=255)
                                                                                               cAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGtt                                                >  2:1601912/1‑151 (MQ=255)
                                                                                                aGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTa                                               <  2:134956/151‑1 (MQ=255)
                                                                                                aGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTa                                               <  1:1568737/151‑1 (MQ=255)
                                                                                                   ggTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTacaa                                            <  2:1734282/151‑1 (MQ=255)
                                                                                                   ggTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTaca                                             <  1:1239954/150‑1 (MQ=255)
                                                                                                      aCTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACat                                         <  1:333043/151‑1 (MQ=255)
                                                                                                        tGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACatat                                       >  2:1606189/1‑151 (MQ=255)
                                                                                                         gAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACata                                        <  2:303467/149‑1 (MQ=255)
                                                                                                           agtcagtcAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTc                                     >  2:200915/1‑150 (MQ=255)
                                                                                                           agtcagtcAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTc                                     >  4:480159/1‑150 (MQ=255)
                                                                                                           agtcagtcAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTc                                     >  1:159190/1‑150 (MQ=255)
                                                                                                           agtcagtcAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCa                                    <  1:1242134/151‑1 (MQ=255)
                                                                                                           agtcagtcAGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCa                                    <  1:342466/151‑1 (MQ=255)
                                                                                                            gtcagtcaGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAg                                   <  2:1641873/151‑1 (MQ=255)
                                                                                                            gtcagtcaGTAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAg                                   <  2:1348343/151‑1 (MQ=255)
                                                                                                              cagtcagtAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGt                                  <  2:519693/150‑1 (MQ=255)
                                                                                                               agtcagtAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATAtt                                      >  4:296658/1‑145 (MQ=255)
                                                                                                               agtcagtAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGttt                                <  1:1488483/151‑1 (MQ=255)
                                                                                                               agtcagtAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGgtt                                <  1:228333/151‑1 (MQ=255)
                                                                                                                gtcagtAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAAGATATTCAGTTTa                               <  1:1318209/151‑1 (MQ=255)
                                                                                                                    gtAAGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACgt                           >  2:1464958/1‑151 (MQ=255)
                                                                                                                      aaGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACgtgg                         <  2:1230169/151‑1 (MQ=255)
                                                                                                                      aaGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACgtgg                         >  1:1559603/1‑151 (MQ=255)
                                                                                                                       aGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACgtggt                        >  2:772401/1‑151 (MQ=255)
                                                                                                                       aGCTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACgtgg                         >  2:536479/1‑150 (MQ=255)
                                                                                                                         cTCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACgtggtg                       >  3:399390/1‑150 (MQ=255)
                                                                                                                          tCGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACGTGGTGa                      >  4:293611/1‑150 (MQ=255)
                                                                                                                           cGATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACGTGGTGAAc                    >  1:1080608/1‑151 (MQ=255)
                                                                                                                            gATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACGTGGTGAACg                   >  1:1156863/1‑151 (MQ=255)
                                                                                                                            gATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACGTGGTGAACg                   <  1:808744/151‑1 (MQ=255)
                                                                                                                            gATAACACGCTGCTGTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACGTGGTCAACg                   >  1:551996/1‑151 (MQ=255)
                                                                                                                                     ctgctgTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACGTGGTGAACGCCGACACGt          >  1:640209/1‑151 (MQ=255)
                                                                                                                                     ctgctgTCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACGTGGTGAACGCCGACACGt          >  2:1596539/1‑151 (MQ=255)
                                                                                                                                           tCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACGTGGTGAACGCCGACACGTTGCggg    >  2:351900/1‑151 (MQ=255)
                                                                                                                                           tCTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACGTGGTGAACGCCGACACGTTGCgg     <  1:689899/150‑1 (MQ=255)
                                                                                                                                            cTAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACGTGGTGAACGCCGACACGTTGCggg    >  2:1531458/1‑150 (MQ=255)
                                                                                                                                             tAACGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACGTGGTGAACGCCGACACGTTGCGGGa   >  2:492139/1‑150 (MQ=255)
                                                                                                                                              aaCGGTGCGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACGTGGTGAACGCCGACACGTTGCGGGaa  <  2:565210/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GGTGGTTGGCGCGACGCCAGACGGTCGTTTTGCCGGAGAACAGCTGGCAGACGGCGGCTTGTCACCCATGCTGGGCCAGGACGCACAAGGGCCAACAGCGGTACTGAAGTCAGTCAGTAAGCTCGATAACACGCTGCTGTCTAACGGTACGTTGCTGAACGTGAAATTCACTCCGGCGACCCTGGAAGGTGAAGCAGGATTACGCAAACTGGCCGACTTCTTACGGGCGTTTACCCAGCTTAAGTTACAACATATTCAGTTTAACGTGGTGAACGCCGACACGTTGCGGGAA  >  NC_000913/4145875‑4146166

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