Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 951,885 T→G 5.6% D557A (GAC→GCC)  pflB ← formate C‑acetyltransferase 1, anaerobic; pyruvate formate‑lyase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913951,8850TG 5.6% 257.5 / 5.1 84D557A (GAC→GCC) pflBformate C‑acetyltransferase 1, anaerobic; pyruvate formate‑lyase 1
Reads supporting (aligned to +/- strand):  ref base T (41/38);  new base G (2/3);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGG  >  NC_000913/951738‑952029
                                                                                                                                                   |                                                                                                                                                
gTCAGAACAGACTGAGTCGGGATAGCGGCACGGTATGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGt                                                                                                                                                  <  4:197024/148‑1 (MQ=255)
 tCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTc                                                                                                                                                 >  2:3370/1‑148 (MQ=255)
 tCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTc                                                                                                                                                 <  1:3370/148‑1 (MQ=255)
   aGAACAGACTGAGTCGGGAAAGCGTCACGGTAGGTGAGCAGTTTCTCAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTAGCCTTCGATTTCGAAGTCGa                                                                                                                                               >  2:325564/1‑148 (MQ=255)
    gAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATa                                                                                                                                             <  8:162854/149‑1 (MQ=255)
     aaCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGt                                                                                                                                                  >  7:195214/1‑143 (MQ=255)
     aaCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGt                                                                                                                                                  <  8:195214/143‑1 (MQ=255)
     aaCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAg                                                                                                                                            >  8:183525/1‑149 (MQ=255)
     aaCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAg                                                                                                                                            >  1:305778/1‑149 (MQ=255)
     aaCAGACTGAGTCGGGATAGCGTCACGGTAGGGGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATACGCCTTCGATTTCGAAGTCGAAAg                                                                                                                                            >  3:89399/1‑149 (MQ=255)
      aCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGc                                                                                                                                           <  2:298614/149‑1 (MQ=255)
           cTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGAc                                                                                                                                      <  4:22326/149‑1 (MQ=255)
            tGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGAcc                                                                                                                                     >  8:161484/1‑149 (MQ=255)
            tGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTATGAATTTTCTTCATGAAACGTTCTACCAGGACAACAGCCAGGTCATCTACACGCGGATAATTGATACCAAACTCCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGAcc                                                                                                                                     >  1:320366/1‑149 (MQ=255)
               gTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTc                                                                                                                                  >  7:277239/1‑149 (MQ=255)
                  gggATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGGCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTc                                                                                                                               <  5:165437/149‑1 (MQ=255)
                      tAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCa                                                                                                                           <  4:66349/149‑1 (MQ=255)
                                              tCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCatat                                                                                                   >  6:47895/1‑149 (MQ=255)
                                              tCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCata                                                                                                    <  1:359413/148‑1 (MQ=255)
                                               cTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCatat                                                                                                   <  7:183525/148‑1 (MQ=255)
                                                tGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATAttt                                                                                                 >  8:11149/1‑149 (MQ=255)
                                                  taTTTTCTTCATGAAACGTTCTACTAGGTCTTCATCCAGGTCATCTAGACGCGTATAATTGTTACCAAATTGCGGGTATTCGCCTTCGATTTCGAAGTCGTTATCCAGATCGTCTTCGTCACGAATCGGTTTAACTTTCGCATAttt                                                                                                 <  2:320366/146‑1 (MQ=255)
                                                  aaTTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATAtt                                                                                                  <  8:289857/146‑1 (MQ=255)
                                                  aaTTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGa                                                                                               <  3:71242/149‑1 (MQ=255)
                                                    ttttCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAACCTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGAtt                                                                                             <  2:98846/149‑1 (MQ=255)
                                                       tcttcATGAAACGTTCTACCAGTTCAACAGCCAGGTAATCTACACGCTGATCATTGTTACCAAACTGCGGGTATTAGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGc                                                                                           <  4:263661/148‑1 (MQ=255)
                                                       tcttcATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATAATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGc                                                                                           <  2:305778/148‑1 (MQ=255)
                                                         ttcATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATAtt                                                                                                  >  7:61005/1‑139 (MQ=255)
                                                         ttcATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATAtt                                                                                                  <  8:61005/139‑1 (MQ=255)
                                                           cATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGAc                                                                                       <  8:277239/148‑1 (MQ=255)
                                                            aTGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGAt                                                                                              >  6:369361/1‑140 (MQ=255)
                                                            aTGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGAt                                                                                              <  5:369361/140‑1 (MQ=255)
                                                            aTGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACa                                                                                      <  7:11149/148‑1 (MQ=255)
                                                            aTGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAg                                                                                     >  1:330998/1‑149 (MQ=255)
                                                             tGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAgg                                                                                    <  4:49691/149‑1 (MQ=255)
                                                              gAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTAGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAgg                                                                                    >  2:47892/1‑148 (MQ=255)
                                                               aaaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGAtt                                                                                             >  6:91675/1‑138 (MQ=255)
                                                               aaaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGAtt                                                                                             <  5:91675/138‑1 (MQ=255)
                                                               aaaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAggg                                                                                   <  7:161484/148‑1 (MQ=255)
                                                                aaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGGAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAg                                                                                 >  1:202949/1‑149 (MQ=255)
                                                                aaCGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAg                                                                                 >  2:7609/1‑149 (MQ=255)
                                                                 aCGTTCTACCCGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCGTCGTCCCGAATCGGTGTAACTTTCGCATATTTGATTGCAGACAGGGAgg                                                                                <  1:325564/149‑2 (MQ=255)
                                                                       tACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTcagcag                                                                          >  2:253868/1‑149 (MQ=255)
                                                                          cAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTcagc                                                                            >  4:222780/1‑144 (MQ=255)
                                                                          cAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTcagc                                                                            <  3:222780/144‑1 (MQ=255)
                                                                          cAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCaa                                                                       >  6:26949/1‑149 (MQ=255)
                                                                          cAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCaa                                                                       >  4:86076/1‑149 (MQ=255)
                                                                              tCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGGAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGa                                                                   <  2:202949/149‑1 (MQ=255)
                                                                              tCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGAGAGGGAGTCAGCAGCAACGGa                                                                   <  1:247262/149‑1 (MQ=255)
                                                                                aaCAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGa                                                                   >  7:238267/1‑147 (MQ=255)
                                                                                  cAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGa                                                                   <  7:191059/145‑1 (MQ=255)
                                                                                  cAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGa                                                                   >  8:191059/1‑145 (MQ=255)
                                                                                   aGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGAc                                                              <  5:26949/149‑1 (MQ=255)
                                                                                       aGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATAtt                                                                                                  <  5:285276/109‑1 (MQ=255)
                                                                                       aGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATAtt                                                                                                  >  6:285276/1‑109 (MQ=255)
                                                                                       aGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                                          >  7:360732/1‑149 (MQ=255)
                                                                                         gTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGa                                                        <  5:47895/149‑1 (MQ=255)
                                                                                               tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCa                                                                                                                           <  1:83900/76‑1 (MQ=255)
                                                                                               tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCa                                                                                                                           >  2:83900/1‑76 (MQ=255)
                                                                                               tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                                          >  8:218280/1‑141 (MQ=255)
                                                                                               tACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGc                                                          <  7:218280/141‑1 (MQ=255)
                                                                                                    gcgGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCATACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACACCACGc                                              >  6:281209/1‑148 (MQ=255)
                                                                                                    gcgGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGc                                              >  5:162288/1‑148 (MQ=255)
                                                                                                           aTTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGc                                      <  8:238267/149‑1 (MQ=255)
                                                                                                               ttACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGa                                   >  4:67351/1‑148 (MQ=255)
                                                                                                                   cAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACg                              >  1:2064/1‑149 (MQ=255)
                                                                                                                   cAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACg                              >  1:8067/1‑149 (MQ=255)
                                                                                                                   cAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACg                              >  6:22733/1‑149 (MQ=255)
                                                                                                                      aCTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCa                           <  1:7609/149‑1 (MQ=255)
                                                                                                                            ggTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCg                     >  6:356354/1‑149 (MQ=255)
                                                                                                                             gTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGt                    >  7:119977/1‑149 (MQ=255)
                                                                                                                               tatCGACTTCGAAATCGAAGGCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTg                   >  4:334552/3‑148 (MQ=255)
                                                                                                                                ttCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGc                  <  8:119977/148‑1 (MQ=255)
                                                                                                                                ttCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGc                  <  5:22733/148‑1 (MQ=255)
                                                                                                                                   gCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAgcg              >  8:215613/1‑149 (MQ=255)
                                                                                                                                    acTTCGAAATCGAAGGCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGc                                                                                           <  8:177809/70‑1 (MQ=255)
                                                                                                                                    acTTCGAAATCGAAGGCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGc                                                                                           >  7:177809/2‑71 (MQ=255)
                                                                                                                                     cTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAgcgc             >  7:263863/1‑148 (MQ=255)
                                                                                                                                         gATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCa        <  8:263863/149‑1 (MQ=255)
                                                                                                                                          aTTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATa                                                      >  6:251983/1‑102 (MQ=255)
                                                                                                                                             tCGAAGGCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCagaga    <  3:334552/149‑1 (MQ=255)
                                                                                                                                              cGAAGTCGATAGCCAGACCGTCTTCGTCACGGATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCagaga    <  5:281209/148‑1 (MQ=255)
                                                                                                                                              cGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGCCAGGGAGTCAGCAGCAACGGACAGCCCAGCGATACCACCCGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCa        >  5:307035/1‑144 (MQ=255)
                                                                                                                                              cGAAGGCGATAGCCAGCCCGTCTGCGGCCCGAATCGGTTTCACTTTCGCATTTTTGATTGCAGACAGGGAGTAAGCAGCACCGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGGCGTGCAGCGCCATCa        <  6:307035/144‑1 (MQ=255)
                                                                                                                                               gAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt          <  6:342076/141‑1 (MQ=255)
                                                                                                                                               gAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAgg  >  3:98367/1‑149 (MQ=255)
                                                                                                                                               gAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGAGCGGATAACGTCACGGTCGTGCAGCGCCATCagctag   >  5:342076/1‑144 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTCATCTACACGCGGATCATTGTTACCAAACTGCGGGTATTCGCCTTCGATTTCGAAGTCGATAGCCAGACCGTCTTCGTCACGAATCGGTTTAACTTTCGCATATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGG  >  NC_000913/951738‑952029

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: